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Abstract: We describe a framework to develop, implement and validate any pertur- 
bative Lagrangian-based particle physics model for further theoretical, phenomenological 
and experimental studies. The starting point is FeynRules, a Mathematica package 
that allows to generate Feynman rules for any Lagrangian and then, through dedicated 
interfaces, automatically pass the corresponding relevant information to any supported 
Monte Carlo event generator. We prove the power, robustness and flexibility of this ap- 
proach by presenting a few examples of new physics models (the Hidden Abelian Higgs 
Model, the general Two-Higgs-Doublet Model, the most general Minimal Supersymmet- 
ric Standard Model, the Minimal Higgsless Model, Universal and Large Extra Dimen- 
sions, and QCD-inspired effective Lagrangians) and their implementation/validation in 
FeynArts/FormCalc, CalcHep, MadGraph/MadEvent, and Sherpa. 



Keywords: 



New Physics, Phenomenology, Monte Carlo simulations 



Contents 



i 



Introduction 

A simple example: from the Standard Model to the Hidden Abelian 
Higgs Model 

2.1 The model 



2.2 From model building to phenomenology 



2.3 Validation of New Physics Models 



2.4 From phenomenology to experiment 



|3|. FeynRules in a nutshell 
|3.1| Model description 
Ol Running FeynRules 



i 



Interfaces 

FeynArts/FormCalc 
|4.1.1| Model framework 
|4.1.2| FeynRules interface 
CalcHep/CompHep 
|4.2.1| Model framework 
|4.2.2| FeynRules interface 
MadGraph/MadEvent 
[4.3.1| Model framework 
4.3.21 FeynRules interface 



4.4 



Sherpa 

4.4.1 Model framework 
4.4.2| FeynRules interface 



Models 

?T| The Standard Model 
|5.1.1| Model descript ion 
|5.1.2| FeynRules implementation 
The general Two-Higgs-Doublet Model 
Model description 
FeynRules implementation 



5.2 



5.2.1 



5.2.2 



5.3 



The most general Minimal Supersymmetric Standard Model 
Model description 
FeynRules implementation 



5.3.1 



5.3.2 



5.4 



The Minimal Higgsless Model 
Model description 
FeynRules implementation 



5.4.1 



10 

n 

12 

II 

15 

16 
|16 
17 
18 
19 
19 
19 

m 

m 



29 

m 



5.4.2 



m 

42 
44 
45 
49 



- 1 - 



5.5 



Extra Dimensional models 

5.5.1 Model description 

5.5.2| FeynRules implementation 



5.6 An access to the low energy world 
5.6.1 Model description 
5.6.2| FeynRules implementation 



6.1 



6.2 



6.3 



Validation of the implemented models 

Strategy for the validation of the implemented models 
The Standard Model 
The general Two-Higgs-Doublet Model 

The most general Minimal Supersymmetric Standard Model 
The Minimal Higgsless Model 
Extra Dimensional Models 
Low-energy effective theories 



6.4 



6.5 



6.6 



6.7 



0. Outlook 

|A|. The FeynRules convention for Standard Model inputs 



Validation tables 

The Standard Model 



B.l 



B.2 The general Two-Higgs-doublet model 



B.3| The most general Minimal Supersymmetric Standard Model 



B.4 The Minimal Higgsless Model 



B.5 Extra Dimensional models 



49 

m 



59 

|9 

m 



64 
67 
70 



1. Introduction 

At the Large Hadron Collider (LHC) discoveries most probably will not be an easy task. 
The typical final states produced at this proton-proton collider running at very high energies 
will involve a large number of jets, heavy-flavor quarks, as well as leptons and missing 
energy, providing an overwhelming background to many new physics searches. Complex 
signal final state signatures will then need a very careful understanding of the detector 
and an accurate modeling of the data themselves. In this process, Monte Carlo (MC) 
simulations will play a key role in describing control data sets and devising robust search 
strategies. 

Already the first step, i.e., establishing "an excess over the Standard Model (SM) 
background", might be very difficult, depending on the type of signature involved At 
this stage, matrix-element-based MC (which give reliable predictions for shapes and can still 
be tuned to some extent to the data) will be used to describe backgrounds and possibly 
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candidates signals. For some specific signals, an accurate prediction of the background 
normalization and shapes, validated via control samples, could be also needed. At the 
same time, accurate measurements and comparisons with the best theoretical predictions 
{e.g., at the next-to-next-to-leading order, resummation calculations, ...) of a set of 
standard-candle observables will also be mandatory to claim a good understanding and 
control of physics and detector effects. Very accurate predictions, possibly including even 
weak corrections, and a reliable estimate of errors (such as those introduced by the parton 
distribution functions) will then be needed. 

Once the presence of excess(es) is confirmed, model building activities will be triggered, 
following both top-down and bottom-up approaches. In each case, tools that are able to 
make predictions for wide classes of Beyond the Standard Model (BSM) physics, as well 
as those that help in building up an effective field theory from the data (such as the so 
called OSET method p), could be employed. Finally, as a theoretically consistent picture 
arises, measurements of the parameters (masses, spin, charges) will be carried out. In this 
case it will be necessary to have at least next-to- leading order (NLO) predictions {i.e., a 
reliable normalization) for the signal processes. As our knowledge about the detector and 
the newly discovered physics scenario gets stronger, more accurate determinations will be 
possible and sophisticated analyses tools could be employed, on the very same lines as 
current top quark analyses at the Tevatron collider, e.g., see Ref. 

As schematically outlined above, Monte Carlo simulations will play a key, though 
different role at each stage of the exploration of the TeV scale, i.e., the discovery and 
identification of BSM Physics, and the measurement of its properties. The realization of 
the need for better simulation tools for the LHC has spurred an intense activity over the 
last years, that has resulted in several important advances in the field. 

At the matrix-element level, these include the development of general purpose event 
generators, such as CompHep/CalcHep @, |, |§, MadGraph/MadEvent 0, |, |9|, |lO| , 
Sherpa |11, 12] and Whizard [l^, as well as high efficiency multiparton generators which 
go beyond the usual Feynman diagram techniques, such as Alpgen [|1^], Helac [^] and 
Comix [^]. As a result, the problem of generating automatically leading-order matrix 
elements (and then cross sections and events) for a very large class of renormalizable 
processes has been solved. Quite amazingly, enormous progress has also been achieved 
very recently in the automatization of NLO computations. First the generation of the real 



corrections with the appropriate subtractions has be achieved in an automatic way [17, 18 



19, M, piT. Then several new algorithms for calculating loop amplitudes numerically have 



been proposed (see, e.g., Ref. [^] for a review) and some of them successfully applied to 



the computation of SM processes of physical interest 24, Pq ]. 

An accurate simulation of a hadronic collision requires a careful integration of the 
matrix-element hard process, with the full parton showering and hadronization infrastruc- 



ture, as efficiently provided by Pythia p6|, p^], Herwig p^, p|] and Sherpa. Here also. 



significant progress has been made in the development of matching algorithms such as 



that by Catani, Krauss, Kuhn and Webber (CKKW) fsg, [S^, ||], Mangano (MLM) H 



and others |34, 36 1, in their comparison ||37| , pq ] and application to SM ||3^, ^ and 
BSM [EO] processes. A breakthrough in merging fixed order calculations and parton show- 
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ers was achieved by Frixione, Webber and Nason |41, 42 1, who showed how to correctly 
interface an NLO computation with a parton shower to avoid double counting and deliv- 
ered the first event generator at NLO, MC@NLO. More recently, a new method along the 
same lines, dubbed POWHEG, has been proposed |Q and applied to Drell-Yan and Higgs 
production Q |4|, . 

The progress in the field of Monte Carlo tools outlined above shows that we are, or will 
be soon, able to simulate all the relevant SM processes at the LHC with an unprecedented 
level of accuracy. It is therefore worth considering the status of the predictions for physics 
Beyond the Standard Model. Quite interestingly, the challenges in this case are of a quite 
different nature. The main reason is that presently there is not a leading candidate for 
BSM, but instead a plethora of models have been suggested, based on very different ideas 
in continuous evolution. The implementation of complex BSM models in existing general 
purpose event generators like those enumerated above remains a long, often painstaking 
and error-prone, process. The derivation of the numerous Feynman rules to describe the 
new interactions, and their implementation in codes following conventions is a very uninter- 
esting and time consuming activity. In addition, the validation of a given implementation 
often relies on a comparison of the obtained analytical and numerical results with those 
available in the literature. Again, due to presence of various conventions, the restricted 
number of public results and the lack of a dedicated framework, such a comparison is often 
done manually, in a partial and not systematic way. Finally, besides a handful of offi- 
cially endorsed and publicly distributed BSM models {e.g., the Minimal Supersymmetric 
Standard Model), many implementations remain private or only used by a restricted set 
of theorists and/or phenomenologists, and never get integrated into the official chain of 
simulation tools used by experimental collaborations. Instead, various "home-made" or 
"hacked" versions of existing MC softwares are commonly used for specific BSM process 
studies, leading to problems in the validation, traceability and maintenance procedures. 

In this work we address the problem of having an efficient framework where any new 
physics model can be developed and its phenomenology can be tested against data. A 
first step in the direction of deriving Feynman rules automatically starting from a model 
Lagrangian has been made in the context of the CompHep/CalcHep event generator 
with the LanHep package [48|. Our aim is to go beyond this scheme and create a general 
and flexible environment where communication between theorists and experimentalists in 
both directions is fast and robust. The desiderata for the new physics phenomenological 
framework linking theory to experiment and vice versa which we provide a solution for are: 



1. General and flexible environment, where any perturbative Lagrangian-based model 
can be developed and implemented. 

2. Modular structure with interfaces to several multi-purpose MC's and computational 
tools. 

3. Robust, easy-to-validate and easy-to-maintain. 

4. Integrable in the experimental software frameworks. 
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5. Full traceability of event samples. 

6. Both top-down and bottom-up approaches are natural. 

This paper is organized in five main sections and various appendices. In the first sec- 
tion, by discussing a simple example, we expose the strategy which we propose to address 
the challenges that model builders, phenomenologists and experimentalists have to face to 
study the phenomenology of a new physics model. In Section ^, we briefly recall how the 
FeynRules package works and present some of the new features recently implemented. 
Section ^ contains a brief description of the various interfaces already available. Section ^ 
contains the information about the models that have already been implemented. In Sec- 
tion ^ we present our strategy to validate BSM model implementations, and illustrate our 
procedure on the already implemented models. Finally, in Section |7[ we discuss the out- 
look of our work. In the appendices we collect some technical information as well as a few 
representative validation tables, which constitute the quantitative results of this paper. 



2. A simple example: from the Standard Model to the Hidden Abehan 
Higgs Model 

From the phenomenological point of view, we can distinguish two classes of BSM mod- 
els. The first class of models consists of straightforward extensions of the SM, obtained 
by adding one (or more) new particles and interactions to the SM Lagrangian. In this 
bottom-up approach, one generally starts from the SM, and adds a set of new operators 
according to some (new) symmetry. The second class of models are obtained in a top-down 
approach, where the fundamental Lagrangian is determined by the underlying global and 
local symmetries, and the SM is only recovered in some specific limit. 

In this section we describe in detail our framework to develop, test, implement and 
validate any perturbative Lagrangian-based particle physics model. We concentrate on 
the case of the bottom-up models, and show how our framework allows to easily extend 
the SM and to go in a straightforward way from the model building to the study of the 
phenomenology. We will comment on the top-down models in subsequent sections. 

2.1 The model 



As an illustration, we use the Hidden Abelian Higgs (HAH) Model, described in Ref. |49|. 
This model can be seen as the simplest way to consistently add a new gauge interaction 
to the SM. More specifically, we consider an SU{3)c x SU{2)l x U{1)y x U{l)x gauge 
theory where all SM particles are singlets under the new gauge group U{l)x- A new Higgs 
field (f) is added that is also a singlet under the SM gauge group and breaks the U{l)x 
symmetry when it acquires its vacuum expectation value (vev), (0) = ^/\/2 . The most 
general Lagrangian describing this model is given by 

-ChAH = -^^HAH.Gaugo + -^^HAH.Fermions + -^HAH.Higgs + -^^HAH.Yukawa, (2-1) 
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with 

/:HAH,Gauge = -\g% " \ W^, " \ B,, B^^^ - \ X,, X^^ + | B,, X^^ 

2 (2.2) 



where ^ denotes the SM Higgs field. The covariant derivative reads 



(2.4) 



B^, = d^,-igsT''Gl-i9''--W^, - ig' Y - igx X X^, (2.3) 
and we define the field strength tensors, 

X^lu =d^Xy — dpX^, 
B^u =d^By — dyB^, 

=d^Wl - d,W; + g' e'^'' Wl , 
Gl,=d^gl-d,gl+gs r'^dUl- 

gs, g, g' and gx denote the four coupling constants associated with the SU{3)c x SU{2)l x 
U{1)y X U{l)x gauge groups while T, a*, Y and X are the corresponding generators and 
represents the totally antisymmetric tensor. We do not write explicitly the terms in 
the Lagrangian describing the matter sector of the theory as they are identical to those of 



the SM described in detail in Section 5.1 



-^^HAH.Fermions = -^^SM.Fermions and ^HAH.Yukawa = -^^SM.Yukawa- (2.5) 

The kinetic mixing term in >CHAH,Gauge induces a mixing between the two U{1) gauge fields 
and thus a coupling between the matter fermions and the new gauge boson. The kinetic 
terms for the U{1) fields can be diagonalized via a GL(2,R) rotation, 



xA^ ( \/i^o\ X, 

bJ [ -X I \B, 



(2.6) 



After this field redefinition, the gauge sector takes the diagonal form 

/:HAH,Gauge = -\ G% G^^ - \ W},, " \ B,, B^^ - \ X,, X^\ (2.7) 



but the covariant derivative now contains an additional term coupling the field X^ to the 
hypercharge, 

D^ = d^-igsT^Gl-ig''--W^ - ig' Y B^ - i{gx X + g' r^Y) X^, (2.8) 
with rj = — X^- When the Higgs fields acquire their vev^. 



""^We work in unitary gauge. 
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the gauge symmetry gets broken to SU{3)c x U{1)em, and we obtain a non-diagonal mass 
matrix for the neutral weak gauge bosons, which can be diagonalized by an 0(3) rotation, 



(2.10) 



sm 



( B\ 








\- 










are given by 


Gw — 
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and 










-2su,r/ 



(2.11) 



with = Mj^/Ml^, where Mx and M^q denote the masses of the gauge bosons before 
the kinetic mixing. 



M'x = e qx 9x and = {g' + g") v'/A, 



(2.12) 



and qx denotes the U{\)x charge carried by (j). The photon remains massless while the 
two other states acquire a mass given by 



M: 



Z.Z' 



:i + + Az) ± V(l-4^2-Az)2 + 44ry2 



(2.13) 



As a result of electroweak symmetry breaking, non-diagonal mass terms for the Higgs 
fields appear that can be diagonalized via an orthogonal transformation. 




Ch ■5/1 ) [ ^1 
-Sh Ch I I /l2 



(2.14) 



where the mixing angles and mass eigenvalues are given by 



tan 29h 



M 



1,2 



— A ' 



(2.15) 



Once the Lagrangian is written down and diagonalized in terms of mass eigenstates, 
one can easily identify the minimal set of parameters which the model depends on. Not 
all the parameters introduced above are independent, as most of them are related by some 
algebraic relations, e.g., the relation between the mass eigenvalues of the gauge bosons, 
Eq. ( 2.13| ), and the fundamental parameters appearing in the Lagrangian. Our choice 
of independent input parameters is given in Table ||. All other parameters appearing in 
-Chah, Gauge and £HAH,Higgs Can be reexpressed in terms of these parameters. Note, however, 
that there are strong experimental constraints from LEP on the masses and the couplings 
of additional neutral gauge bosons to fermions, which need to be taken into account when 
building the model. In Ref. [49| it was pointed out that r/ = 0.01 is still allowed. In general. 



in order to determine a benchmark point that takes into account the direct and indirect 
experimental constraints, it is required to perform (loop) computations for several physical 
observables. We will comment more on this in the next subsection. 
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JT d,i dllifc: Lei 


Oy lilUUl 


Value 


U[l)x coupling constant 


ax 


1/127.9 


Kinetic mixing parameter 


V 


U.Ui 


Z pole mass 


Mz 


91.188 GeV 


Z' pole mass 


Mz> 


400.0 GeV 


SM Higgs quartic coupling 


X 


0.42568 


Abelian Higgs quartic coupling 


P 


0.010142 


Abelian/SM Higgs interaction 


K 


0.0977392 



Table 1: Input parameters for the Hidden Abelian Higgs model. Other SM input parameters are 
not shown. 

Let us note that, although Eq. ( |2.2| ) is a very simple extension of the SM, from a more 
technical point of view an implementation of the HAH model in a matrix-element gener- 
ator is already not trivial. In this case, it is not sufficient to start from the existing SM 
implementation and just add the vertices contained in >CHAH,Higgs) because mixing in the 
gauge and scalar sectors implies that all SM vertices involving a Higgs boson and/or a Z 
boson need to be modified. For example, although there is no direct coupling between the 
abelian Higgs field (p and the matter fermions, all the Yukawa couplings receive contribu- 
tions from the two mass eigenstates hi and /i2, weighted by the mixing angle 9^, resulting 
in an almost complete rewriting of the SM implementation. In the next subsection we will 
describe how this difficulty can be easily overcome and the phenomenology of the Hidden 
Abelian Higgs model studied. 

2.2 From model building to phenomenology 

The starting point of our approach is FeynRules (see Section ^). Since in this case we 
are interested in a simple extension of the SM, it is very easy to start from the Feyn- 
Rules implementation of the SM which is included in the distribution of the package and 
to extend the model file by including the new particles and parameters, as well as the HAH 
model Lagrangian of Eq. (|2.2| ). Note that, at variance with the direct implementation into 
a matrix-element generator where we need to implement the vertices one at the time, we 
can work in FeynRules with a Lagrangian written in terms of the gauge eigenstates and 
only perform the rotation to the mass eigenbasis as a second step. This implies that it is 
not necessary to modify jCHAH.Fermions since the new fields only enter through >CHAH,Gauge 
and ^HAH.Higgs- 

Several functions are included in FeynRules to perform sanity checks on the La- 
grangian {e.g., hermiticity) . The diagonalization of the mass matrices can be easily per- 
formed directly in Mathematica(R)^ and FeynRules allows us to easily obtain the Feyn- 
man rules for the model. As already mentioned, not only the Feynman rules of the Higgs 
and gauge sectors are modified with respect to the Standard Model, but also the interaction 
vertices in the fermionic sector change due to the mixing of the scalars and the neutral 

^Mathematica is a registered trademark of Wolfram Research, Inc. 
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Branching ratios for hi 


hi bb 


87.7% 


hi ^ cc 


8.1% 


hi r"*" T~ 


4.2% 


Branching ratios for /12 


/12 hi hi 


14.2% 


/i2 ^ bb 


0.1% 


h2^ ZZ 


26.4% 


/i2 ^ W+ W- 


59.3% 



Table 2: Dominant decay channels 



Branching ratios for Z' 


Z'^jj 


44.1% 


Z' hi hi 


0.9% 


Z ^ /l2 /l2 


0.570 


Z' tt 


12.5% 


Z' bb 


5.2% 


Z' ^ r+ r- 


14.9% 


Z' £+ £- 


29.7% 


Z' ^vv 


9.5% 


z' w+ w- 


2.0% 



the new particles in the HAH model. 



weak bosons. The vertices obtained in this way can already be used for pen & paper work 
during the model building, and to compute simple decay rates and cross sections. Since 
FeynRules stores the vertices in Mathematica, it is easy to use them directly for such 
computations. 

After this preliminary study of our model where the mass spectrum of the theory was 
obtained and basic sanity checks have been performed, typically the model is confronted 
with all relevant direct and indirect constraints coming from experiment. This is a necessary 
step to find areas of parameters space which are still viable. Once interesting regions in 
parameter space are identified, the study of the collider phenomenology of the model, e.g.^ 
at the LHC, can start with the calculation of cross sections and decay branching ratios. Let 
us consider first the the calculation of the decay widths of both SM and new particles. Using 
the FeynArts implementation of the new model obtained via the FeynRules interface, 
it becomes a trivial exercise to compute analytically all tree-level two-body decays for the 
Higgs bosons and the Z' boson (alternatively, one could calculate them numerically via 
e.g., CalcHep or MadGraph/MadEvent). The results for the branching ratios of the 
dominant decay modes are shown, for the benchmark scenario considered here, in Table |2|. 
Once decay widths are known, cross sections can be calculated. However, in many cases 
it is insufficient to have only predictions for total cross sections, as a study of differential 
distributions, with possibly complicated multi-particle final states, is necessary to dig the 
signal out of the backgrounds. Furthermore, even a parton-level description of the events 
might be too simplified and additional radiation coming from the colored initial and final- 
state particles, as well as effects coming from hadronization and underlying events need to 
be accounted for. For this reason, phenomenological studies are in general performed using 
generators which include (or are interfaced to) parton shower and hadronization Monte 
Carlo codes. The parton level events for the hard scattering can be generated by the 
general purpose matrix-element (ME) program, and those events are then passed on to the 
parton shower codes evolving the parton-level events into physical hadronic final states. 
However, similar to FeynArts/FormCalc, for new models the ME generators require the 
form of the new vertices, and different programs use different conventions for the vertices. 
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making it difficult to export the implementation from one ME generator to another. To 
solve this issue, FeynRules includes interfaces to several ME generators that allow to 
output the interaction vertices obtained by FeynRules directly in a format that can be 
read by the external codes. For the moment, such interfaces exist for CalcHep/CompHep, 
MadGraph/MadEvent and Sherpa. It should be emphasized that some of these codes 
have the Lorentz and/or color structures hardcoded, something that limits the range of 
models that can be handled by a given MC. In this respect (and others) each of MC tools 
has its own strengths and weaknesses: having several possibilities available maximizes the 
chances that at least one generator is able to efficiently deal with a given model and the 
case in which several MC tools can be used, as most of the examples discussed in this 
paper, allows for a detailed comparison and robust validation of the implementation. 

For the sake of illustration, we used the MadGraph/MadEvent implementation of 
the HAH model to generate signal events for the 55 — > /i2 — *■ hihi 7766 signal proposed 
in Ref. [^ ] as a signature of this model at the LHC. Using the same set of cuts, and the 
same smearing method as in Ref. [50|, we have been able to easily generate final state 
invariant mass distributions for both signal and background events. The result can be seen 
in Fig. 11], which compares well to the Fig. 5 of Ref. [ 50 1 . 



0.040 



Invariant mass distribution 




300 350 400 

M,,,, (GeV) 



500 



Figure 1: Invariant mass distribution for the four particle final state ^jbb, both for the gg 
/i2 hihi jjbb signal (plain) and the main SM backgrounds (dashed) events at the LHC. All 
simulation parameters and analysis cuts are identical to those listed in Ref. |5(|]. 



2.3 Validation of New Physics Models 

In the previous subsection we discussed how the implementation of the model in Feyn- 
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Rules allows to go all the way from model building to phenomenology without having to 
deal with the technicalities of the various ME generators. In this section we argue that our 
approach does not only allow to exploit the strength of each ME generator, but it also has 
a new power in the validation of BSM models by directly comparing various ME genera- 
tors among themselves. Since the various ME generators use different conventions for the 
interaction vertices it becomes hence possible to compare the implementations of the same 
model into various different matrix-element generators, and thus the different tools, among 
themselves. Furthermore, in many cases generator specific implementations of BSM models 
already exist, at least for restricted classes of models, in which case the FeynRules model 
can be directly validated against the existing tool A, and then exported to any other tool 
B for which an interface exists. In the same spirit, any BSM model should be able to 
reproduce the SM results for observables which are independent of the new physics. 

Let us illustrate this procedure through the example of the HAH model. We start by 
implementing our model into CalcHep, MadGraph/MadEvent and Sherpa by means 
of the corresponding interfaces. We then compute the cross sections for all interesting two- 
to-two processes for this model, and check that we obtain the same numbers from every 
ME generator. Note that, since in this case we have modified the scalar sector of the SM, 
we pay particular attention to the unitarity cancellations inherent to the SM in weak boson 
scattering, showing in this way that our implementation does not spoil these cancellations. 
Since the different codes used for the computation of the cross-sections all rely on different 
conventions for the interaction vertices, we hence demonstrated that our model is consistent 
not only by checking that the cancellations indeed take place in all the implementations, 
but we also get the same results for the total cross section after the cancellation, a strong 
check very rarely performed for a general BSM model implementation. We will comment 
on the validation of more general models in subsequent sections. 

Finally, let us comment on the fact that a robust implementation of a BSM model 
into a code does not only require the model to be validated to a very high-level, but the 
implementation should also be clearly documented in order to assure its portability and 
reproducibility. For this reason FeynRules includes an output to T^jX which allows to 
output the content of the FeynRules model files in a human-readable format, including 
the particle content and the parameters which define the model, as well as the Lagrangian 
and the Feynman rules. 

2.4 Prom phenomenology to experiment 

Our approach does not only apply to phenomenological studies at the theory level, but 
it allows in principle to continue and pass the model to an experimental collaboration 
for full experimental studies. In general, the experimental softwares used to simulate the 
detector effects have very strict requirements regarding sanity checks for new modules {e.g., 
private Monte Carlo programs) to be included in the software, and a very long and tedious 
validation procedure is needed. However, if we will come to the point where we have 
to discriminate between various competing models at the LHC, this approach becomes 
extremely inefficient due to the large number of tools to be validated. In our approach this 
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validation can be avoided, streamlining in this way the whole chain all the way from model 
building to the experimental studies and vice versa. 

Let us again illustrate this statement through the example of the HAH model. Since 
this model is now implemented in FeynRules, we can easily pass it into various ME 
generators using the translation interfaces, and we demonstrated the validation power 
inherent to this approach in the previous section. These models can then be used in a ME 
generator in the same way as any other built-in model, without any modification to the 
original code, i.e., without creating a private version of the ME generator. If the model is 
validated it can easily be passed on to the experimental community, which can then read 
the FeynRules output and use it in their favorite ME generator already embedded in 
their software framework. 

By following this procedure tedious validations for each model implementation in a 
given MC are avoided. In addition, the portability and the reproducibility of the exper- 
imentally tested models is guaranteed. As at the origin of the chain is solely the Feyn- 
Rules model file, all the information is concentrated in one place, and thus everybody can 
reproduce all the results at any stage, from the model building to the collider signatures, 
by starting from the very same file. In addition, since the FeynRules model file contains 
the Lagrangian of the model, it is very easy to go back to the model file and understand 
its physical content, a step which might be very difficult working with manually created 
model files for the ME generators, written in an often rather cryptic programming language 
hiding the essential physics. In this way it becomes very easy to use later the very same 
model information, and to reproduce analyses by just changing benchmark points or by 
adding a single new particle/interaction to the same model. 



3. FeynRules in a nutshell 

FeynRules is a Mathematica package which allows to derive Feynman rules directly 



from a Lagrangian [51|. The user provides the Lagrangian for his/her model (written in 
Mathematica) as well as all the information about the particle and parameter content of 
the model. This information uniquely defines the model, and hence is enough to derive all 
the interaction vertices from the Lagrangian. FeynRules can in principle be used with 
any model which fulfills basic quantum field theoretical requirements {e.g., Lorentz and 
gauge invariance), the only current limitation coming from the kinds of fields supported 
by FeynRules (see below). In particular it can also be used to obtain Feynman rules for 
effective theories involving higher-dimensional operators. In a second step, the interaction 
vertices obtained by FeynRules can be exported by the user to various matrix-element 
generators by means of a set of translation interfaces included in the package. In this 
way the user can directly obtain an implementation of his/her model into these various 
tools, making it straightforward to go from model building to phenomenology. Presently, in- 
terfaces to CalcHep/CompHep, FeynArts/FormCalc, MadGraph/MadEvent and 
Sherpa are available. In the following we briefly describe the basic features of the package 
and the model files, the interfaces to the matrix-element generators being described in 
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S Scalar fields 

F Fermion fields (Dirac and Major ana) 

V Vector fields 

U Ghost fields 

T Spin two fields 

Table 3: Particle classes supported by FeynRules. 



Section ^. For more details on both the FeynRules package as well as the interfaces, we 
refer the reader to the FeynRules manual and to the FeynRules website [H, |52|. 

3.1 Model description 

The FeynRules model definition is an extension of the FeynArts model file format and 
consists of the definitions of the particles, parameters and gauge groups that characterize 
the model and the Lagrangian. This information can be placed in a text file or in a 
Mathematica notebook or a combination of the two as convenient for the user. 

Let us start with the particle definitions. Following the original FeynArts conven- 
tion, particles are grouped into classes describing "multiplets" having the same quantum 
numbers, but possibly different masses. Each particle class is defined in terms of a set 
of class properties, given as a Mathematica replacement list. For example, the up-type 
quarks could be written as 

F[l] == { ClassName -> q, 

ClassMembers -> {u, c, t}, 

Self Conjugate -> False, 

Indices -> {Index [Generation] , Index [Colour] } , 

Flavor Index -> Generation, 

Mass -> {Mq, 0, 0, {MT, 174. 3», 

Width -> {Wq, 0, 0, {WT, 1.508», 

QuantumNumbers -> {Q -> 2/3}, 

PDG -> {2, 4, 6» 

This defines a Dirac fermion (F) represented by the symbol q. Note that the antiparticles are 
automatically declared and represented by the symbol qbar. This class has three members, 
u, c, and t (ubar, cbar, and tbar for the antiparticles, respectively), distinguished by a 
generation index (whose range is defined at the beginning of the model definition) . These 
fields carry an additional index labelled Colour. The complete set of allowed particle classes 
is given in Table ^ Additional information, like the mass and width of the particles, as 
well as the U{1) quantum numbers carried by the fields can also be included. Finally, some 
more specific information not directly needed by FeynRules but required by some of the 
matrix-element generators {e.g., the Particle Data Group (PDG) codes |^^) can also be 
defined. A complete description of the particle classes and properties can be found in the 
FeynRules manual. 
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A Lagrangian is not only defined by its particle content, but also by the local and 
global symmetries defining the model. FeynRules allows to define gauge group classes 
in a way similar to the particle classes. As an example, the definition of the QCD gauge 
group can be written 

SU3C == { Abelian -> False, 

GaugeBoson -> G, 

CouplingConstant -> gs, 
StructureConstant -> f , 
Representations -> {T, Colour}} 

where the gluon field G is defined together with the quark field during the particle dec- 
laration. The declaration of abelian gauge groups is analogous. FeynRules uses this 
information to construct the covariant derivative and field strength tensor which the user 
can use in his/her Lagrangian. 

The third main ingredient to define a model is the set of parameters which it depends 
on. In general, not all the parameters appearing in a Lagrangian are independent, but 
they are related through certain algebraic relations specific to each model, e.g., the re- 
lation cosOw = Mw/Mz relating at tree-level the masses of the weak gauge bosons to 
the electroweak mixing angle. FeynRules therefore distinguishes between external and 
internal parameters. External parameters denote independent parameters which are given 
as numerical inputs to the model. An example of a declaration of an external parameter 
reads 

aS == {ParameterType -> External, 

Value -> 0.118} 

defining an external parameter aS with numerical value 0.118. Several other properties 
representing additional information needed by matrix-element generators are also available, 
and wc refer to the FeynRules manual for an extensive list of parameter class properties. 
Internal parameters are defined in a similar way, except that the Value is given by an 
algebraic expression linking the parameter to other external and/or internal parameters. 
For example, the cos 6yj parameter definition could read 

cw == -[ParameterType -> Internal, 

Value -> MW/MZ} 

Note that it is also possible to define tensors as parameters in exactly the same way, as 
described in more detail in the manual. 

At this point, we need to make a comment about the conventions used for the different 
particle and parameter names inside FeynRules. In principle, the user is free the choose 
the names for the gauge groups, particles and parameters at his/her convenience, without 
any restriction. The matrix-element generators however have certain information hard- 
coded {e.g., reference to the strong coupling constant or electroweak input parameters). 
For this reason, conventions regarding the implementation of certain SM parameters has 
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been established to ensure the proper translation to the matrix-element generator. These 
are detailed in the manual and recalled in Appendix 

In complicated models with large parameter spaces, it is sometimes preferable to re- 
strict the model to certain slices of that parameter space. This can be done as usual by 
adjusting the parameters to lie at that particular point in parameter space and doing the 
desired calculation. However, sometimes these slices of parameter space are special in that 
many vertices become identically zero and including them in the Feynman diagram calcu- 
lation can be very inefficient. In order to allow both the general parameter space and a 
restricted parameter space, we introduce the model restriction. A model restriction is a 
Mathematica list containing replacements to the parameters which simplify the model. 
For example, in the SM the CKM matrix has non-zero matrix elements, but it is sometimes 
useful to restrict a calculation to a purely diagonal CKM matrix. Rather than creating a 
new implementation of the SM with a diagonal CKM matrix, a restriction can be created 
and used when desired. The following statement restricts the SM to a diagonal CKM 
matrix 

CKM[i_, i_] :> 1, 

CKM[i_, j_] :> /; i != j. 

When this restriction is applied, all vertices containing off diagonal CKM elements vanish 
identically and are removed before passing it on to a matrix-element generator. The result is 
that these vertices never appear in Feynman diagrams and the calculation is more efficient. 
Several restriction files can be created corresponding to various different slices of parameter 
space. The user then selects the restriction file that they are interested in and applies it to 
the model before running the translation interfaces. 

3.2 Running FeynRules 

After having loaded the FeynRules package into Mathematica, the user can load the 
model and the model restrictions via the commands 

LoadModelE < model file 1 > , < model file 2 > , ... ]; 
LoadRestriction [ < restriction file > ] ; 

where the model can be implemented in as many files as convenient or it can be imple- 
mented directly in the Mathematica notebook in which case the list of files would be 
empty. The restriction definitions can also be placed in a file or directly in the Mathe- 
matica notebook. The Lagrangian can now be entered directly into the notebook^ using 
standard Mathematica commands, augmented by some special symbols representing spe- 
cific objects like Dirac matrices. As an example, we show the QCD Lagrangian, 

L = -1/4 FS[G, mu, nu, a] FS[G, mu, nu, a] +1 qbar . Ga[mu] . DC[q, mu] ; 

^Alternatively, the Lagrangian can also be included in the model file, in which case it is directly loaded 
together with the model file. 
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where FS[G, mu, nu, a] and DC[q, mu] denote the SU{3)c held strength tensors and 
covariant derivatives automaticahy defined by FeynRules. At this stage, the user can 
perform a set of basic checks on the Lagrangian (hermiticity, normahzation of kinetic 
terms, . . . ), or directly proceed to the derivation of the Feynman rules via the command 

verts = FeynmanRules [ L ]; 

FeynRules then computes all the interaction vertices associated with the Lagrangian L 
and stores them in the variable verts. The vertices can be used for further computations 
within Mathematica, or they can be exported to one of the various matrix-element gen- 
erators for further phenomenological studies of the model. The translation interfaces can 
be directly called from within the notebook, e.g., for the FeynArts interface, 

WriteFeynArtsOutput [ L ] ; 

This will produce a file formatted for use in FeynArts. All other interfaces are called in 
a similar way. As already mentioned, let us note that, even if FeynRules is not restricted 
and can derive Feynman rules for any Lagrangian, the matrix-element generators usually 
have some information on the Lorentz and color structures hardcoded, and therefore they 
are much more limited in the set of vertices they can handle. Each matrix-element generator 
has its own strengths, and in the FeynRules approach the same model can be easily 
exported to various codes, exploiting in this way the strength of each individual tool. In 
practice, the interfaces check whether all the vertices are compliant with the structures 
supported by the corresponding matrix-element generator. If not, a warning is printed and 
the vertex is discarded. Each interface produces at the end a (set of) text file(s), often 
consistently organized in a single directory, which can be read into the matrix-element 
generator at runtime and allows to use the new model in a way similar to all other built-in 
models. For more details on the various interfaces, we refer to Section ^ and to the manual. 

4. Interfaces 

In this section we provide a concise description of the FeynRules interfaces to several 
matrix-element generators and symbolic tools available to perform simulations/calculations 
from Lagrangian-based theories. The most important features of the general structure 
of the new physics models in the codes and their limitations are emphasized. Complete 
description of the options and more technical details can be found in the FeynRules user's 
manual, available on the FeynRules website. Interfaces to other codes, once available, 
will be included in the main release of the package and documented in the user's manual. 

4.1 FeynArts/FormCalc 

FeynArts is a Mathematica package for generating, computing and visualizing Feynman 
diagrams, both at tree-level and beyond [^]. For a given process in a specific model, 
FeynArts starts by generating all the possible topologies, taking into account the number 
of external legs and internal loops associated to the considered case, together with the set of 
constraints required by the user, such as the exclusion of one-particle reducible topologies. 
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This stage is purely topological and does not require any physical input. Based on a 
pre-defined library containing topologies without any external leg for tree-level, one-loop, 
two-loop and three-loop calculations, the algorithm successively adds the desired number 
of external legs. Then, the particles present in the model must be distributed over the 
obtained topologies in such a way that the resulting diagrams contain the external fields 
corresponding to the considered process and only vertices allowed by the model. Finally, 
a Mathematica expression for the sum over all Feynman diagrams is created. 

The second step in the perturbative calculation consists in the evaluation of the am- 
plitudes generated by FeynArts. This can be handled with the help of the Mathemat- 
ica package FormCalc which simplifies the symbolic expressions previously obtained in 
such a way that the output can be directly used in a numerical program [^]. FormCalc 
first prepares an input file that is read by the program Form which performs most of the 



calculation |5£, 57, 5£, 60|. The Lorentz indices are contracted, the fermion traces are 
evaluated, the color structures are simplified using the SU (3) algebra, and the tensor reduc- 
tion is performed^. The results are expressed in a compact form through abbreviations and 
then read back into Mathematica where they can be used for further processing. This 
allows to combine the speed of Form with the powerful instruction set of Mathematica. 

4.1.1 Model framework 

The FeynArts models have a very simple structure which can be easily extended to 
include BSM models. In particular, the current distribution of FeynArts contains al- 
ready several models, including a complete implementation of the Standard Model, as well 
as a Two-Higgs-Doublet Model and a completely generic implementation of the Minimal 
Supersymmetric Standard Model. 

The FeynArts models are separated into two files: 

• The generic model file: This file is not specific to any model, but it contains the 
expressions for the propagators and the Lorentz structures of the vertices for generic 
scalar, fermion and vector fields. Note that since this file is not specific to any model, 
different BSM models can be related to the same generic model file. 

• The classes model file: This file is dedicated to a specific model, and contains the 
declarations of the particles and the analytic expressions of the couplings between the 
different fields. This information is stored in the two lists M$ClassesDescription 
for the particle declarations and M$ClassesCouplings for the couplings. 

FeynArts requires all the particles to be grouped into classes, and as a consequence also 
all the classes couplings must be given at the level of the particle classes. If this is the case, 
the number of Feynman diagrams generated at runtime is much smaller, which speeds up 
the code. Since the FeynRules model files are an extension of the FeynArts classes 
model files, the explicit structure of the particle class definitions is very similar to the 
FeynRules particle classes discussed in Section |^. 



*Let us note that, in order to function correctly, FormCalc requires the ampUtude to be given in 
Feynman gauge. 
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4.1.2 FeynRules interface 

FeynRules includes an interface that allows to output the interaction vertices derived 
from the Lagrangian as a FeynArts model file. Note however that at the present stage 
only the classes model file is generated by FeynRules, the generic model file is hardcoded. 
The generic model file used by FeynRules generated models,f eynrules .gen, is included 
in the FeynRules distribution in the Interf aces/FeynArts subdirectory and needs to be 
copied once and for all into the Models directory of FeynArts. f eynrules . gen is based 
on the corresponding Lorentz.gen file included in FeynArts, with some extensions to 
higher-dimensional scalar couplings as those appearing in non-linear sigma models (see 
Section ^^ ). 

The FeynRules interface to FeynArts can be called within a Mathematica note- 
book via the command 

WriteFeynArtsDutput [ L ] ; 

FeynRules then computes all the vertices associated with the Lagrangian L, and checks 
whether they all have Lorentz structures compatible with the generic couplings in the 
generic coupling file, and if so, it extracts the corresponding classes coupling. It not, 
a message is printed on the screen and the vertex is discarded. At this point we should 
emphasize that in order to obtain FeynArts couplings at the level of the particle classes, it 
is necessary that the Lagrangian is also given completely in terms of particle classes. If the 
interface encounters a Lagrangian term which violates this rule, it stops and redefines all the 
classes such that all particles live in their own class. It then starts over and recomputes all 
the interaction vertices, this time for a Lagrangian where all particle classes are expanded 
out explicitly. In this way a consistent FeynArts model file is obtained which can be used 
with FeynArts. It should however be noted that the generation of the diagrams can be 
considerably slower in this case, which makes it desirable to write the Lagrangian in terms 
of particle classes whenever possible. 

The model file produced by FeynRules has the usual FeynArts structure. Besides 
the lists which contain the definitions of the particle classes and the couplings, the Feyn- 
Rules generated model files contain some more information, which can be useful at various 
stages during the computation: 

- M$ClassesDescription: This is in general a copy of the corresponding list in the 
original FeynRules model file. 

- M$ClassesCouplings: Each entry in the list represents a given interaction between 
the particle classes, together with the associated coupling constant, represented by 
an alias gcxx, xx being an integer. Note that at this level FeynRules does not 
compute the counterterms necessary for loop calculations, but they should be added 
by the user by hand. 

- M$FACouplings: A replacement list, containing the definition of the couplings gcxx 
in terms of the parameters of the model. 
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Furthermore, several other replacement lists (M$ExtParams, M$IntParams, M$Masses) are 
included, containing the values of the parameters of the model, as well as the masses and 
widths of all the particles. 

4.2 CalcHep/CompHep 

The CalcHep |^, |6| and CompHep [Q, |] software automate the tree-level Feynman dia- 
gram calculation and production of partonic level collider events. Models with very general 
Lorentz structures are allowed and general color structures can be incorporated via auxil- 
iary fields. Vertices with more than four particles are not supported at this time. In this 
subsection, we will describe the model file structure and how the FeynRules interface to 
CalcHep and CompHep works. 

4.2.1 Model framework 

Models in CalcHep and CompHep are essentially comprised of four files: 

• prtclsN .mdl: a list of all the particles in the model along with information about 
the particles that is necessary for calculation of Feynman diagrams. 

• varsN.mdl: a list of the independent (external) parameters in the model along with 
their numerical value. 

• funcN.mdl: a list of the dependent (internal) parameters of the model along with 
their functional definition. These definitions can contain any standard mathematical 
functions defined in the C code. 

• IgrngN.mdl: a list of all the vertices in the model. It includes the specification of 
the particles involved in the vertex, an overall constant to multiply the vertex with 
and the Lorentz form of the vertex. 

Note that the letter N in the names of the files is an integer which refers to the number of 
the model. 

4.2.2 FeynRules interface 

The CalcHep/CompHep interface can be invoked with the command 

WriteCHOutput [ L ] 

where L denotes the Lagrangian. 

When invoked, this interface will create the directory M$ModelName with -CH appended 
if it does not already exist. It will create the files prtclsN .mdl, varsN.mdl, funcN.mdl 
and IgrngN.mdl. It will then derive the Feynman rules with four particles or less and write 
them to IgrngN.mdl. It will simplify the vertex list by renaming the vertex couplings as 
xl, x2, x3, etc. and write the definitions of these couplings in funcN.mdl along with the 
other internal parameters. 
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Although CalcHep and CompHep can calculate diagrams in both Feynman and uni- 
tary gauge, they are much faster in Feynman gauge and it is highly recommended to imple- 
ment a new model in Feynman gauge. However, if a user decides to implement the model in 
unitary gauge, he/she should remember that according to the way CalcHep/CompHep 
were written, the ghosts of the massless non-abelian gauge bosons must still be imple- 
mented. (In particular, the gluonic ghosts must be implemented in either gauge for this 
interface.) 

One major constraint of the CalcHep/CompHep system is that the color structure 
is implicit. For many vertices {e.g., quark-quark-gluon) , this is not a problem. However, 
for more complicated vertices, there may be an ambiguity. For this reason, the writers of 
CalcHep/CompHep chose to split them up using auxiliary fields. Although this can be 
done for very general vertices, it is not yet fully automatized in FeynRules. Currently, 
only the gluon four-point vertex and squark-squark-gluon-gluon vertices are automatically 
split up in this way. 

The model files are ready to be used and can be directly copied to the CalcHep/ 
CompHep models directories. The default format for this interface is the CalcHep format. 
A user can direct this interface to write the files in the CompHep format by use of the 
CompHep option. The user who writes CompHep model files should note one subtlety. 
If the model is written to the CompHep directory and if the user edits the model inside 
CompHep and tries to save it, CompHep will complain about any C math library functions 
in the model. Nevertheless, it does understand them. We have checked that if the model 
works in CalcHep, it will work in CompHep and give the same results. 

CalcHep has the ability to calculate the widths of the particles on the fly. By default, 
this interface will write model files configured for automatic widths. This can be turned 
off by setting the option CHAutoWidths to False. This option is set to False if CompHep 
is set to True. 

This interface also contains a set of functions that read and write the external param- 
eters from and to the CalcHep variable files (varsN.mdl). After loading the model into 
FeynRules, the external parameters can be updated by running 

ReadCHExtVars [ Input -> < file > ] 

This function accepts all the options of the CalcHep interface plus the option Input 
which instructs FeynRules where to find the CalcHep variable file. The default is 
varsN.mdl in the current working directory. If reading a CompHep variable file, then the 
option CompHep should be set to true. After reading in the values of the variables in the 
CalcHep file, it will update the values in FeynRules accordingly. 

The current values of the external parameters in FeynRules can also be written to a 
CalcHep external variable file (varsN.mdl) using 

WriteCHExtVars [ Output -> < file >] 

This can be done to bypass writing out the entire model if only the model parameters are 
changed. 
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4.3 MadGraph/MadEvent 

The MadGraph/MadEvent v4.4 software [0, |8|, |, |l^ allows users to generate tree-level 
amplitudes and parton-level events for any process (with up to nine external particles). It 



uses the Helas library |61, 62] to calculate matrix elements using the helicity formalism 
in the unitary gauge. Starting from version 4, users have the possibility to use several 
pre-defined BSM models, including the most generic Two-Higgs-Doublet Model and the 
Minimal Supersymmetric Standard Model, but can also take advantage of the UsrMod 
interface to implement simple Standard Model extensions. 

The existing scheme for new model implementations in MadGraph/MadEvent has 
two major drawbacks. First, users need to explicitly provide algebraic expressions for the 
coupling values used by MadGraph to calculate amplitudes. Second, the first version of 
the UsrMod interface only works for models extending the existing Standard Model by 
adding a limited set of new particles and/or interactions. This renders difficult any attempt 
to modify existing BSM models, or to generalize models previously implemented with this 
method. 

The current version of Mad Graph rely on a new clearly defined structure for all 
model libraries. New model libraries can be generated via the corresponding interface 
to FeynRules, which generates all the required code files automatically. Finally, a new 
version of the UsrMod scripts exists which can be used complementary to FeynRules 
for simple extensions of existing models. All these three new frameworks are introduced 
and described in the present section. 

4.3.1 Model framework 

All model libraries supported in the latest versions of MadGraph/MadEvent now have 
the same structure. They are composed of a set of text and Fortran files grouped in a 
single directory, stored in the Models subdirectory of the root MadGraph/MadEvent 
installation: 

• particles.dat: a text file containing a list of all particles entering the model and 
the corresponding properties (name, spin, mass, width, color representation, PDG 
code, . . . ) 

• parain_card . dat: a text file containing the numerical values of the necessary external 
parameters for a specific model. The parameter card has a format compliant with 
the SUSY Les Houches Accord (LHA) and is dependent on the physics model. One 
should pay attention to the fact that some of these parameters are related one to each 
other {e.g., the masses and the widths are generally related to more fundamental La- 
grangian parameters). If possible, this file should also contain (preferably at the end) 
a list of Les Houches QNUMBERS blocks describing properties of non-SM particles 
to facilitate the interface of matrix-element and parton-shower based generators, as 



proposed in Ref. [33|. 



intparam_def inition. inc: a text file containing all the algebraic expressions re- 
lating internal parameters to external and/or internal parameters. There are two 
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different kinds of internal parameters. Indeed, most of the expressions can be com- 
puted once and for all, but in some cases where the parameter depends on the scale of 
the process {e.g., the strong coupling constant), it might be desirable to re-evaluate 
it at an event-by-event basis. 

• interactions . dat: a text file containing a list of all interactions entering the model. 
Each interaction is characterized by an ordered list of the involved particles, the name 
of the corresponding coupling, the corresponding type of coupling (for coupling order 
restrictions) and possible additional switches to select particular Helas routines. 

• couplingsXX . f (where XX can be any integer number): these files contain the alge- 
braic expressions for the couplings, expressed as Fortran formulas. By convention, 
the file couplingsl .f contains all expressions which should be re-evaluated when an 
external parameter {e.g., the renormalization scale) is modified on an event-by-event 
basis. The files couplingsXX . f where XX is greater than 1 contain all expressions 
which should be only re-evaluated if the default external parameter values are ex- 
plicitly read from the LHA parani_card . dat parameter card. The actual number of 
these files may vary, but a single file should be small enough to be compiled using 
standard Fortran compilers. The full list of these files should be included in the 
makefile through the makeinc . inc include file. 

• input . inc and coupl . inc: FORTRAN files containing all the necessary variable dec- 
larations. All parameters and couplings can be printed on screen or written to file 
using the routines defined in param.write . inc and coupl_write . inc, respectively. 
If needed, the latter can also be printed in a stricter format using routines defined in 
helas_couplings . f , so they can be used by external tools {e.g.. Bridge []64| ). 

Additional Fortran files, which are not model dependent, should also be provided in 
order to build the full library. Most of them simply include one or more of the above files, 
except lhajread.f which contains all the routines required to read the LHA format. A 
makefile allows the user to easily compile the whole package, to produce a library or a 
test program called testprog which can be used to quickly check the library behavior by 
producing a standard log output. 

The UsrMod v2 framework 

The UsrMod v2 framework has been designed as the successor of the widely-used 
original UsrMod template described in Ref. |^]. Taking advantage of the fixed structure 
we just defined, it provides the user with two new possibilities. First, any pre-existing 
model can be used as a starting point. This of course includes all models described in 
the present paper and soon part of the MadGraph/MadEvent distribution, but also all 
future models. This gives a natural framework for building simple extensions following a 
bottom-up approach, i.e., by adding successively new particles and new interactions and 
testing their implications at each step. Second, the possible modifications are no longer 
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restricted to the addition of new particles/interactions, but any alteration of the model 
content (including particle removal, modification of existing properties, . . . ) allowed in the 
context of MadGraph/MadEvent is supported in an user-friendly way. 

The UsrMod v2 approach can advantageously replace the full FeynRules package 
when only minor modifications to an existing MadGraph/Mad Event model are neces- 
sary, e.g., in order to study the phenomenology of a specific new particle and/or interaction, 
or when the use of the FeynRules machinery is not possible, e.g., if a Mathematica 
license is not available. However, when possible, we believe the use of FeynRules should 
be favored over the UsrMod v2 for the consistent implementation of full new models in 
MadGraph/Mad Event, especially due to the extended validation possibilities available 
in this context. 

Prom the implementation point of view, the UsrMod v2 package consists in various 
Python scripts (with one single main script) and works on any platform offering support 
for this programming language. The actual implementation of a new model is decomposed 
into four distinct phases: 

1. Saving: the model directory used as a starting point should be copied to a new 
location. The USRMOD script should be then run a first time to create a content 
archive used as a reference to identify the forthcoming modifications. 

2. Modifying: the particles . dat, interactions . dat and ident_card . dat files can be 
modified to arbitrarily add, remove or modify the particle, interaction and parameter 
content. 

3. Creating: the USRMOD script should be then run a second time to actually modify all 
the model files to consistently reflect the changes applied in the previous phase. 

4. Adjusting: the couplingsXX . f filc(s) can finally be edited, if necessary, to add or 
modify the relevant coupling expressions. The param_card.dat file can also be edited 
to modify default values of external parameters. 

At any time, the archive file created during the first phase can be used to restore the 
initial content of all model files. Several archive files can also be simultaneously saved into 
the same directory to reficct, for example, the successive versions of a single model. Finally, 
the intrinsic structure of the UsrMod v2 package favors various technical (not physical) 
consistency checks in the output files to minimize as much as possible the compilation and 
runtime errors. 

4.3.2 FeynRules interface 

The MadGraph/Mad Event interface can be called from the FeynRules package using 
the 

WriteMGOutput [ L ] 
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routine described in the FeynRules documentation^, where L is the name of the model 
Lagrangian. Since MadGraph/MadEvent currently only supports calculations in the 
unitary gauge, all the Goldstone and ghost fields are discarded in the particles.dat 
output, which is directly generated from the model description. After expanding all possible 
field indices {e.g., associated to flavor), an exhaustive list of non-zero vertices is generated 
and output as interactions.dat. If possible, the relevant coupling is extracted and, in 
case it does not already exist, stored in a new coupling variable of the form MGVXXX in 
a couplingsXX . f file. All the other required model-dependent files are finally generated, 
including the parain_card.dat where the default values (which are also the default value 
for the reading routines in paramjread. inc) are set as specified in the FeynRules model 
file, and where the QNUMBERS blocks correctly reflect the new particle content. All the 
produced files, together with the relevant model independent files are stored in a local 
directory model_nanie_MG, ready to be used in MadGraph/MadEvent. As mentioned 
previously, the testprog test program can be compiled and run to check the consistency 
of the created library. 

The two main restrictions of the MadGraph/MadEvent interface are related to the 
allowed Lorentz and color structures of the vertices. As already mentioned, even though 
FeynRules itself can deal with basically any interaction involving scalars, fermions, vec- 
tors and spin-two tensors, the Helas hbrary, used by MadGraph/MadEvent to build 
and evaluate amplitudes is more restricted. In the case no correspondence is found for a 
specific interaction, a warning message is displayed by the interface and the corresponding 
vertex is discarded. If this particular vertex is required for a given application, the user has 
still the possibility to implement it manually following the Helas library conventions and 
to slightly modify the interface files to deal with this addition. In case the vertex structure 
is also not present in MadGraph/MadEvent, a more involved manual modification of 
the code is also required. The second limitation of the present interface comes from the fact 
the color factor calculations are currently hardcoded internally MadGraph. While Feyn- 
Rules can deal with fields in any representation of the QCD color group, MadGraph 
itself is basically limited to the color representations appearing in the Standard Model and 
Minimal Supersymmetric Standard Model, e.g., a color sextet is not supported. 

Let us mention that work to alleviate both limitations is already in progress. The 
FeynRules package could, for example, be used to generate automatically missing Helas 
routines, while a more open version of the MadGraph matrix-element generator, e.g., 
taking advantage of a high-level programming environment, could advantageously deal 
with arbitrary color structures 

4.4 Sherpa 



Sherpa |11, [T^ is a general-purpose Monte Carlo event generator aiming at the complete 
simulation of physical events at lepton and hadron colliders. It is entirely written in 
C++ featuring a modular structure where dedicated modules encapsulate the simulation 
of certain physical aspects of the collisions. 

^Several options are available when running the interface. We refer for a detailed review to the Feyn- 
Rules manual. 
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Figure 2: Schematic view of tlie Sherpa's MODEL module, hosting the particle and parameter 
definitions of physics models as well as corresponding interaction vertices. 



The central part is formed by the hard interaction, described using perturbative methods. 
The respective generator for matrix elements and phase-space integration is AmegicH — h 
|65], which employs the spinor helicity formalism ||6^ , |6^ in a fully automated approach 
to generate processes for a variety of implemented physics models, see Sec. 1.4. 1| . Phase- 
space integration is accomplished using self-adaptive Monte Carlo integration methods 
||,|9|,0,0. 

The QCD evolution of partons originating from the hard interaction down to the hadroniza- 
tion scale is simulated by the parton-shower program Apacic+-|- |72]. It accounts for 
parton emissions off all colored particles present in the Standard Model and the Minimal 
Supersymmetric Standard Model. New shower generators have recently been developed 
in the framework of Sherpa 74 1. They account for QCD coherence and kinematics 



effects in a way consistent with NLO subtraction schemes, which makes them preferred 
over Apacic+-|-. They will be employed in future versions of Sherpa using an improved 
merging prescriptions for matrix elements and showers [^6| . 

An important aspect of Sherpa is its implementation of a general version of the CKKW 
algorithm for merging higher-order matrix elements and parton showers 
been validated in a variety of processes |3£ 



75, Ut 



], p. It has 
and proved to yield reliable results in 
comparison with other generators p7| , p8[ . 

Furthermore, Sherpa features an implementation of the model for multiple-parton inter- 



actions presented in Ref. [77|, which was modified to allow for merging with hard processes 



of arbitrary final-state multiplicity and eventually including CKKW merging [78|. Further 



more Sherpa provides an implementation of a cluster- fragmentation model |79], a hadron 
and tau decay package including the simulation of mixing effects for neutral mesons pO||, 



and an implementation of the YFS formalism to simulate soft-photon radiation [ 81 1 . 



4.4.1 Model framework 

Physics model definitions within Sherpa are hosted by the module MODEL. Here the particle 
content and the parameters of any model get defined and are made accessible for use 
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within the Sherpa framework. This task is accomphshed by instances of the basic class 
Model_Base. Furthermore the interaction vertices of various models are defined here that 
in turn can be used by AmegicH — h to construct Feynman diagrams and corresponding 
helicity amplitudes^. The corresponding base class from which all interaction models are 
derived is called Interaction_Model. A schematic overview of the MODEL module is given 
in Fig. ||. 

The list of currently implemented physics models reads: the Standard Model including 



effective couplings of the Higgs boson to gluons and photons |82|, an extension of the SM 



by a general set of anomalous triple- and quartic gauge couplings |83, U], the extension 



of the SM through a single complex scalar [M^], the extension of the Standard Model by 



a fourth lepton generation, the SM plus an axigluon [SC], the Two-Higgs-Doublet Model, 
the Minimal Supersymmetric Standard Model, and the Arkani-Hamed-Dimopoulos-Dvali 
(ADD) model of large extra dimensions |87, 88|, for details see Ref. [12|. Besides rou- 



tines to set up the spectra and Feynman rules of the models listed above corresponding 
helicity-amplitude building blocks are provided within AmegicH — h that enable the evalu- 
ation of production and decay processes within the supported models. In particular this 
includes all the generic three- and four-point interactions of scalar, fermionic and vector 
particles present in the SM and Minimal Supersymmetric Standard Model plus the eff'ective 
operators for the loop-induced Higgs couplings and the anomalous gauge couplings. The 
implementation of the ADD model necessitated the extension of the helicity formalism to 
interaction vertices involving spin- two particles [p9| . 

A necessary ingredient when dealing with the Minimal Supersymmetric Standard Model 
are specific Feynman rules for Majorana fermions or fermion number violating interactions. 
To unambiguously fix the relative signs amongst Feynman diagrams involving Majorana 
spinors the algorithm described in Ref. is used. Accordingly, the explicit occurrence 
of charge-conjugation matrices in the Feynman rules is avoided and instead a generalized 
fermion flow is employed that assigns an orientation to complete fermion chains. This 
uniquely determines the external spinors, fermion propagators and interaction vertices 
involving fermions. 

The implementation of new models in Sherpa in the traditional way is rather straight- 
forward and besides the public model implementations shipped with the Sherpa code 
there exist further private implementations that were used for phenomenological studies, 



cf. Ref. |91, 93]. From version Sherpa-1.2 on Sherpa will support model implementa- 
tions from FeynRules outputs - facilitating the incorporation of new models in Sherpa 
further. 

4.4.2 FeynRules interface 

To generate FeynRules output to be read by Sherpa, the tailor-made FeynRules rou- 
tine 

WriteSHOutput [ L ] 



'Note that within Sherpa Feynman rules are always considered in unitary gauge. 



-26- 



has to be called, resulting in a set of ASCII files that represent the considered model through 
its particle data, model parameters and interaction vertices^. 

To allow for an on-the-flight model implementation from the FeynRules outputs, in- 
stances of the two basic classes Model_Base and Interaction_Model_Base are provided 
dealing with the proper initialization of all the particles and parameters, and the interac- 
tion vertices of the new model, respectively. The actual C-|--|- classes for these tasks are 
called FeynRules_Model and Interaction_Model_FeynRules, see Fig. |2|. 
The master switch to use a FeynRules generated model within Sherpa is 

MODEL = FeynRules 

to be set either in the (model) section of the Sherpa run card or on the command line once 
the Sherpa executable is called. Furthermore the keywords FRJ' ARTICLES, FR_IDENTFILE, 
FR_PARAMCARD, FR_PARAMDEF and FR.INTERACTIONS, specifying the names of corresponding 
input files, need to be set. The actual format and assumed default names of these input 
cards will be discussed in the following: 

• FR_P ARTICLES specifies the name of the input file listing all the particles of the theory 
including their SM quantum numbers and default values for their masses and widths, 
default name is Particle . dat. An actual particle definition, e.g., for the gluon, looks 
like 

kf Mass Width 3*e Y SU(3) 2*Spin maj on stbl m_on IDName TeXName 
21 0. .0 8 2 -1 1 1 G G 

Hereby kf defines the code the particle is referred to internally and externally, typi- 
cally its PDG number [^]. The values for Mass and Width need to be given in units 
of GeV. The columns 3*e and Y specify three times the electric charge and twice 
the weak-isospin. SU(3) defines if the particle acts as a singlet (O), triplet (3) or 
octet (8) under SU{3)c. 2*Spin gives twice the particle's spin and maj indicates if 
the particle is charged (O), self-adjoint (-1) or a Majorana fermion (l). The fiags on, 
stbl and m_on are internal basically and define if a particle is considered/excluded, 
considered stable, and if its kinematical mass is taken into account in the matrix- 
element evaluation. IDName and TeXName indicate names used for screen outputs and 
potential I^T^X outputs, respectively. 

• In FR_IDENTFILE all the external parameters of the model get defined, default file 
name is ident_card.dat. Names and counters of corresponding parameter blocks to 
be read from FR_PARAMCARD are listed and completed by the actual variable names 
and their numerical types, i.e., real R or complex C. Besides, variable names for all 
particle masses and widths are defined here. To given an example, the section defining 
the electroweak inputs of the SM may look like 

''Note again that Feynman rules have to be considered in unitary gauge. 
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SMINPUTS 1 aEWMl R 
SMINPUTS 2 Gf R 
SMINPUTS 3 aS R 
CKMBLOCK 1 cabi R 



• In the file specified through FR_PARAMCARD the numerical values of all elementary pa- 
rameters, particle masses and decay widths are given, default file is param_card.dat. 
Following the example above the electroweak inputs of the SM can be set through: 



Block SMINPUTS 

1 1.2790000000000E+02 # aEWMl 

2 1.1663900000000E-05 # Gf 

3 1.1800000000000E-01 # aS 
Block CKMBLOCK 

1 2.2773600000000E-01 # cabi 



• FR_PARAMDEF gives the file name where all sorts of derived parameters are defined, 
default param_definition.dat. Such variables can be functions of the external pa- 
rameters and subsequently other derived quantities. A few examples for the case of 
the SM again might read: 



aEW = pow (aEWMl, -1.) R 

G = 2.*sqrt(aS)*sqrt(M_PI) R 

CKMll = cos (cabi) C 



Electroweak coupling constant 
Strong coupling constant 
CKM-Matrix ( CKMll ) 



The parameter definitions get interpreted using an internal algebra interpreter, no 
additional compilation is needed for this task. All standard C-|— |- mathematical 
functions are supported, e.g., sqr, log, exp, abs. For complex valued parameters, 
e.g., CKMll, the real and imaginary part can be accessed through Real (CKMll) and 
Iinag(CKMll), the complex conjugate is obtained through Conjugate (CKMll). 

• The keyword FR_INTERACTIONS (lnteractions.dat) specifies the input file contain- 
ing all the vertices of the considered model in a very simple format: 

VERTEX 21 21 21 # G G G 

1 G # right-handed coupling 

2 G # left-handed coupling 

3 F [1,2, 3] # colour structure 

4 Gauge3 # Lorentz structure 

The keyword VERTEX signals the start of a new Feynman rule followed by the PDG 
codes of the involved particles. Note, the first particle is always considered incoming 
the others outgoing. Counters number 1 and 2 indicate the right and left-handed 
coupling of the vertex rule, the right and left-hand projector being given by Pr/l = 



-28- 



2(1 ± 75), respectively. Couplings are given in terms of the elementary and derived 
parameters. Counter number 3 explicitly gives the color structure of the interaction 
in terms of the SU{3)c structure constants or generators. The spin structure of the 
vertex is given under 4, identified through a keyword used by Sherpa to relate a 
corresponding sub-amplitude to the correct helicity-amplitude building block. 

Sherpa's interface to FeynRules is designed to be as general as possible, it is, how- 
ever, by construction restricted in two ways. 

The functional form of the model parameters, and respectively the couplings, is limited 
by the capabilities of the algebra interpreter that has to construe them. This limitation, 
however, might be overcome by using an external code to calculate the needed variables 
and redefining them as external giving their numerical values in FR_PARAMCARD. 
More severe limitations originate from the restricted ability of Sherpa/ AmegicH — h to han- 
dle new types of interactions. Only three and four-point functions can be incorporated. For 
the color structures only the SU{3)c objects 1, Sij, 6a,b, T^j, f"'^'^ and products of those, 
e.g., y^^'^ycde^ supported. Lorentz structures not present in the SM or the Minimal 
Supersymmetric Standard Model are currently not supported by the interface. Further- 
more, Sherpa cannot handle spin-3/2 particles. QCD parton showers are only invoked for 
the colored particles present in the SM and the Minimal Supersymmetric Standard Model. 
Hadronization of new colored states is not accomplished, they have to be decayed before 
entering the stage of primary hadron generation. 

5. Models 

In this section we briefly present the implementation of the Standard Model and other 
several important New Physics models in FeynRules. Our aim is to show that very 
complete and sophisticated implementations are possible and that FeynRules offers a 
very natural and convenient framework where models can be first developed (from the 
theoretical point of view) and then tested/constrained against experimental data. Since 
the main focus is the implementation procedure, the actual model descriptions, as well as 
the information about values of parameters, are kept to the minimum. More exhaustive 
information about each of the following models, all of which being publicly available, can 
be found on the FeynRules website. 

5.1 The Standard Model 
5.1.1 Model description 

As it serves as basis to any new bottom-up implementation, we briefly describe here the 
Standard Model implementation. The SM of particle physics is described by an SU{'3>)c x 
SU{2)l X U{\)y gauge theory, where the electroweak symmetry is spontaneously broken 
so that the fundamental fermions and the weak gauge bosons acquire a mass. The particle 
content of the SM is summarized in Table ^. The Lagrangian can be written as a sum of 
four parts, 

^SM = -^^SM, Gauge + -^SM , Fermions , Yukawa- 

(5.1) 



-29- 



r Oil LiLic 


spill 


XVtJpi cooiiLcxLlUllo 


rO /,,0 ;0 \T 

,/0 

^iR 

"■iR 


1 /o 

1/2 
1/2 
1/2 

1 /9 
i/2 

1/2 


(1,2,-1/2) 
(3,2,1/6) 

(3, i, Z/6) 

(3,1,-1/3) 


$ 





(1,2,1/2) 




1 
1 
1 


(1,1,0) 
(1,3,0) 
(8,1,0) 



Table 4: The SM fields and their representations under the Standard Model gauge groups SU{3)c x 
SU{2)l X U{1)y. 



The pure gauge sector of the theory reads 



(5.2) 



where the SM field strength tensors are defined following the conventions introduced in 
Eq. (|2.4|). The Lagrangian describing the matter fermions can be written as 



/:SM,Fermions = 0° Q! + ^ ilp + tIP U% + < ilp 4, + ^ „ (5.3) 

where Z)^ denotes the SU{3)c x SU{2)l x U{1)y covariant derivative, and we use the con- 
ventions of Eq. (12.31). The superscript refers to the gauge eigenstates. Note in particular 
that explicit mass terms for the matter fermions are forbidden by gauge symmetry. The 
Higgs field is described by the Lagrangian 



/:SM,Higgs = i?^ <l>t Z)^ <!> - ^2 <J> _ A (cDt CD) 



(5.4) 



li jj? < 0, then the Higgs field acquires a vev that breaks the electroweak symmetry 
spontaneously. Expanding the Higgs field around its vev. 



(5.5) 



we generate mass terms for the Higgs boson H and the electroweak gauge fields. The mass 
eigenstates for the gauge bosons are the W and Z bosons, as well as the photon, which 
remains massless. The relations between those fields and the original SU{2)l x U{\)y 
gauge fields are 



A.. 



B„ 



(5.6) 
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where we introduced the weak mixing angle 



= cos^iu = — — . (5.7) 
Mz 

The interactions between the fermions and the Higgs field are described by the Yukawa 
interactions 

/:sM,Yukawa = -u'Ir Q° l> - S^r y% Q] - 4 h] $ + h.c. , (5.8) 

where <^ = ia^ After electroweak symmetry breaking the Yukawa interactions generate 
non-diagonal mass terms for the fermions that need to be diagonalized by unitary rotations 
on the left and right-handed fields. Since there is no right-handed neutrino, we can always 
rotate the leptons such that the mass matrix for the charged leptons becomes diagonal 
and lepton flavor is still conserved. For the quarks however, the diagonalization of the 
mass matrices introduces flavor mixing in the charged current interactions, described by 
the well-known CKM matrix. 

5.1.2 FeynRules implementation 

The SM implementation in FeynRules is divided into the four Lagrangians described 
in the previous section. In particular, one can use the dedicated functions for the field 
strength tensors and the covariant derivative acting on the left and right-handed fermions. 
Matrix-element generators however need as an input the mass eigenstates of the particles, 
and therefore it is mandatory to rotate all the gauge eigenstates into mass eigenstates 
according to the prescriptions discussed in the previous section. This can be done very 
easily in FeynRules by writing the Lagrangian in the gauge eigenbasis, and then let- 
ting FeynRules perform the rotation into the mass eigenstates (note, at this point, that 
FeynRules does not diagonalize the mass matrices automatically, but this information 
has to be provided by the user). However, as the SM Lagrangian is the starting point for 
many bottom up extensions, the actual implementation was performed directly in terms of 
the fermion mass eigenstates. The benefit is a slight speed gain due to the rotations in the 
fermion sector. The default values of the external parameters in the model file are given 
in Table |^ in Appendix [B.l| . 

Three restriction files for the SM implementation are provided with the default model 
distribution: 

• Massless . rst: the electron and the muon, as well as the light quarks (u, d, s) are 
massless. 

• DiagonalCKM.rst: the CKM matrix is diagonal. 

• Cabibbo.rst: the CKM matrix only contains Cabibbo mixing. 

Another particularity of the SM implementation is that it was performed both in uni- 
tary and in Feynman gauge. The model file contains a switch FeynmanGauge, which, if 
turned to False, puts to zero all the terms involving ghost and Goldstone fields. The 



- 31 - 



default value is False. 



Possible extensions 

The SM is at the basis of almost all BSM models, and thus the number of possible 
extensions of the SM implementation is basically unlimited. A first extension of this model 
was presented in Section ^ with the HAH model, based on the simplest possible extension 
of the gauge sector of the SM. Other possibilities are the addition of higher-dimensional 



operators compatible with the SM symmetries (see Section 5.6) or the inclusion of right- 
handed neutrinos via see-saw models. 

5.2 The general Two-Higgs-Doublet Model 

The Two-Higgs-Doublet Model (2HDM) has been extensively studied for more than twenty 
years, even though it has often been only considered as the scalar sector of some larger 
model, like the Minimal Supersymmetric Standard Model or some Little Higgs models 
for example. The general 2HDM considered here already displays by itself an interesting 
phenomenology that justifies its study like, for example, new sources of CP violation in 
scalar-scalar interactions, tree- level flavor changing neutral currents (FCNC's) due to non- 
diagonal Yukawa interactions, or a light pseudoscalar state and unusual Higgs decays (see 
Ref. m). 



5.2.1 Model description 

The 2HDM considered here is based on two SU{2)l doublets 0i and (f)2 with the same 
hypercharge Y = +1. If one imposes only gauge invariance, the most general renormalizable 
Lagrangian is composed of four parts, 

'C2HDM = 'C2HDM, Gauge + '^^2HDM,Fcrmions + 'C2HDM,Higgs + '^^2HDM,Yukawa- (5.9) 

The gauge and fermion sectors of the model are identical to the SM, 

-^^2HDM,Gauge = -^^SM , Gauge 

and -C2HDM ,Fermions ,Fermions- 

(5.10) 

The Lagrangian of the Higgs sector differs from the SM, and can be written 

/:2HDM,Higg« = Df, 4 01 + 4 cj)2 - V{<Pu<P2), (5.11) 

and the scalar potential reads, in the notation of Ref. |95], 

V{(t>l,(t>2) =/Xl0l(/'l +/X24</'2 + (/i3<Al'A2 +h.C.) +Ai (^^l^lY + X2 (4<^2)^ 

+ A3 (0l<Ai) {4^2) + A4 {4i(t>2) (4</'i) (5-12) 



+ 



A50l(/'2 + Ae^I^i + X742h) (0i'/>2) + h.c. 



where /ii^2 and Ai^2,3,4 ai'e real parameters while /i3 and X^^j are a priori complex. We 
assume that the electromagnetic gauge symmetry is preserved, i.e., that the vevs of and 
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(/>2 are aligned in the SU (2) l space in such a way that a single SU (2) l gauge transformation 
suffices to rotate them to the neutral components, 



and (02) = ^ , J , (5.13) 





with vi and V2 two real parameters such that v\ + V2 = v"^ = {\/2Gp) ^ and V2/V1 = tan/3. 
The most general form for the Yukawa interactions of the two doublets reads 

QIV2 r, 



2HDM, Yukawa 



V 

E2V2 



V 

.2, 



[{Ae(l)l+Te(p2)eR], 



(5.14) 



with = ia (j)* and where the 3x3 complex Yukawa coupling matrices Aj and Fj are 
expressed in the fermion physical basis, i.e., in the basis where the fermion mass matrices 
are diagonal. We choose as free parameters the Fj matrices, while the other Yukawa 
couplings, the A, matrices, are deduced from the matching with the observed fermion 
masses. Conventionally, the two indices a and h of the elements of the Yukawa matrices 
(Fj)^^ and (Aj)^^^ refer to the generations of the SU{2)l doublet and singlet, respectively. 

5.2.2 FeynRules implementation 

The 2HDM Lagrangian implemented in FeynRules is composed of Eqs. ( ^.12 ), ( 5.14| ), 



together with the canonically normalized kinetic energy terms for the two doublets and the 
other SM terms. An important feature of this model is the freedom to redefine the two 
scalar fields (j)i and 02 using arbitrary U{2) transformations 




= Uy^lj ' U(^U{2) (5.15) 
since this transformation leaves the gauge-covariant kinetic energy terms invariant. This 




notion of basis invariance has been emphasized in Ref. [95| and considered in great detail 
more recently in Refs. Since a given set of Lagrangian parameter values is only 

meaningful for a given basis, let us take advantage of this invariance property to select 
the Higgs basis (by defining the additional file HiggsBasis.fr) where only one of the two 
Higgs fields acquires a non-zero vev, i.e., 

{H'i) = ^ and {Hl) = 0. (5.16) 

Let us note that the Higgs basis is not univocally defined since the reparametrization 
H2 — > e^"H2 leaves the condition Eq. ( 5.16D invariant, so that the phase of H2 can be 



fixed in such a way that A5 becomes real. Other basis choices can in principle be easily 
implemented as different extension files for the main Lagrangian file Lag.fr. 

The minimization conditions for potential of Eq. ( ^.12| ) read, in the Higgs basis defined 
in Eq. (|5T^ ): 

1^3 = -^6y , (5.17) 
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which reduces the number of free parameters in the most general 2HDM to ten (seven 
real parameters, three complex ones and three minimization conditions). Besides the usual 
three massless would-be Goldstone bosons, the physical spectrum contains a pair of charged 
Higgs bosons with mass 

= -y- + /^2, (5.18) 
and three neutral states with the squared mass matrix 



/ 2Xiv^ Re{XQ)v^ -Im(A6)v2 

ReiXev^) + (A4/2 + As)?;^ | . (5.19) 

y -lm{XQ)v^ m^± + (A4/2 - X^y 



The symmetric matrix is diagonalized by an orthogonal matrix T. The diagonalization 
yields masses rrii for the three physical neutral scalars of the model (where the index i 
refers to mass ordering), 

= Tdiag(mf , m^, m^)r^ . (5.20) 



The doublet components are related to these physical states through 

(5.21) 



Re(i/?) \ 
Im(FO) / 



V2 
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The Yukawa couplings of the model are fully determined by the Fj matrices in Eq. ( 5.14| ) 



since the Aj are, by definition, fixed to the diagonal fermion mass matrices in the Higgs 
basis. 

In the current implementation of the 2HDM into FeynRules, the user has to provide 
numerical values for all the Aj parameters in the basis of Eq. ( ^.16| ), together with the 
charged Higgs mass mjj±. The other parameters of the potential, such as the 
then deduced using; Eqs. ( CT) and (|5l^ ). As a consequence, the orthogonal matrix 



T must be calculated externally. This, together with the change of basis required if the 
user wants to provide potential parameters and Yukawa coupling values in bases different 
from this of Eq. ( ^.16| ), can be done using the TwoHiggsCalc calculator introduced in 
Ref. Q which has been modified to produce a parameter file compatible with the present 
implementation. This calculator can also be used to calculate the required Higgs boson 
tree- level decay widths. 

5.3 The most general Minimal Supersymmetric Standard Model 

Most present supersymmetric models are based on the four-dimensional supersymmetric 
field theory of Wess and Zumino 1 9£ ] . The simplest model is the straightforward supersym- 



metrization of the Standard Model, with the same gauge interactions, including i?-parity 
conservation, and is called the MSSM [|lO[l| , 101]. Its main features are to link bosons with 



fermions and unify internal and external symmetries. Moreover, it allows for a stabilization 
of the gap between the Planck and the electroweak scale and for gauge coupling unification 
at high energies, provides the lightest supersymmetric particle as a dark matter candidate 
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and appears naturally in string theories. However, since super symmetric particles have 
not yet been discovered, supersymmetry must be broken at low energies, which makes the 
superpartners heavy in comparison to their Standard Model counterparts. 

Supersymmetric phenomenology at colliders has been extensively investigated for basic 
processes at leading order [|lO|, |o|, |l05|, |lO|, [lO^, |o|, |lO|, |ll^, |lT|, |ll3| 



and next-to-leading order pl4| , pl5| , p6| , pT7| , |rT|, |120| of perturbative QCD. More 
recently, for some processes, soft-gluon radiation has been resummed to all orders in the 
strong coupling constant and the results have been matched with the next-to-leading order 
calculations |121, 122, |123| , 124, 125, |126[| . However, even if those calculations are useful for 
inclusive enough analyses, they are not suitable if we are interested in a proper description 
of the full collider environment, for which Monte Carlo event generators are needed. For 
a couple of years, all the multi-purpose generators already mentioned contain a built-in 
version of the MSSM. The model files for FeynArts/FormCalc are described in Ref. 
1 12^ , [T28| , for CalcHep in Ref. gH], for MadGraph/MadEvent in Ref. |30l, and for 
Sherpa in Ref. il3l|1 . The Sherpa and FeynArts/FormCalc implementations keep 
generic mixing in the scalar sector while the other generators rely on a simplified model 
with only helicity mixing for third generation sfermions. 

Our MSSM implementation in FeynRules is the most general one in a sense that 
it is keeping all the flavor-violating and helicity-mixing terms in the Lagrangian and also 
all the possible additional CP-violating phases. This yields thus 105 new free parameters 
1 132], and in order to deal in a transparent way with all of those, our implementation will 
follow the commonly used universal set of conventions provided by the Supersymmetry 



Les Houches Accord (SLHA) [133, |134j] , except for some minor points. We will dedicate a 
complementary paper to a complete description of the model 135|. 



5.3.1 Model description 
Field content 

Each of the Standard Model quarks and leptons is represented by a four-component Dirac 
spinor f^, where i stands for the generation index and the superscript denotes inter- 
action eigenstates. It has two associated scalar superpartners, the sfermion and the 
antisfermion being related to the two-component holomorphic Weyl fermion Xfi and 
antifermion Xfi^ respectively. Let us recall that we relate the Dirac fermion representations 
to the Weyl ones by 

/° = 1^^^; j . (5.22) 

For clarity, we will use in the following the left-handed component f^- and the right-handed 
component f^^- of the Dirac fermion and not the Weyl fermions x fi a-nd Xfi- To preserve 
the electroweak symmetry from gauge anomaly and in order to give masses to both up-type 
and down-type fermions, the MSSM contains two Higgs doublets Hi, together with their 



-35- 



Name 


Particle 


spin 


Superpartner 


spin 


Representations 


(s) quarks 


\AUt A.CLI J 


1/2 







fS 2 1/6) 




X.Wi 


1/2 

1 


R 





(3,1,-2/3) 






1/2 


K 





(3,1,1/3) 


(s)leptons 


iXvi XliV 


1/2 







(1,2,-1/2) 




Xii 


1/2 









Higgs(inos) 


(H+ H^f 







1/2 


(1,2,1/2) 











1/2 


(1,2,-1/2) 


i3-boson, bino 




1 




1/2 


(1,1,0) 


VF-bosons, winos 




1 




1/2 


(1,3,0) 


gluon, gluino 


9 


1 


^9 


1/2 


(8,1,0) 



Table 5: The MSSM fields and their representations under the Standard Model gauge groups 
S'[/(3)c X SU{2)l X U(1)y- The quarks and leptons are denoted in terms of two-component Weyl 
spinors. 




fermionic partners, the higgsinos ipHi, 

Finally, the spin-one vector bosons of the Standard Model will be associated to Majorana 
fermions, the gauginos ipB, i'w'' ^'^d ^9- The names and representations under the Stan- 
dard Model gauge groups SU{2>)c x SU{2)l x U{1)y of the various fields are summarized 
in Table |5[ 

The full MSSM Lagrangian can we written as 

-^MSSM = -CmSSM, Gauge + -^MSSMTermions + -^MSSM.Yukawa 

,Scalar kinetic + ^MSSM.Scalar FDW (5.24) 
,Ino kinetic + -^MSSM .Ino Yukawa ,Ino mix 



Starting from the expression of the Lagrangian in the gauge-eigenstate basis of fields given 
above, we diagonalize the non-diagonal mass matrices arising after electroweak symmetry 
breaking and provide transformation rules allowing to re-express the Lagrangian in the 
physical basis. 



Supersymmetry-conserving Lagrangian In order to have more compact notations. 
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we introduce the SU (2) ^-doublets of left-handed fermions and sfermions, 



u 





Qi = \ ::\ ^ L'i={ : \ , Q^={ " and ^={:\. (5.25) 




The pure gauge sector and the matter fermion section of the MSSM are identical to the 
Standard Model, 

■^MSSM,Gauge = ■^SM,Gauge and /^MSSM,Fermions = .CsM,Fermions- (5.26) 

The Lagrangian for the scalar sector can be divided into one purely kinetic part, 

>CMSSM,Scalar kinetic = D^^HI D^^Hu + D^H^ D^'Hd + D^.V''^ Dt'Ll + D^JP^^ Dt^f^R, 



(5.27) 



and one part derived from the D-terms, the F-terms and the superpotential W. Indeed, 
auxiliary F- and D-fields must be added to the theory in order to preserve supersymmetry 
when considering off-shell states, thereby keeping the total number of fermionic and bosonic 
degrees of freedom equal. Solving their equations of motion leads, for a set of n Dirac 
fermions {'i/'n} and 2ra associated supersymmetric scalar partners {0Ln,^^^}, to 



■^MSSM.Scalar FDW = " WlW^^ - Wj^Wj^- 

1 r 

~ 2 
1 

~ 2 L- 



W'l^ i^tPL^j + 2<^^ i^iPL^Pj + wHj, ^^Pl^^ + h.c, 

g {4T''4>Li) - g (4r>ffi)l [ - g {4>'^T-4>Lj) + g (4^r>fl,-) 



(5.28) 



where ip^ is the field charge-conjugated to ip and the derivatives of the superpotential W 
are given by 

dW dW 
Wl=^ and W'r 



""^^-dHidHj ' dHid<^nj ''''' d<Pnid^ni 

Let us note that in the case of Major ana fermions, there is only one associated scalar field. 
In the framework of the MSSM, the superpotential reads 

Wmssm = (y")/' (Q°e//,)-4^ {y%' iQ'jeHd)-f^^ (y')/' meHa)+iiH^eHa, (5.30) 

where y", y'* and y' denote the 3x3 Yukawa matrices, /i the Higgs off-diagonal mass- 
mixing, and e the SU{2)l invariant tensor. We could add to this superpotential other 
gauge-invariant and renormalizable terms. 
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with the Yukawa-hke couphngs A, A' and A", and a slepton-Higgs off-diagonal mass term k. 
Those couphngs would however violate either the lepton number L or the baryon number 
B, as well as the individual lepton flavors. Moreover, they allow for various 5-violating 
or L-violating processes which have never been seen experimentally. We could just forbid 
these terms by postulating B and L conservation, but neither B nor L are fundamental 
symmetries of nature since they are violated by non-perturbative electroweak effects [136| . 
Therefore, an alternative symmetry is rather imposed, the -R-parity |137], defined by 



R=i-1) 



3B+L+2S 



(5.32) 



S being the spin of the particle, forbidding any term different from those in Eq. ( 5.3C| ). 
All the Standard Model particles have thus a positive i?-parity while the superpartners 
have a negative one. In this paper, we will only describe the implementation of the R- 
parity conserving MSSM. However, the ii-parity violating extension of our implementation 



is straightforward, available and described in Ref. |135]. We can now write the remaining 
part of the scalar Lagrangian, 



-^MSSM.Scalar FDW 



_i_ U' ^-.O^t /-.u\ J /)0 
+ H B^ u^^ [y )i Qj 



-\l2\^{\Hu\^ + \Bd\^) 
+ [HieQ''^] [y"V"]/ [Q'.eH^] + [BleQ^^^ [y'^V]/ [Q'.eH,] 



+ h.c. 



u 



R 



1 ,,1+1 i TO 



-[L'^Hy%'pRi\ [pR'Hy\'L[ 



l\ I ?01 



dU i JO 
j (J-Ri 



[fSHy\'L^ 

g0.t qO + ^ lO^^ 10 + | # t 



[L''Hy''vpRi\ [d'^'Hy^'Q'i 



,d\ I r\Ol 



HlHu 



+ ^ ,-}0*t ,-.0 _ 9^ m JO 
+ — Uj^ Ujii — o(d a, 



''R ^Ri 



J m TO 

9 l-R I'Ri 



9 



g0it^fcg0_5^0*t ^klO 



9 fj\ „k 



^B\a-H,-^-Bla^B^ 



-Oif rjia ~0 



"R 



-a JO 
^Ri 



(5.33) 



From the W^^ terms of the superpotential, bilinear in the fermionic fields, we can also 
generate the Yukawa couplings between the matter fermions and the Higgs fields. 



C 



MSSM.Yukawa 



-U 



°'^y")/ Q° e Bu + 4*(y^)/ Q', e B, + /^^(y')/ L'^ e B, + h.c. (5.34) 



;iOif,rd\ i nO 



lOu.M j tO . 
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The terms of the MSSM Lagrangian containing higgsino and gaugino fields can be divided 
into three parts; pure kinetic terms, 



-CMSSMJno kinetic = 2 V'B ^0 ^pB + '^'tpWl' '^W' + ^'^g'' ^-0 V'g<' + 2 V'H„ '^Hu + -ji^Ha iJp tpHa, 

(5.35) 

Yukawa interactions obtained from the superpotential terms Wij, 



Ano Yukawa = - [4^ ^ J (y")/ Q° + [4 6 J (Y^)/ Q° + [Z^ e J (y')/ L'j 



(5.36) 



and additional super symmetry-conserving gauge-like interactions which have no counter- 
part in the Standard Model and which are not taken into account through the covariant 
derivatives, 



>CMSSM,Ino mix =V2g'[- ^Q^^^ {Qf i^s) + {Lf ^Pb) 



+ V2g' 
_ _9_ 

_ _g_ 

V2 



\hI {{pj,^ ^b) + {rn, V'b) 



-l-h.c. 

(5.37) 



Supersymmetry-breaking Lagrangian 

As stated in the introduction, the masses of the superpartners must be considerably 
larger than those of the Standard Model particles. Realistic supersymmetric models must 
hence include supersymmetry breaking, which is expected to occur spontaneously at some 
high scale. The Lagrangian density then respects supersymmetry invariance, but the vac- 
uum state does not. Moreover, in order not to introduce quadratic divergences in loop- 
calculations, supersymmetry has to be broken softly. In practice, since we do not know the 
supersymmetry-breaking mechanism and the corresponding scale, we will add all possible 
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terms breaking supersymmetry explicitly at low-energy ||13. 
1 r 



MSSM,Soft 



MiTpB 1pB + M2 V'vy* V'VKfe + ^3 "03" '(Pg'^ 



+ 



~0i\ 



U 



Tu)/ Q," eHu + (Td)/ Q^j eHd + ^ (Ti)/ ^ eH^ + h.c 



R 



(5.38) 



The first line of Eq. ( ^.38 ) contains gaugino mass terms, the second and third lines the 
sfermion mass terms, where m^, m?, m?, m?, m? are 3x3 hermitian matrices in gener- 

Q L d ^ 

ation space, the fourth line mass terms for the Higgs fields, and the fifth line the trilinear 
scalar interactions, Tu, Tj, and Te being also 3x3 matrices in generation space. Let 
us note that the additional Higgs mass terms are required in order to break electroweak 
symmetry spontaneously. We remind the reader that for the i?-parity violating MSSM, 
additional soft supersymmetry-breaking terms must also be included, 



C 



MSSM, i? Soft 



+ 



\t% LO e I] + T'\ l\ e <F^^ + uf d^'^ + h.c. 



(5.39) 



Particle mixing 

In order to break the electroweak symmetry to electromagnetism, the classical Higgs 
potential, i.e. all the Lagrangian terms quadratic or quartic in the Higgs fields, must have 
a non-trivial minimum. Due to gauge invariance, one of the charged Higgs vev can be 
rotated away, which yields a zero vev for the other charged Higgs. Again, we can use 
gauge-invariance to choose the remaining vevs real and positive, so that we can replace the 
neutral Higgs fields in £mssm by 

^ ^ and flS ^ '-i±3, (5.40) 

where Vu and are the two vacuum expectation values of the neutral Higgses and and 
are complex scalar fields. Let us note that those relations are only valid in unitary gauge, 
which is the case we are dealing with here. We can then extract mass matrices for the gauge 
bosons and W'^, diagonalize them, and derive the physical mass-eigenstates, the photon 
A, the W and Z boson. The transformation rules relating the mass and interaction bases 
are given by Eq. (^]^). The weak mixing angle 9y^ and the physical masses Mz and Mw 
are defined by 



^2 



^°^'^- = ^2^ ' ^^ = ^^^V^« + ^1 and Mw = yvl + vj. (5.41) 



-40- 



In the Higgs sector, three out of the eight real degrees of freedom of the two doublets are the 
would-be Nambu-Goldstone bosons becoming the longitudinal modes of the weak bosons, 
while the five others mix to the physical Higgses, h^, H^, and H^. The diagonalization 
of several mass matrices leads to the transformation rules 

h^. = cos ah^ + sin a + zcos B and /i? = —sin ahP + cos a + isin 3 A^, 

(5.42) 

H+ = cosPH+ and H~=sm(3H-, 

where tan /? is the ratio of the two vevs Vu and v^- 

Collecting the terms of >Cmssm bilinear in the Majorana fermions ips, V'lys, ipfjo and 
wo , we can extract the neutral gaugino-higgsino mass matrix Y which can be diagonalized 

d 

through a unitary matrix A^, 

N* Y N^'^ = diag (m^o , m^o , m^o , m^o ) . (5.43) 
This matrix relates the four physical neutralinos to the interaction-eigenstates, 

(x? xl xl7 = N{^B i>w^ i>Hl i^Hof- (5.44) 

Similarly, we can collect the terms yielding the charged gaugino-higgsino mass matrix X 
and diagonalize it through the two unitary matrices U and V, 

U* XV-^ = diag{m-±,m-±). (5.45) 
Those matrices relate the interaction-eigenstates to the physical charginos xf^ according to 
ixt xtf = V{iJw+ i^H+f and (xr X2f = U{i^w- i^n-f- (5-46) 



Let us note that the fields iI>\y± obtained after rotating ijjiyi and V'vk^ as in Eq. ( p^ ) 
for the W boson. 

Within the Standard Model, the only source of flavor violation arises through the rota- 
tion of the up-type and down-type quark interaction-eigenstates basis {(fli,u^i, f^ffij ^ffi} to 
the basis of the physical mass eigenstates {dnjUu, dm^Ufn} in which the Yukawa matrices 

and y'^ are diagonal. These diagonalizations require four unitary matrices 

d^Li = dLi , dffi = Uddji, uli = Vu ULi , u%i = UuUR, (5.47) 

and render the charged-current interactions proportional to the CKM matrix, 

^CKM = V^Vd . (5.48) 

The leptonic sector is however diagonalized with the help of only two matrices, since the 
neutrinos are assumed massless, 

lLi = VilLi, l%i = UilR, uf = iy,. (5.49) 

Let us note that the matrices Uu, U^, Vi and Ui do not appear in the rotated Lagrangian. 
In the sfermion sector, we define the super-CKM basis [139] as the basis in which the 
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sfermion interaction-eigenstates undergo the same rotations as their fermionic counter- 
parts. As a consequence, the squark charged-current interactions are also proportional 
to the CKM matrix. However, the fermion and sfermion fields can be misaligned due to 
possible off-diagonal mass terms in the supersymmetry-breaking Lagrangian >CMSSM,Soft- 
The diagonalization of the four mass matrices M?, Mj, Mj and M? requires thus the 
introduction of three 6x6 matrices i?", R'^ and i?' and one 3x3 matrix i?^, 

diag (mi^ , . . . , mij = fl" M? R^^ , diag {m\ , • • • , m| ) = i?^ M? R"^ , (5.50) 
diag(m?^,...,m?3) = R" M? R"^ , diag (m? , . . . , m|) = R^ Uf R^K (5.51) 

These matrices relate the physical mass-eigenstates to the interaction-eigenstates through 

{ui,U2,U3,U4,U5,Uq)'^ = {ul, Cl, h, Ur, Cr, tn)'^ , 
(^1,^2,^3, (4, 4, 4)^ = R''- {dL,SL,bL,dR,SR,bR)'^, 
(1^1,1^2,1^3) =R {l'e,T^fj.,T^T) , 
([1, [2, ^"3, k, k, k)^ = R^ ieL,P'L,fL, BR, p,R, trY. 

The less general built-in implementation of the MSSM in MadGraph/MadEvent and 
CalcHep can easily be recovered by setting all the off-diagonal elements of these general 
mixing matrices to zero, except for the third generation flavor-conserving and helicity- 
mixing elements. In a framework of ii-parity violating scenarios, additional mixings be- 
tween charged leptons and charginos, neutrinos and neutralinos, Higgses, sleptons and 
sneutrinos can arise. Moreover, the sneutrino fields can acquire vevs, since they are not 
protected by lepton number conservation, unlike in the i2-parity conserving MSSM where 



it is a conserved quantum number. We refer the reader to Ref. |135] for more information. 

5.3.2 FeynRules implementation 

Current implementation 

We describe here the implementation of the most general i?-parity conserving MSSM 
in FeynRules. We divide the complete Lagrangian in six pieces, 

.Gauge ,Fermions + -^MSSM .Scalars ,_ 

, (5.53) 

,Inos 



related to Eq. (|5.24| ) through 



-^MSSM, Scalars = -^MSSM, scalar kinetic, 

'CMSSM.Inos = -^^MSSMJno kinetic + -^^MSSM.Ino mix, (5.54) 
-^^MSSM.FD = -^^MSSM, Scalar FDW + -^MSSM .Yukawa ,Ino Yukawa • 

-^MSSM.Scaiars Contains the scalar kinetic terms and gauge interactions, all the F-terms and 
Z)-terms are embedded in -Cmssm,fd, and we have grouped all the non- Yukawa terms in- 
volving neutralinos and charginos >CMSSM,inos- The scalar sector is implemented in terms 
of gauge-eigenstates while the fermion and gaugino/higgsino sectors are directly imple- 
mented in terms of mass eigenstates. For clarity and for generalization purposes, the 
FeynRules implementation is splitted into 18 files. 
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• Gauge . f r: the definition of tlie gauge groups and the indices, 

• PrmExt . f r: the external parameters, following the SLHA conventions, 

• Prmlnt . f r: all the internal parameters, derived from the external ones, 

• PrmAux.fr: auxiliary parameters, such as identity matrices, 

• FldFer.fr: physical fermionic states (quarks, leptons, charginos, neutralinos, gluino), 

• FldVec.fr: physical gauge bosons (gluon, photon, W and Z bosons), 

• FldSca.fr: physical scalar fields (sfermions and Higgses), 

• FldAux.fr: gauge-eigenstates (sfermions, Higgses, SU{2)l x U{1)y bosons), 

• FldGst.fr: gluonic ghost (required in order to have CalcHep running properly), 

• LagGau.fr: /^Cauge, defined in Eq. (^), 

• LagFer.fr: ^Permions, defined in Eq. (U), 

• LagSca.fr: £scalars, defined in Eq. ( |5.27] ), 

• LagFer.fr: £inos) defined in Eqs. ( p.35 ) and (5.37), 

• LagFD.fr: the definition of the superpotential, 

• LagFDgen.fr: the D-terms and a generic routine deriving automatically the F-terms 
from the superpotential, allowing for the calculation of /^fd, 

• LagBrk.fr: the soft supersymmetry-breaking Lagrangian £soft, from Eq. ( [5.38| ), 

• LagGst.fr: the ghost Lagrangian (only QCD ghosts are supported so far), 

• SUSY.fr: the main file collecting all the different pieces. 

Before running the model, the user has to provide values for the four switches FeynmanGauge, 
$CKMDiag, $sWScale and $svevScale. 

• The switch FeynmanGauge is not used in the present implementation and has to be set 
to False, since the whole model is implemented in unitary gauge. However, further 
developments will allow for various gauges. 

• The switch $CKMDiag allows for a CKM matrix different from the identity or not, 
depending on its True or False value. The CKM matrix could also be forced to the 
identity through a restriction file. 

• The two switches $sWScale and $svevScale can be set to the values "weak" or 
" susy " , regarding the scale at which the electroweak parameters and the vevs of the 
Higgs fields will be evaluated. 
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Figure 3: A schematic "moose" diagram of the Three-Site model. The circles represent gauge 
groups. The two circles on the left are SU{2) gauge groups while the one on the right is a U{1) 
gauge group. 



Possible extensions 

A large number of super symmetric models can be built from the most general MSSM. 
As stated above, we can take the proper limit in the various mixing matrices and get back 
to the commonly studied constrained MSSM scenarios. Moreover, other gauges are planned 
to be implemented {e.g., Feynman gauge or non-linear gauges). Other possible extensions 
consist in the addition of new particles, such as an additional singlet (the Next-to-Minimal 
Supersymmetric Standard Model) or right-handed neutrinos (and the corresponding sneu- 
trinos). 



5.4 The Minimal Higgsless Model 



To date, we have not detected a single fundamental scalar field and we do not know 
if we ever will. Electroweak symmetry breaking may instead occur as the result of the 
condensation of a strongly coupled bound state as in technicolor theories [140, 141] or it 
may occur as the result of boundary conditions in a compactified extra dimension as in 
so-called Higgsless theories [|142| , |143(1 . These theories are closely related [144] and each 



contains a tower of vector resonances which are responsible for unitarizing WW and WZ 
scattering in the absence of a fundamental scalar field (such as the Higgs) |145| , 146 , 147 ]. 
Unitarity constrains the mass of the lowest of these resonances to be below ~ 1.2 TeV, 



making it discoverable (or excludable) at the LHC ]|148| , |149| , |150(] . 

Although, it is impossible to implement the entire tower of resonances, it is also unnec- 
essary as the low energy phenomenology is dominated by the lowest modes. Deconstruction 



] 151 , 152 ] gives a consistent, gauge invariant way of implementing only a subset of the res- 
onances in a low energy effective theory. The minimal deconstructed Higgsless model 
contains all the non-scalar fields of the Standard Model plus the first new resonances, the 
W' , Z' and heavy partners of the fermions and is called the Minimal Higgsless Model or 
Three-Site Model |153| ]. A schematic "moose" diagram of the Minimal Higgsless Model is 
presented in Fig. ^ 



Precision electroweak constraints [154, 155, 156, 157] are satisfied in Higgsless models 
by allowing the fermions to delocalize into the bulk in a certain "ideal" way ]1l58|,[l59|,|T60|]. 
The consequence of this is that the heavy partners of the W and Z, the W' and Z', 
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are fermiophobic and have very small or vanishing couplings to the light SM fermions. 
Nevertheless, these heavy vector resonances can be discovered (or excluded) at the LHC 
|14^, [T49|, [150|. 

5.4.1 Model description 
Gauge Sector 

The gauge group of the Three-Site Model is 

G = SU{'i)QCD X SU{2)q X SU{2)i x [7(1)2, (5.55) 

where SU{2)q is represented by the leftmost circle in Fig. ^ and has coupling g, SU{2)i is 
represented by the center circle in Fig. ^ and has coupling g and U{1)2 is represented by 
the rightmost dashed circle in Fig. ^ and has coupling g' . We define 

X = - and t = — = (5.56) 
9 9 c 

where + (? = 1. 

The Lagrangian for the Minimal Higgsless Model can be written as a sum of six parts, 

C-MEM = -^^MHM, Gauge + -^^MHM, DS + -^MHM, GF + -^MHM, Ghost + -^^MHM, ip + -^^MHM, -02 • 

(5.57) 

The kinetic and self interaction terms for the gauge bosons is given by the usual gauge 
invariant terms: 

^MHM.Gauge = - ^ G^^ - - Fq^^^ Fj^^^ - - F^^^ F^ - -^F^^^ F^" . (5.58) 

The horizontal bars in Fig. |^ represent non-linear sigma models Sj which come from 
unspecified physics at a higher scale and which give mass to the six gauge bosons other 
than the photon. This is encoded in the leading order effective Lagrangian term 



f 



(5.60) 



[{D^J^oV Di'^o + (Df^^iY D^'^l , (5.59) 

where 

D^Eo = 9^So + igWo^'Eo - ig^oWi^f,, 

The non-linear sigma models can be written in exponential form 

= e*2-^/^, (5.61) 

which exposes the Goldstone bosons that become the longitudinal components of the mas- 
sive gauge bosons. tTj and Wj are written in matrix form and are 

2 3 J ^ 2 J ^ and ^2 = , 5.62 
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where j is or 1. 

The mass matrices of the gauge bosons can be obtained by going to unitary gauge 
(Ej 1) and are, 



M| = 



I -X 



-X 2 



and Ml = ^ 



( x^ -X ^ 



v 



—X 2 —xt 
-xt xH"^ 



(5.63) 



for the charged and neutral gauge bosons respectively where 



(5.64) 



and the photon is massless. After diagonalizing the gauge boson mass matrices, we find 
that the other masses are given by 



Mw 
Mw' 

for the charged gauge bosons and 

Mz -- 
Mz' -- 



Mg 



Mg 



2 + a;2 _ ^4 + x^, 



\/2 + a;2 + a/4 + a;-^, 



(5.65) 



Mg 
2 

Mg 



y/2 + X^{l+t^) - A, 

^2 + x^{l + t^)+A, 



(5.66) 



where 



A = ^y4 + x\l-t'^)^, 



(5.67) 



for the neutral gauge bosons. We note that x can be obtained from the ratio of the charged 
gauge boson masses 

^2 _ ( MwV _ 2 + x^- yiT^ 

^~\Mw'J ~2 + x2 + V4+^' 



R 



(5.68) 



The couplings can be determined in terms of the electric charge e, x and t. 



1 _ 1 1 1 
e 9 9 9 

1 

9" 



( 1 + + ^ ) , 



(5.69) 



x^ X 



2*2 



{g'f =e''{l + f + xH^). 
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Gauge Fixing Sector 

As we mentioned previously, tiie liorizontal lines in Fig. ^ represent non-linear sigma fields. 
Although tree level calculations can be done in unitary gauge, there are times when a 
different gauge is useful. Many calculations with gauge bosons in the external states can 
be computed more simply using the equivalence theorem and replacing the massive gauge 
bosons with the Goldstone bosons that they eat. Another case where a gauge different 
from unitary gauge is advantageous is in CalcHep, where the time of computation of 
processes is dramatically decreased when using Feynman gauge. For this reason, we have 
implemented this model in both Feynman and unitary gauges. 

We must first determine the Goldstone bosons that are eaten by the gauge bosons. We 
do this using the Lagrangian of Eq. ( ^.59 ). Expanding the non-linear sigma field, we obtain 



the mixing terms between the gauge bosons and the Goldstone bosons. After inserting the 
eigenwave functions of these fields, we obtain 

>CMHM,7riy = 7;9f(vlixViv - Vw) + vl{vw - 5xtVw)) {dfj.TT, W^} 

f ^ ' (5.70) 

+ -g^(fi,{xv%^,-vw')^v\,{v\^,-bxtvl^,)) {d^Ti^W'^] 

where 5 is 1 if the gauge boson is neutral but otherwise. 

The gauge fixing function is constructed to fix the gauge and cancel the mixing of the 
Goldstone bosons and gauge bosons. For each site, the gauge-fixing term is 



(5.71) 
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= d 
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^ J f ( „ns\ 
1^9 J {-T^ I ), 



where by tt"'^ we mean just the neutral sector of vri, namely 




(5.72) 



^MHM,GF = --Tr^Gg + GI + G2) . (5.73) 



With this definition, the gauge fixing Lagrangian is 

1 

Ghost Sector 

The ghost Lagrangian terms are obtained by multiplying the BRST transformation of the 
gauge fixing term on the left with the antighost. To do this, we must find the BRST 
transformations of the gauge fixing terms. We begin by writing the infinitesimal BRST 
transformation of the fields in the gauge fixing term. 



^BRST Wf^j = -{df,Cj + igj [Wf,j ,Cj]), (5.74) 
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for the gauge bosons where cj is the ghost for site j in matrix notation 



\ 2 3 V2 3 \ and C2 = P , 5.75 

where j is or 1. The BRST transformation to quadratic order in the Goldstone bosons is 



1 i 1 1 r 

(5.76) 



so that 



Go 


= d 


5 Wn 




= d 




G2 


= d 


^brstW2 



l^afi^RST'^O ), 

l^afi^BRST'^l - '^fliiST^o), (5-77) 



"bast"! y 



The ghost Lagrangian is 

-Cmhm, Ghost = — Tr^co^B^gyGo + ci^^^^^Gi + C2Sgj^gj.G2^ + h.c. (5.78) 
Fermion Sector 

The vertical hues in Fig. ^ represent the fermionic fields in the theory. The vertical lines on 
the bottom of the circles represent the left-handed chiral fermions while the vertical lines 
attached to the tops of the circles are the right-handed chiral fermions. Each fermion is a 
fundamental representation of the gauge group to which it is attached and a singlet under 
all the other gauge groups except U{1)2- The charges under C/(l)2 are as follows: If the 
fermion is attached to an SU{2) then its charge is 1/6 for quarks and —1/2 for leptons. If 
the fermion is attached to U{1)2 its charge is the same as its electromagnetic charge: for 
neutrinos, —1 for charged leptons, 2/3 for up type quarks and —1/3 for down type quarks. 
The usual gauge invariant kinetic terms are used, 

^MHM.V' = ii'LQJp^LO + i'ipLllplpLl + ilpRlplpRl + i4'R2ptpR2- (5.79) 

The fermions attached to the internal site {SU {2)i) are vectorially coupled and Dirac 
masses are thus allowed. We have taken these masses to be Mp. The symmetries also allow 
various linkings of fermions via the non-linear sigma fields. We have assumed a very simple 
form, inspired by an extra dimension and represented by the diagonal lines in Fig. ^. The 
left chiral field at site j is linked to the right chiral field at site j + 1 through the non-linear 
sigma field at link j. The mass parameter for these diagonal links is taken to be eiMp and 
€rMf for the left and right links respectively. All together, the masses of the fermions and 
the leading order interactions of the fermions and non-linear sigma fields are given by 

-CmHM.VS = -Mf EL^LoSo^iil + 4>L11pRl + i>LieRT.i'tpR2 , (5.80) 
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where e/, is the same for ah fermions but e/j is a diagonal matrix which distinguishes flavors. 
For example, for the top and bottom quarks we have 




= r-^ " . (5.81) 



The mass matrix can be obtained by going to unitary gauge and diagonalized by a 
biunitary transformation. By doing this, we find the following masses. 



7f 



Mfo =^-^\/l + el + el-C, 



where 



(5.82) 



C=^{l + el + eir-4elel. (5.83) 
5.4.2 FeynRules implementation 

The FeynRules implementation was initially based on a LanHep implementation [148| . 



It was translated to FeynRules syntax and slightly modified to fit the requirements of 
the FeynRules package and interfaces. All symmetries, fields and parameters were imple- 
mented according to the definitions of the last subsection. The independent variables that 
the user can adjust are the electromagnetic and strong couplings, the masses of the Z, W, 
W' and SM fermions (where not set explicitly to zero) and the scale of the heavy fermions 
Mp. The scale of the Z pole mass and Fermi constant are also implemented as required by 
some Monte Carlo codes, but they are not used directly in this model. A change in these 
last two parameters will not affect the value of the other parameters. 

Two gauges were implemented in this model file. A variable FeynmanGauge, similar 
to the SM case, was created to switch between the two. When the switch is set to True, 
Feynman gauge is chosen and the Lagrangians contain the Goldstone bosons eaten by the 
gauge bosons and the ghosts. If, on the other hand, it is set to False, all Goldstone and 
ghost terms are set to zero. 

5.5 Extra Dimensional models 

One popular approach to solve the hierarchy problem of the Standard Model is to ex- 
tend space-time to higher dimensions |161| , |162| ]. In this framework, the usual four- 
dimensional space is contained in a four-dimensional brane embedded in a larger structure 
with additional dimensions, the bulk. Moreover, gravitational and gauge interactions 
unify close to the only fundamental scale of the theory, the weak scale. In theories with 
Large Extra Dimensions (LED) [ p7| , |8^ , |163| |, the gravitational interactions are the only 
ones propagating into the bulk, which dilutes their coupling strength and make it appear 
weaker inside the four-dimensional branes. As a consequence, the graviton field is ac- 
companied by a tower of massive Kaluza-Klein states. In scenarios with Universal Extra 



Dimensions |161], each field of the Standard Model possesses a tower of excitations with 
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the same quantum numbers, but different masses. Even if none of these states has been ob- 
served so far, TeV-range excitations could be detected at the present Tevatron or the future 



Large Hadron Colhder. In a minimal Universal Extra Dimension scenario (MUED) [164| , 
one has one single flat additional dimension, y, which is spatial and compactified on a 
S1/Z2 orbifold of radius R. Momentum conservation in the extended space-time gener- 
ates a conserved quantum number, the Kaluza-Klein parity, which implies that different 
Kaluza-Klein modes cannot mix and that the lightest excitation could be a candidate for 



dark matter |165, 166 , 167 1 



5.5.1 Model description 

Large Extra Dimension Model In a LED theory with one compact and space-like 
additional dimension, all the Standard Model fields are confined into a four-dimensional 
brane and there are only Kaluza-Klein excited states for the graviton. The Kaluza-Klein 
states have higher masses than the standard graviton but both behave the same, i.e., they 
couple gravitationally to the Standard Model fields. We start by specifying the generic 



effective Lagrangian of an unbroken gauge theory |168] 



Gauge Fermions Scalars • 

The Lagrangian of each sector can be expressed as 



n V 



3(n + 2) 



(n) rpi fj, 



(5. 



(5.85) 



K = V16vrGiv with Gn being the four-dimensional Newton constant, h^'^ the n-th gravi- 
ton Kaluza-Klein mode in four dimensions and c/)^"^ its scalar component in the fifth di- 
mension. The graviton couples to the energy momentum tensor, 



5C 



6g' 



(5.86) 



For a generic unbroken gauge theory, the various energy-momentum tensors T^j^ for an 
unbroken gauge theory read^. 



T: 



Scalars 



pFermions 



+ 



(5.87) 



rpGsMgC 



4 '^ 2 ^ 



+ [-FPF^p + m\A^A, 



where $ is a (complex) scalar, ^' a fermion, and A^ a vector field. and F^jy denote 
the usual covariant derivatives and field strength tensors. Prom this general model, we can 



*We work in unitary gauge. 



-50- 



derive a realistic Large Extra Dimensional theory containing all the Standard Model fields. 
We choose to work in the unitary gauge also for the gravitational part, which eliminates 
the non-physical degrees of freedom absorbed by the massive fields, and we re-write the 
Lagrangian as, 

/:LED., = -f J^fa^'^^-^T/,), (5.88) 



where ^^^("^ is the n-th Kaluza-Klein mode of the physical graviton and j now denotes 
explicitly the Standard Model sectors. The energy- momentum tensor of the Standard 
Model can be written as a sum of four parts, 



rp ^^Fermions _i_ ^riHiggs _i_ rpGaugc , rri Yukawa 



(5.89) 



The energy-momentum tensor of the fermionic sector reads, 

fTiFormions 

= -Vp-u [{Qi i YDp Qi + Li i YDp Li + um i YDp um + dm i YDp dm + Im i YDp Im] 
- \dp {Qi iYQi + Lii Y Li + um i Y um + dm i Y dm + Im i Y ^Ri) 

■7^ {Qi i If^Diy Qi + Lii "fp.Du Li + um i l/iD^ um + dm i l^Dv dm + Im i 7fiDu Im) 



+ 



--{dp. {Qi i^uQi + Lii Li + Um i iv um + dm i iv dm + Im i Iv Im) + (a* ^ i^) 



(5.90) 



where Z)^ denotes the SU{3)c x SU{2)l x U{1)y covariant derivative. Similarly, the 
energy-momentum tensor for the gauge sector reads. 



jiGauge 



pv 



-r]pu 



BPB^p-W^pPW'^p-GlPGlp. 



(5.91) 



Finally, the energy-momentum tensors for the sectors involving a Higgs boson read. 



jiHiggs 
pu 



n Yukawa 



-Vpv 

with l> = ia"^ . 



Dp^^DP^ + ^2$t$ _ A ( 



to 






+ 



-UiR Vij Qj ^ - dm yf • Qj - km y\. Lj $ + h.c. 



(5.92) 



Minimal Universal Extra Dimension Model 

We consider a theory in five dimensions, the fifth one being spatial and compactified 
on a S1/Z2 orbifold of radius R, i.e., the points y and —y are identified, where y is the extra 
coordinate. This symmetry is essential to define chiral fermions. Unlike LED models, in 
UED models all fields have access to the extra dimension and depend on the fifth coordinate 



y. We split the most general Lagrangian in four pieces 161 , 169 1, 

-^MUED = -^MUED, Gauge + -^MUED, Fermions + -^MUED, Higgs + -^MUED, Yukawa- (5.93) 
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The gauge sector is described by the field strength tensor terms 

^MUED,Gauge = --BmnB^^ - -^MN^k^^ - ^'^MnG^^ ■ (5-94) 

We define the subscript M = (/x, 5), where /x is the usual four-dimensional Lorentz index 
and 5 the fifth dimension index. The fermionic sector is decomposed into its leptonic and 
quark part, 

.^MUED.Fermions = .^MUED, Leptons + .^MUED, Quarks, (5.95) 

with 

>CmUED, Leptons =iL(x'^, y)r*^DM^(x'^, 2/) + y)r^I?M^(.T^, 2/), 

>Cmued, Quacks =iO(x^ y)T^DMQ{x^', y) + iU{x^, y)T^DMU{x^, y)+ (5.96) 
iD{x^,y)T^DMD{x^,y). 

Q (x^, y) and L {x^^, y) denote SU{2)l fermion doublets and U {x^^, y), D {x^, y), E {x^^, y) 
are the up-type quark, down-type quark and charged lepton singlet fields. The gamma 
matrices in five dimensions are defined as 

r^^ = {^/^,i^5|^ (5.97) 

while the five-dimensional covariant derivative Dm is 

DM = dM-i Ygf^Bu - \i - i gf^T^Glj. (5.98) 

The hypcrcharge Y , the Pauli matrices and the color matrices are the generators 
of the gauge groups, while the five- dimensional coupling constants are related to the four- 
dimensional ones through 

gf = ^g,. (5.99) 

Finally, the Higgs Lagrangian is given by 

-CMUED.Higgs = [i^M^(:r^y)]t [Z)^i7(x^ y)] + M'i^^(x^ y) /^(a;^ y) 

2 (5.100) 



A 



H^x^,y)H{x^',y) 



and the Yukawa sector describing the interactions between the fermions and the Higgs field 

by 

>CMUED,Yukawa = - A„Q(.T^y)[/(.x^y)i^(,x^y) -AdQ(x^y)L>(x'^,y)i^(x^y) 

(5.101 

-XeQ{x'',y)E{x'',y)H{x'',y), 
where H = ia^ H*{x^, y), and Aj are the usual Yukawa matrices. 
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5.5.2 FeynRules implementation 

Large Extra Dimension Model 

The FeynRules implementation for the LED model described is based on the La- 
grangian defined in Eq. ( |5.84| ). The theory contains, besides the graviton, the full set of 
SM fields, together with their respective coupling to the graviton field via the stress tensors. 
For the graviton we restricted the implementation to the lowest mode, i.e., only the mass- 
less graviton is taken into account. An extension of this model to include Kaluza-Klein 
excitations of the graviton, as well as to any BSM extension of the SM, e.g., the HAH 
model described in Section ^, is straightforward. 



Minimal Universal Extra Dimension Model 

To implement the MUED model in FeynRules, we start with the most general five- 
dimensional Lagrangian described above. Then, FeynRules derives the effective four- 
dimensional Lagrangian and the corresponding Feynman rules automatically. This is 
achieved by expanding the five-dimensional fields and imposing the dimensional reduction 



by integrating out the extra coordinate y |170|| , 
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+V2Y^ [PRi^Pix] cos (^) + PLi^^^\x) sin ^ ^ 

> n=l 




(5.102) 



where A5{x^,y)), H{x^,y), ^{x'^,y) and '4j{x'^,y) denote a five-dimensional 

gauge field, a Higgs field, a fermionic doublet and a fermionic singlet field, respectively. 
To integrate out the extra dimensional coordinate, we follow the integration procedure 
described in Refs. pll, p|, p70| , using orthogonality relations. Let us note that if one 
would like to add other pieces to the Lagrangian, it is sufficient to deal with five-dimensional 
expressions, together with the appropriate definitions of the field expansions. 

In our implementation we are considering Kaluza-Klein excitations only up to the first 
mode. The inclusion of the next Kaluza-Klein states is of course straightforward, each 
considered Kaluza-Klein mode has to be defined in the model file. For a given field, we 
identify the zeroth mode as the Standard Model particle and define as new particles the 
Kaluza-Klein excitations. The zeroth mode particle content is thus identical to the SM 
one, see Table while the particle content of the first modes is summarized in Table |^. 



-53- 



Particle description 


Spin 


Representations 


Lepton doublet 


Lf' = 




1/2 


(1,2,-1/2) 


Lepton singlet 






1/2 


(1,1,-1) 


Quark doublet 


^ = 


I ' \ 


1/2 


(3,2,1/6) 


Up-quark singlet 






1/2 


(3,1,2/3) 


Down-quark singlet 






1/2 


(3,1,-1/3) 








1 


(1,1,0) 


Gauge bosons 






1 


(1,3,0) 








1 


(8,1,0) 


Higgs 









(1,2,1/2) 



Table 6: A summary of the first mode Kaluza-Klein particles. 



As for the SM particles, the new gauge-eigenstates mix to physical states 



W, 



1(1) 



2(1) 



^(1) 



T^3(l) 



(5.103) 



The first mode analog of the electroweak mixing angle is assumed to be very small. There- 



fore the first modes of the Z boson and the photon are simply W^i 
tree-level mass M^^) of the first gauge boson excitations are. 



3(1) 



and B, 



(1) 



The 



M« = V^ + (m(0))^ 



(5.104) 



where m^^^ is the zeroth mode mass. In the fermion sector, there might be a mixing between 
SU{2)l doublets and singlets. However, this mixing is expected to be highly suppressed. 
Therefore we ignore it in our implementation. As in the SM, the zeroth mode fcrmions get 
their mass from the Yukawa couplings after the Higgs get its vev. For the Kaluza-Klein 
first modes, the tree-level mass terms follow the structure 



^ _pW fin) ^ 




where F^^^ and denote the doublet and the singlet fields, respectively. The diagonal 
contributions come from the kinetic terms in the fifth dimension and the off-diagonal ones 
are induced by the Higgs vev. After the diagonalization of the mass matrices, the masses 
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of the physical eigenstates are obtained, 



( 



n 



) 



2 



(0) 

7 
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m: 



R 



+ m 



(5.105) 



Possible extensions 

As we have mentioned before, we are only considering expansions of the fields up to the 
first Kaluza-Klein excited states. However one could easily incorporate the next modes 
by expanding the fields up to the n*^ component. It is straightforward to incorporate the 
second mode Kaluza-Klein excitations, by following the example of the first modes. 

5.6 An access to the low energy world 

Effective theories allow us to have predictions about the masses and the interactions at low 
energy without having to know everything about the fundamental theory. We only need 
to know the (exact or approximate) symmetries of the fundamental theory. One of the 
most successful examples is Chiral Perturbation Theory (xPT), the effective theory of the 
strong interactions. Most of the processes involving pseudoscalar mesons happen at scales 
at which pertubative QCD cannot be trusted anymore. Despite our resulting incapacity 
to compute their properties directly from the fundamental Lagrangian, the shape of the 
effective Lagrangian can be inferred using the global approximate chiral symmetry of QCD. 
The other famous example, especially in this pre-LHC era, are all the effective models 
developed to solve the hierarchy problem where the Higgs boson is a pseudo-Goldstone 
boson of a new strong sector. These models were built on the fundaments laid down by 
QCD. However, in this case, much less is known, neither the fundamental nor all the low 
energy degrees of freedom. This is probably the origin of the number of models available 
going from Technicolor to Extra Dimensions. However the low energy effects of many of 
them, called Strongly Interacting Light Higgs (SILH) Models, can be described by the same 
effective Lagrangian up to some coefficients |171| ]. 

5.6.1 Model description 
xPT at the lowest order 

The leading effective Lagrangian invariant under the SU{nF)L x SU{nF)R x U{l)v-, where 
ni? is the number of massless quark flavors, is [|172|] 



where U transforms as [/ — > giUg]^. The first part of the Lagrangian is due to the axial 
[/(I) breaking required by the rj' mass |136| ] and is at the source of the strong CP problem. 
The second term is the usual non-linear sigma Lagrangian and the last one takes into 
account the quark masses. Only the first three flavors are sufficiently light compared to the 




(5.106) 
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confinement scale of QCD of about 1 GeV to consider the chiral symmetry as approximate, 
so we set ni? = 3 in the SM. At low energy, the symmetry of the strong interactions has to 
be broken spontaneously to its vectorial subgroup U{nF)v in the limit where the number 
of colors (Nc) is large. The unitary matrix U can thus be developed as a function of the 
pseudo-Goldstone bosons around its vacuum [173|, 



U 



fc=i ^ •'^ 



(5.107) 



with 



TT 



\ V2K- 



3'/ 



V2: 



IT ' 



-vr3 + 

V2K0 



3'/ 



(5.108) 



The coefficients are partially fixed by unitarity. 
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cH \ , 

2 2 16 



(5.109) 



The remaining free parameters, b, c, . . . , can be used to check the computation of any phys- 
ical quantity which should be independent of these quantities [174, 175], or they can be 
fixed to obtain the most suited form for U. For example, the common form U = exp ^ ^^^ ^ 

requires b = ^ and c = j^. The mass matrix of the three neutral states, vr^, rj^ and r]^ is 
not diagonal. Consequently, the physical eigenstates are obtained by a rotation. 



The SILH Model 

The Strongly- Interacting Light Higgs Model is an effective theory of a possible strong sector 
responsible for the electroweak symmetry breaking (EWSB). The aim of this model is to 
disentangle if the EWSB is due to a strong sector or not by testing the interactions of the 
Higgs and the gauge bosons. 

The SILH requirements are that the new strong sector has a spontaneously broken 
symmetry at low energy and that the Higgs boson is an exact Goldstone boson when the 
SM interactions are switched off. As for the mass matrix of the light quarks in QCD, the 
SM interactions are included order by order as small breaking parameters and induce a 



-56- 



mass term for the Higgs boson. The httle Higgses and Holographic composite Higgs models 
are examples of such theories. 

At scales lower than the resonances, the effects of the model are described by a set of 
dimension-six operators involving the Higgs boson, 



-H^HhHfR + h.c.^ 



/2 2 2 2 

+ + (5.110) 



where H'^ D = H'^D^H — [D^H"^)!!, Qp is the coupling of the new strong sector, Mp 
is the mass of the new heavy states and / = Mp/ gp. For completeness, the model also 
contains the dimension-six Lagrangian for the SM vectors, 

^vect = -^^^{D'Wp,r{Dpwn' - ^i(9'Bp,)(9pBn - ^^.{D^'Gp^nDpcn" 



+ii^^^^-^^^'^'^^"'' + iM^/^^^^r^G^^- • (5-111) 



{f'^''^^^ ^-P^ + 16^2 jv,^ 

The Cj are free parameters of the model and can thus be fixed to reproduce the low energy 
behavior of a more specific model. 



5.6.2 FeynRules implementation 

xPT in FeynRules 

xPT is implemented in FeynRules at the lowest order ( 5.106| ), namely at the order . In 
the implemented model, the Lagrangian is developed up to the ^ (vr^) ^ with arbitrary co- 
efficients h and c. Consequently, vertices depend linearly on scalar products of the momenta 
and contain up to six scalars. The Lorentz structure of these vertices are not included in 
the default generic FeynArts file which contains only renormalizable structures. Let us 
note that at present this implementation only works with the FeynArts interface. 

Since the isospin limit = rrid = m is taken in the FeynRules model file, only r/g 
and rjQ mix. The mass eigenstates are given by 







' COS 9 


— sin0 


II 


( 7?M 






i sin0 
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vr vr 
'4' +4 



(5.112) 



The isospin breaking can be easily added. The major modification is to extend the 2x2 
matrix in ( |5.112 ) to a 3 x 3 mixing matrix since the three neutral states mix in this case. 

®A11 the higher order terms have been removed from the Lagrangian to save computation time. 
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However, the effects are small and thus usually neglected. For the same reason, the mass 
matrix of the quarks is assumed to be real, but, the phase allowed by the U{1)a breaking 
can easily be added in the definition of the quark mass matrix. 



The lowest order Lagrangian, Eq. (5.106), can reproduce the experimental data within 



70. These discrepancies can be solved with next order corrections. The inclusion of the 
next order operators either (p^,0) |176] or {jp', l/Nc) 17^, |17| in the La grangian is 



also straightforward since U is already defined in the model file. In the same way, the weak 
and electromagnetic interactions could also be added. However, the Lorentz structure of 
the vertices is hardcoded in the FeynArts interface, so the new structures generated by 
these new operators need to be added in the interface and in an associated FeynArts 
generic file, as it was already done for the lowest order Lagrangian. 

SILH in FeynRules 

The SILH implementation is based on the SM implementation but restricted to the unitary 
gauge. The two Lagrangians of Eqs. ( ^.110| ) and ( 5.111 ) have been added to the SM one 



with all the q considered as free external parameters. The first one contributes to the 
kinetic term of the Higgs and of the gauge bosons, so the physical fields are renormalized 
versions of the bare fields appearing in the Lagrangian^'^ , 

H = Hbare (^1 + ' 



i,bai~e y 



M2 o2\ (5.113) 



The gauge couplings have also to be redefined to obtain canonical kinetic terms for the 
gauge bosons. 



g = gtare 1 + C^y 



¥w_, 9" g' 



g' = g'bare U + Cfi tan^ 6^-^ - ^c^TT^^ , (5.114) 



gS = gs,bare 1 + ^Cgj—yT— , . 

Finally, the vev and the masses of the Higgs boson and the fermions are also corrected 



''(cw + cb)M^ /Mp is supposed to vanish to avoid non-diagonal kinetic terms. 
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compared to their SM values, 



2 1^ 



1 - 



1 - C ( C// + |C6 



(5.115) 



All of these redefinitions have been done at the first order in ^ = jj. Any result is thus 
valid only at this order. It should be noted that all of these non-renormalizable interactions 
cannot be transferred directly to any other HEP tool for the moment. 



6. Validation of the implemented models 

In this section, wc first review some generic features of the various validation procedures 
used to assess the model robustness and then move to the results of the validation process 
for the Beyond the Standard Model theories presented in the previous section. 

6.1 Strategy for the validation of the implemented models 

In order to validate the implementation of a model inside FeynRules, a very first natural 
check is to compare the obtained Feynman rules, using directly the Mathematica output, 
against the ones found in the literature. Subsequently, using the existing interfaces to 
symbolic tools such as FeynArts, we can go further in the validation procedure with 
analytical checks of some observables, confronting the results obtained with the help of 
the model files generated by FeynRules to the corresponding expressions found in the 
literature. 

The next step in our validation procedure regards the model files generated through 
the different interfaces between FeynRules and matrix-clement generators. For a given 
model, we calculate predictions using the set of Monte Carlo tools interfaced to Feyn- 
Rules both with the model files generated by FeynRules and with the built-in (stock) 
implementations of the considered model, if they exist of course. After consistently fixing 
the set of external parameters to the values corresponding to a chosen benchmark point, 
results for various quantities, such as decay widths, total cross sections or (unintegrated) 
squared matrix elements at a given phase-space point, are computed and confronted. For 
a phase-space integrated observable a, we evaluate all the possible quantities 

'^-^ = '^ \a.J ^ (6-1) 

where a""'^ refers to predictions obtained with two specific generator and a given model file. 
The Aaft's quantify possible discrepancies between different implementations. For example, 
^MG-FR, CH-ST represents the discrepancies between predictions obtained with CalcHep, 
using the built-in implementation of the considered model (CH-ST), and with MadGraph, 
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using the FEYNRuLES-generated model files (MG-FR). For unintegrated squared matrix 
elements |Mp, we generalize this quantity to 



A'^'^ = max < 



a,b 



phase space points 



^a\2 _ 1^612 



|M«|2 + IM'-P 



(6.2) 



summing over a given number of phase-space points. 

For each model presented in the previous section, the input parameters used in the 
comparison and some numerical examples are given in Appendix ^ The complete list of 
results can be found on the FeynRules website [^]. 

6.2 The Standard Model 

All the Feynman rules obtained with FeynRules were checked against the expressions 
given in the literature, while the total cross sections for a set of 35 key processes have 
been evaluated in CalcHep, MadGraph/MadEvent and Sherpa and compared to 
the existing stock versions. A selection of processes, together with the set of external 



parameters used for this check, is given in Appendix B.l. Note in particular that Mad- 
Graph/MadEvent and Sherpa work in unitary gauge, whereas CalcHep allows for 
both unitary and Feynman gauges. We therefore also demonstrated explicitly the gauge 
invariance of our implementation. 

6.3 The general Two-Higgs-Doublet Model 

The validation of the present implementation is done to by comparing 2 — > 2 matrix 
elements, for 10 different phase-space points, between the FeynRules unitary-gauge im- 
plementations in MadGraph/MadEvent, SHERPAand CalcHep, and the existing stock 
implementation in MadGraph/MadEvent described in Ref. |^] and distributed with the 
code. Since all the vertices appearing in the 2HDM are also present in various BSM models 
presented in this work, such as the MSSM, this limited comparison is clearly sufficient to 
fully validate our implementation. 

The selected benchmark point includes non-trivial values for all the Aj parameters 



of Eq. (5.12), as well as Yukawa couplings for the second and third generation fermions 
leading to FCNC effects. Agreement, claimed if the quantity A^'^ defined in Eq. ( |6.2| ) is 
smaller than 0.1%, has been found for all the tested processes involving various combina- 
tions of scalars, fermions and gauge bosons both in the initial and final states. Further- 
more, a selection of 185 2 — > 2 cross sections have been computed using CalcHep, Mad- 
GRAPH/MADEvENTand Sherpa, and compared to the existing MadGraph/MadEvent 
implementation. In all cases the results agree within 1% between the different codes. Some 



examples are shown in Appendix B^. In addition to these checks, the behavior of various 
cross sections at high energies for processes involving new scalars has also been verified. In 
each case, the cross section behavior is in agreement with unitarity expectations. 

6.4 The most general Minimal Supersymmetric Standard Model 

We have compared the Feynman rules computed by FeynRules to those which can be 
found in the literature, both for the general MSSM |179|, |180|] and for a constraint MSSM 
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where all the scalar mixings are neglected |101, 181], and we have found agreement for all 
the vertices. Then, we have re-calculated all tree-level squark and gaugino hadroproduction 
helicity amplitudes in the case of general (and possibly complex) scalar mixing with the 
help of FeynArts/FormCalc and the model file generated by FeynRules. The results 



have been compared to the analytical formulas given in Refs. |111, 112| and we found 
complete agreement. 

To validate the FEYNRuLES-generated model files for the various Monte Carlo gener- 
ators, we compared the results obtained in the very particular limit where CP symmetry 
and flavor are assumed to be conserved within the whole model, the CKM matrix appearing 
in the charged-current interactions being thus neglected as well. In addition, in the scalar 
sector, the flavor-conserving helicity mixings are neglected for first and second generation 
sfermions. In this scenario, built-in implementations exist in MadGraph/MadEvent 
and CalcHep. As a benchmark scenario, we choose the typical minimal supergravity 
point SPS la [182|, defined by the input parameters given in Appendix B^. For Sherpa, 
the validation process is on-going and extensions to other programs that will be interfaced 
to FeynRules in the future are foreseen. Let us note also that, since CalcHep is not 
able to deal with files following the SLHA conventions, we had to modify the built-in files 
in order to correctly fix the free parameters of the model. 

We start by evaluating all the 320 two-body decay widths corresponding to kinemat- 
ically allowed decays in our scenario in order to check the norm of various three-point 
vertices. We find a complete agreement for both SM and MSSM processes. In particular, 
we have confronted the MG-FR implementation to the already validated MG-ST imple- 
mentation 1 131], and the agreement between the two predictions was evaluated through 
the quantity Amg-ST,mg-fr defined in Eq. (^), which has been found to be smaller than 
1%, the differences being hence associated to a pure Monte-Carlo statistical error. Some 
examples can be found in the Table 15 of Appendix B.5. 

Subsequently, we have investigated the total cross section for 636 key 2^2 pro- 
cesses with the help of both MadGraph/MadEvent and CalcHep, and with both 
FEYNRuLES-generated and built-in model files. In order to properly compare the different 
implementations and unitary cancellations, we manually set all the widths of the particles 
to zero. We checked that the six quantities A^b defined in Eq. (^) for any set of two 
predictions a and h are below the percent level. This check allows not only to verify the 
absolute value of all the three-point vertices and some four-point ones, but also to start 
being sensitive to the relative sign of a large part of those vertices. This is also the very first 
systematic comparison between predictions obtained with MadGraph/MadEvent and 
with CalcHep for the MSSM. We have considered the production of either a pair of any 
Standard Model particles, or a pair of any supersymmetric particles, from any Standard 
Model initial state. Even though most of these channels are not really phenomenologically 
relevant because most of the initial states are impossible to realize at a collider, they allow 
for a sensitive check of each three-point vertex, and a large part of the four-point ver- 
tices. Let us note that the remaining untested vertices, such as the four-scalar interactions 
in >C Scalar FDW ™ Eq- ( ^-33 ), are irrelevant for 2 — > 2 leading-order calculations with a 
Standard Model initial state. We have considered two different cases, one where we have 
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fixed the energies of the initial particles to 600 GeV, and one to 1 TeV. We found that 
the six possible quantities are below the percent level comparing any two of the four 
implementations, MG-FR, MG-ST, CH-FR and CH-ST, except for 

W^W- ^W^W- and ht^ ZW- (6.3) 

where the discrepancies are due to the presence of singularities in the matrix elements, 
which makes any of the considered Monte Carlo generators unable to correctly evaluate 
the total cross section, and except for about 15 processes where the built-in implementation 
of the MSSM in CalcHep gives results different from those of all the other implementations 
due to mistakes in the model files for several triple scalar couplings. For the two processes 
in Eq. ( |6.3| ), we have evaluated the (unintegrated) squared matrix elements |Mp at given 
phase-space points for all models and computed the quantity A^'^ of Eq. ( |6.2p . Summing 
over 100 different random phase-space points, we have found that A^"^ is below 1%, which 
is sufficient to conclude that the four implementations agree. Some examples of numerical 
results can be found in Tables |l^, ^ and 18 in Appendix |B.3| . Let us recall that the 



complete list is available on the FeynRules webpage. 

Finally, in order to test every sign of each three-point vertex and a large part of the 
four-point vertices, we have evaluated squared matrix elements for 100 random phase-space 
points and calculated the quantity A^'^ defined in Eq. (6.2). We have investigated 2708 



processes, relative to the production of two supersymmetric particle plus one Standard 
Model particle at a center-of-mass energy of 2 TeV. Using MG-FR and MG-ST, we have 
found a complete agreement for each process, i.e., was below 0.1 %. Ah the results 
can be found on the FeynRules website. The comparison with matrix elements calculated 
by CH-FR and CH-ST is devoted to a further study. 

6.5 The Minimal Higgsless Model 

Our FeynRules implementation was compared to the LanHep implementation, both in 
unitary and in Feynman gauge. It was run in Feynman gauge in CalcHep and CompHep 
and in unitary gauge in CalcHep, MadGraph/MadEvent and Sherpa. The parame- 



ters used for this validation and the final particles cuts are given in Appendix B^. In all 
cases, the cross section was calculated, compared, and agreement to better than 1% was 
found. 

Some examples of the obtained cross sections are given in the various tables of Ap- 



pendix B.4. In Table 20, a selection of strong processes are presented. In Tables 21 and 



22, charged and neutral electroweak processes are presented, respectively. 



6.6 Extra Dimensional Models 
Large Extra Dimension Model 

The Feynman rules for LED obtained by FeynRules were explicitly checked with those 
available in the literature. Unlike for the models presented above, we cannot use any 
matrix-element generator to compute cross sections or decay rates, because the interfaces 
linking FeynRules to Monte-Carlo generators are not yet defined to work for a theory 
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with spin- two particles. Therefore we have chosen to vahdate our LED model implementa- 
tion via analytical expressions. We found complete agreement both for the Feynman rules 



of the generic LED theory which can be found in Refs. |168, 183| and for those of the full 
LED implementation. Let us notice that from the generic LED model validated above, we 
can extrapolate the results of the latter to guess those of the QCD and electromagnetic 
part of the full LED model. This check was performed and conclusive. 

Universal Extra Dimension Model 

For the MUED model, we start by confronting the Feynman rules to those available in 
the literature. In a second step, we have calculated cross sections for 2^2 processes 
and compared the results obtained with the help of MadGraph, Sherpa and CalcHep, 
and we confronted our results to those obtained with the help of an existing CalcHep 



implementation |16£]. The input parameters for the benchmark point which we have 



chosen for our numerical validation are given in Appendix B.5. Let us note that the only 
Higgs field which we have considered is the Standard Model one, even though its first 
Kaluza-Klein excitation is also implemented. The comparison has been carried out for 



the 118 processes in total |165| , [166 , |167| , [184| . First, using MadGraph and considering 



Standard Model processes, we have calculated squared matrix elements at given phase- 
space points and confronted the results obtained with the FEYNRuLES-generated model 
file to those obtained with the built-in Standard Model implementation. Subsequently, 
using the FEYNRuLES-generated model file for MadGraph, Sherpa and CalcHep and 
the existing CalcHep implementation of the MUED model, we have compared total cross 
sections for the 118 chosen processes at a center-of-mass energy of 1400 GeV. We have 
found agreement for each of them, and some examples are shown in Tables 24 and 25 in 



Appendix B.5. 



6.7 Low-energy effective theories 
xPT at lowest order 

The xPT model was checked analytically. The expressions of the amplitudes computed with 
pen & paper work of Ref. [|185| ] have been compared to the one obtained using FeynRules 
and FeynArts. The loop amplitudes have been integrated using cutoff regularization and 
only the terms quadratic in the cutoff A have been computed. More precisely, we computed 
the one-loop corrections to the two-point functions n — tt, rj — rj, rj' — rj' and r] — rj' . The 
momentum independent part of each of these amplitudes, but the last one, corresponds to 
mass corrections and is 6-independent as it should. For example, the renormalization of 
the pion wave function is 

A^ 
(4vr/)^ 



Z = 6{l-8b) ^ ^ , where ttr = Vl + Ztt, (6.4) 



while its mass correction vanishes^^. The tj' rjinr decay amplitude has also been com- 
puted at tree-level and at one-loop. In this last result, also the c-dependence cancels. 



All the other results can be found in Ref. [185| 
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Eventually, a total of about 60 diagrams were computed with both methods and perfect 
agreement was found. 



SILH Model 

The check of the SILH model consists in an analytic comparison between the decay widths 
computed in Ref. [171] and computations based on the vertices given by FeynRules. We 
used the same simplifications, i.e., ct = and {cw + CB)Myy /M^ = and neglect £vect- 
For tree-level decay widths of the Higgs into two fermions, both implementations leads 
to exactly the same results. For decays into a gauge boson pair, the contribution to the 
Feynman rules from higher-dimensional operators read, e.g. for the hWW vertex. 



chw 



[n 



2v (4^ I'l Pi Pi Pi 



(6.5) 



where pi, p2 and ps denote the momenta of the Higgs boson and the two W bosons 
respectively. As a consequence, the corrections are not just proportional to the SM decay 
widths since the vertices have a more complicated kinematic structure. We find 



r (/i ^ W+W-) =r{h^ w+w 



+(.chw 



I SM 

inn 



1 



9p 



(47r)2 167rz;2 



« + 2M2,) (l-4^) 



1/2 



(6.6) 



r(/i ^ ZZ) = T{h^ ZZ) 



SM 



1 - ( c/f - ^ {cw + tan^ 9^cb] 
9p 



+S.{.CHW + tan 



9 ruH 
(47r)2 167rt;2 



2, + 2M|)(l-4^j (6.7) 



in agreement with Ref. |186| ]. In the situation where we have the hierachy of couplings 
g < gp < An, the second term in Eqs. ( ^^ ) and ( |6.7| ) is suppressed parametrically with 
respect to the first one and could thus be neglected, leading to just a rescaling of the SM 
decay widths [171]. Let us note however that in the case where the ratio g'^/{4:Tr)^ is not 
too small, this additional term could have a numerical impact on the decay rates. 



7. Outlook 



We have described a new framework where BSM physics scenarios can be developed, studied 
and automatically implemented in Monte Carlo or symbolic calculation tools for theoretical 
and experimental investigations. The main purpose of this work has been to contribute 
to streamlining the communication (in both directions) between the theoretical and the 
experimental HEP communities. 

The cornerstone of our approach is the Mathematica package FeynRules where any 
perturbative quantum field theory Lagrangian, renormalizable or not, can be written in a 
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straightforward way and the corresponding Feynman rules obtained automatically. All the 
relevant information can then be passed through dedicated interfaces to matrix element 
generators for Feynman diagram calculations at the tree level or at higher orders. The 
scheme itself looks very simple and is in fact not a new idea. The novelty, however, lies in 
several technical and design aspects which, we believe, constitute a significant improvement 
over the past. 

First, the use of Mathematica as a working environment for the model develop- 
ment gives all the flexibility that is needed for symbolic manipulations. The many built-in 
features of Mathematica, such as matrix diagonalization and pattern recognition func- 
tions, play an important role in building not only robust interfaces for very different codes 
but also to open up new possibilities. For instance, besides implementing BSM models, 
new high-level functionalities/applications can be easily developed by the users themselves, 
made public and possibly included in subsequent official releases. In other words, the code 
is naturally very open to community contributions. A typical application that could take 
advantage of this open structure is the (semi-) automatic development of model calculators 
inside the FeynRules package itself, including mass spectrum and decay width calcula- 
tions. 

Second, the interfaces to MC codes, all of them quite different both in philosophy, 
architecture and aim, offer the possibility of testing and validating model implementations 
at an unprecedented level of accuracy. It also maximizes the probability that a given 
model might be dealt with by at least one matrix element generator. For example, purely 
symbolic generators such as FeynArts/FormCalc can be used for tree- level (or even 
loop) calculations which can then be compared or extended to numerical results from 
CalcHep, MadGraph/MadEvent or Sherpa. In this respect we note that one of 
the current major and common limitations of the matrix-element generators (but not of 
FeynRules) is connected to use of a fixed library of Lorentz and/or gauge structures for 
the vertices. These libraries are in principle extendable, but at present this is done by hand 
and it entails a tedious and often quite long work. The automatization of this part (and 
possibly the reduction of higher-dimensional operators to rcnormalizable ones) through 
FeynRules would be the final step towards full automatization of any Lagrangian into 
Monte Carlo codes. Work in this direction is already in progress. 

With such a framework in place, we hope that several of the current problems and 
drawbacks in new physics simulations faced by the experimental groups will be alleviated 
if not completely solved. First the need of dedicated codes for specific models, which 
then call for long and tedious validations both from the physics point of view as well 
as in their interplay with collaboration softwares, will be greatly reduced. General pur- 
pose tools, from Herwig and Pythia to Sherpa, MadGraph/MadEvent, Whizard, 
CalcHep/CompHep (and potentially also Helac and Alpgen), several of which have 
been successfully embedded in frameworks such as Athena (ATLAS) and CMSSW (CMS), 
offer several ready-to-go solutions for any FeynRules based models. Reducing the pro- 
liferation of highly dedicated tools will greatly simplify the maintenance and reliability of 
the software and more importantly the reproducibility of the MC samples in the mid and 
long terms. In addition, we believe an effort towards making new Lagrangians in Feyn- 
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Rules and the corresponding benchmark points pubhcly available (for example by the 
proponents of the model themselves at the same time of the release of the corresponding 
publication) would certainly be a great advantage for the whole HEP community. 

It is our hope that FeynRules will effectively facilitate interactions between theorists 
and experimentalists. Until now there has not been a preferred way to link the two commu- 
nities. Various solutions have been proposed and used, most of them plagued by significant 
limitations. The best available proposal so far has been to use parton-level events in the 
Les Houches Event File (LHEF) format [187]. This is a natural place to cut a line given 
that theorists and phenomenologists can generate events through various private or public 
tools, and then pass them for showering and hadronization to codes such as Pythia or 
Herwig. These codes are not only already embedded in the experimental software (for the 
following detector simulations), but also have been (or will be) tuned carefully to reproduce 
control data samples. 

While we think this is still a useful approach that should be certainly left open and 
supported, we are convinced the framework we propose is a promising extension. The deep 
reason is that, in our approach, there is no definite line between theory and experiment. 
On the contrary it creates a very extended region where the two overlap and work can be 
done in the same common framework. This leaves much more freedom in where exactly 
the two ends meet and what kind of checks and information can be exchanged. As a result 
there are several practical advantages that come for free. 

As an example, we remind that the LHEF format only standardizes the information 
on the events themselves and some very basic global properties, but any information on 
the physical model (i.e., the explicit form of the Lagrangian) or the parameter choices is 
in general absent. This is of course due to implementation differences among various codes 
which severely compromise any standardization attempt at this level. It is clear that, in 
the long run, this might lead to serious problems of traceability of the MC samples and 
various ambiguities in understanding experimental analyses (such as placing of exclusion 
limits). This problem is of course completely overcome by the approach advocated here, 
since models are now fully and univocally defined. In this sense FeynRules itself offers 
the sought for standardization. 

Another, and maybe even more striking example is that, within FeynRules, model 
building and/or refinement can in principle be done in realtime together with the related 
experimental analyses. One could imagine, for example, that if the TeV world is as rich 
as we hope and as data start showing hints for new particles or effects, a large number of 
competing Lagrangians (and not only benchmark points as used in the typical top-down 
SUSY analyses) could be easily and quickly implemented and readily confronted to data. 
This could be done in a virtuous loop, where theorists and experimentalists "meet" at a 
convenient point of the simulation chain. In other words, various top-down and bottom- 
up studies can fit naturally in this framework, partially addressing often reported worries 
about the actual possibilities to extract precise information on BSM physics from LHC 
data. 

Finally, let us also mention a more long-term advantage of the proposed framework. 
Automatic NLO calculations for SM processes are now clearly in sight, and, in this con- 



-66- 



text, it might be reasonable to ask wlictlicr those developments can be extended to BSM 
processes. We believe the answer is yes, and, since this generalization will probably rely 
on the simultaneous implementation of the model characteristic in different codes {e.g., 
dealing with different parts of the calculation like real and virtual corrections, or analytic 
and numerical results) our approach might also naturally play a crucial role in this context. 
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A. The FeynRules convention for Standard Model inputs 

There are several parameters and particles that have special significance in Fcynman di- 
agram calculators. Examples of these are the strong and electromagnetic couplings, the 
names of the fundamental representation and the structure constants of the strong gauge 
group and so on. For this reason, we have chosen to fix the names of these objects at the 
FeynRules level. Adherence to these standards will increase the chances of successful 
translation. 

The strong gauge group has special significance in many Feynman diagram calculators. 
A user who implements the strong gauge group should adhere to the following rules. The 
indices for the fundamental and adjoint representations of this gauge group should be called 
Colour and Gluon respectively. Furthermore, the names of the QCD gauge boson, cou- 
pling constant, structure constant, totally symmetric term and fundamental representation 
should be given by G, gs, f , dSU3 and T as in the following example: 

SU3C == { 

Abelian -> False, 

GaugeBoson -> G, 

StructureConstant -> f , 
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SymmetricTensor -> dSU3, 
Representations -> {T, Colour}, 
CouplingConstant -> gs 
} 

In addition, the strong coupling constant and its square over Att should be declared in the 
parameter section in the following form: 

\ [Alpha] S == { 



Par amet erType 


-> 


External, 


Value 


-> 


0.118, 


ParameterName 


-> 


aS, 


BlockName 


-> 


SMINPUTS , 


Interact ionOrder 


-> 


{QCD, 2>, 


Description 


-> 


"Strong c 



gs == { 

ParameterType 
Value 

ParameterName 
Interact ionOrder 

} 



-> Internal , 

-> Sqrt [4 Pi \ [Alpha] S] , 
-> G, 

-> {QCD, 1} 



Note that as is given as the external parameter and gs as the internal parameter. The 
description of as may be edited, but it should be remembered that, for the Monte Carlo 
programs that run the strong coupling constant, the value of as should be set at the Z 
pole. For calculation programs that do not run the strong coupling, on the other hand, it 
should be set according to the scale of the interaction. A description may also be added to 
the parameter gs- 

The electromagnetic interaction also has special significance in many Feynman diagram 
calculators and we outline the following standard definitions. The electric coupling constant 
should be called ee, the electric charge should be called Q. The declaration of the electric 
charge should follow the following conventions for naming: 



\ [Alpha] EWMl == { 
ParameterType 
Value 

ParameterName 
BlockName 



-> External, 
-> 127.9, 
-> aEWMl, 
-> SMINPUTS, 



InteractionOrder -> {QED, -2}, 



Description 
}, 

\ [Alpha] EW == { 

ParameterType 
Value 



-> "alpha_EM inverse at the Z pole. 



-> Internal, 

-> l/\ [Alpha] EWMl , 
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InteractionDrder -> {QED, 2}, 
ParameterName -> elEW, 
}, 

ee == { 

ParameterType 
Value 



-> Internal, 

-> Sqrt [4 Pi \ [Alpha] EW ] , 



InteractionDrder -> {QED, 1> 
} 



As for the strong coupling, the description of may be edited , but it should be 

remembered that for calculation programs that run the electric coupling, it should be set 
at the Z pole. For programs which do not run it, the electric coupling should be set at the 
interaction scale. Again, a description may be added to the definition of \ [Alpha] EW and 
ee. 

The Fermi constant and the Z pole mass are very precisely known and are often used in 
calculators to define coupling constants and the scale where couplings are run from. They 
should be included in the SMINPUTS block of the Les Houches accord and should be defined 
by at least the following: 



Gf == { 

ParameterType -> External, 

Value -> 1.16639 * 10" (-5), 

BlockName -> SMINPUTS, 

InteractionDrder -> {QED, 2}, 

Description -> "Fermi constant" 

}, 

ZM == { 

ParameterType -> External, 

Value -> 91.188, 

BlockName -> SMINPUTS, 

Description -> "Z pole mass" 
} 



Moreover, the weak coupling constant name gw and the hypercharge symbol Y are 
used by some calculators and the user is encouraged to use these names where appropri- 
ate. The masses and widths of particles should be assigned whenever possible. If left 
out, FeynRules will assign the value 1 to each. Finally, particles are also identified 
by a PDG number. The user is strongly encouraged to use existing PDG codes in their 
model wherever possible. If not included, a PDG code will be automatically assigned by 
FeynRules beginning at 6000001. 

^^The reason for choosing as the external input parameter, and not asw itself is only to be compliant 
with the Les Houches Accord. 
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Parameter Symbol Value 



Inverse of the electromagnetic coupling a^y^{Mz) 127.9 

Strong coupling as{Mz) 0.118 



Fermi constant 


Gp 


1.16639e-5 GeV 


Z pole mass 


Mz 


91.188 GeV 


c quark mass 


Vflc 


1.42 GeV 


b quark mass 


nib 


4.7 GeV 


t quark mass 


mt 


174.3 GeV 


r lepton mass 


rrir 


1.777 GeV 


Higgs mass 


niH 


120 GeV 


Cabibbo angle 


9c 


0.227736 



Table 7: Input parameters for the SM. 



B. Validation tables 

In this Appendix we report the main results of our work, i.e., the validation tables. In 
general, the tables list quantities (such as decay widths or cross sections) that have no 
direct phenomenological interest but they are physical, easily reproducible and provide 
an exhaustive check of the (complex) values of all the couplings of the model. In several 
instances, other powerful checks (such as gauge invariance, unitarity cancellation at high 
energy, and so on) have been performed that are not presented here. When possible, 
comparisons between the so-called "stock implementation", i.e. implementations already 
available in the Monte Carlo tools, have been made as well as comparisons between different 
Monte Carlo's also in different gauges. All numbers quoted in this section are expressed in 
pico-barn and correspond to a collision in the center-of-mass frame. In all cases agreement 
to better than 1% was obtained. 

B.l The Standard Model 

In this section we give the results for the 35 cross sections tested for the SM between 
CalcHep, MadGraph/MadEvent and Sherpa for a total center of mass energy of 550 
GeV. A pt cut of 20 GeV was applied to each final state particle. The set of external 
parameters used for the test is given in Table 0. A selection of processes is shown in 
Tables | and |. 
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Lepton and weak boson processes in the Standard Model 



Process 


MG-FR MG-ST 


CH-FR CH-ST 


SH-FR SH-ST 


e"*" e~ — >■ e"*" e~ 
e"*" e" — /x"^ /x~ 
e~ — fe I^e 


7.341e+2 7.343e+2 
3.721e-l 3.720e-l 
4.914e+l 4.913e+l 


7.342e+2 7.342e+2 
3.719e-l 3.719e-l 
4.915e+l 4.915e+l 


7.343e+2 7.343e+2 
3.720e-l 3.720e-l 
4.915e+l 4.915e+l 


r+ r- ^ 
r+r- ^ 
r+r- ^ ZA 
t'^t" ^ AA 


5.370e+0 5.360e+0 
3.186e-l 3.180e-l 
2.005e+0 2.007e+0 
2.782e+0 2.780e+0 


5.368e+0 5.368e+0 
3.182e-l 3.182e-l 
2.006e+0 2.006e+0 
2.779e+0 2.779e+0 


5.368e+0 5.368e+0 
3.183e-l 3.183e-l 
2.006e+0 2.006e+0 
2.779e+0 2.779e+0 


ZZ ^ ZZ 

w+w- ^ zz 

HH ZZ 
HH W+W- 


1.960e+0 1.959e+0 
2.726e+2 2.729e+2 
6.268e+l 6.266C+1 
9.449e+l 9.450e+l 


1.961e+0 1.961e+0 
2.726e+2 2.726e+2 
6.266C+1 6.266e+l 
9.447e+l 9.447e+l 


1.961C+0 1.961C+0 
2.726C+2 2.726e+2 
6.266C+1 6.266e+l 
9.448e+l 9.448e+l 



Table 8: Cross sections for a selection of SM production processes. The built-in SM implementation 
in MadGraph and CalcHep and Sherpa are denoted MG-ST, CH-ST and SH-ST, respectively, 
while the FEYNRuLES-generated ones MG-FR, CH-FR, SH-FR. The center-of-mass energy is fixed 
to 550 GeV, and a pr cut of 20 GeV is applied to each final state particle. 
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Quark and gluon processes in the Standard Model 



Process 


MG-FR 


MG-ST 


CH-FR 


CH-ST 


SH-FR 


SH-ST 


GG^ GG 


1.177e+5 


1.178e+5 


1.177e+5 


1.177e+5 


1.177e+5 


1.177e+5 


uu GG 


2.021e+2 


2.021e+2 


2.021e+2 


2.021e+2 


2.021e+2 


2.021e+2 


uu W+ W~ 


1.772e+0 


1.774e+0 


1.774e+0 


1.774e+0 


1.774e+0 


1.774e+0 


uu — ^ Z Z 


1.936e-l 


1.933e-l 


1.935e-l 


1.935e-l 


1.935e-l 


1.935e-l 


uu ^ Z A 


3.380e-l 


3.382e-l 


3.381e-l 


3.381e-l 


3.381e-l 


3.381e-l 


uu^ AA 


1.833e-l 


1.833e-l 


1.832e-l 


1.832e-l 


1.832e-l 


1.832e-l 


uu ^ ss 


9.864e+0 


9.861e+0 


9.868e+0 


9.868e+0 


9.869e+0 


9.869e+0 


ud ^ cs 


3.531e-l 


3.531e-l 


3.531e-l 


3.531e-l 


3.531e-l 


3.532e-l 


us ^ cd 


1.019e-3 


1.019e-3 


1.019e-3 


1.019e-3 


1.019e-3 


1.019e-3 


tt^GG 


6.522e+l 


6.527C+1 


6.528e+l 


6.528e+l 


6.528e+l 


6.528e+l 


ti^ ZA 


i.aiic+o 


1.311C+0 


1.312e+0 


1.312e+0 


1.312e+0 


1.312e+0 


tt^ AA 


8.844C-2 


8.846C-2 


8.849C-2 


8.849e-2 


8.848e-2 


8.848e-2 


tt ^ uu 


1.621e+l 


1.618e+l 


1.619e+l 


1.619e+l 


1.619e+l 


1.619e+l 


ti^W+W- 


1.713e+l 


1.713e+l 


1.713e+l 


1.713e+l 


1.714e+l 


1.714e+l 


tt zz 


1.253e+0 


1.254e+0 


1.253e+0 


1.253e+0 


1.253e+0 


1.253e+0 



Table 9: Cross sections for a selection of SM production processes. The built-in SM implementation 
in MadGraph, CalcHep and Sherpa are denoted MG-ST, CH-ST and Sherpa, respectively, 
while the FEYNRuLES-generated ones MG-FR, CH-FR and SH-FR. The center-of-mass energy is 
fixed to 550 GeV, and a pt cut of 20 GeV is applied to each final state particle. 
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B.2 The general Two-Higgs-doublet model 

In this section we give a selection of the resuhs for the 185 cross sections tested for the 
2HDM between CalcHep, MadGraph/MadEvent and Sherpa for a total center-of- 
mass energy of 800 GeV. A px cut of 20 GeV was applied to each final state particle. The 
set of external parameters used for the test is given in Table 10. A selection of processes 



is shown in Tables 11, 12 and 13 
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Parameter Symbol Value 



Inverse of the electromagnetic coupling 




127.934 


Strong coupling 


as{Mz) 


0.1172 


Fermi constant 


Gf 


1.16637e-5 GeV-^ 


Z pole mass 


Mz 


91.18876 GeV 


c quark pole mass 


rric 


1.25 GeV 


h quark pole mass 


nib 


4.2 GeV 


t quark pole mass 


mt 


174.3 GeV 


T lepton mass 


rrir 


1.777 GeV 


Cosine of the Cabibbo angle 


cos 9c 


0.974589144 


c quark "Yukawa" mass 




0.6 GeV 


h quark "Yukawa" mass 


^Yuk 


3.0 GcV 


t quark "Yukawa" mass 


^Yuk 


175 GeV 


Potential parameters 


Al,2,3 


1.0 




A4 


0.5 




As 


0.4 




Re(A6) 


0.3 




Re(A7) 


0.2 


Charged Higgs mass 


mH± 


300 GeV 


Yukawa parameters (real) 


(rd)2,2 


0.4 GeV 




(rd)2,3 


0.2 GeV 




(rd)3,3 


5 GeV 




(rn)2,2 


2 GeV 




(rM)2,3 


1 GeV 




(rM)3,3 


100 GeV 




(r/)2,2 


0.1 GeV 




(r/)2,3 


0.5 GeV 




{^1)3,3 


3 GeV 



Table 10: Input parameters for the 2HDM. All parameters that are not quoted have a zero value. 
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Lepton processes in the Two-Higgs-Doublet Model 



Procoss 


:mg-fr 


^IG-ST 


CII-FE 


SII-FR 


r+ 

1 


/ 




1.917e-5 


1.916e-5 


1.916e-5 


1.916e-5 


r+ 




— y hi hn 
^ /(-l iti2 


2.043e-3 


2.043e-3 


2.043e-3 


2.043e-3 


r+ 


7" 




2.043e-3 


2.041e-3 


2.042e-3 


2.043e-3 


T+ 




— >■ h^ hn 


2.349e-5 


2.349e-5 


2.349e-5 


2.348e-5 


r+ 




— > hn ho 


2.260e-4 


2.259C-4 


2.259C-4 


2.259e-4 


r+ 




— > h'i hri 


7.242e-4 


7.240e-4 


7.241e-4 


7.240e-4 


r+ 




/i+ h~ 


1.345e-2 


1.345e-2 


1.345e-2 


1.345e-2 


r+ 


1 


^ Z hi 


5.818e-4 


5.824C-4 


5.824C-4 


5.821e-4 


r+ 

1 


/ 


Z hn 


8.012e-3 


8.016e-3 


8.012e-3 


8.011e-3 


r+ 


7" 


Z h-i 


4.844e-3 


4.842e-3 


4.842e-3 


4.842e-3 


r+ 


'y- 


— y hi 


2.915e-3 


2.915e-3 


2.914e-3 


2.913e-3 


r+ 


T~ 


— > "y h.o 


8.059e-5 


8.059e-5 


8.060e-5 


8.058e-5 


r+ 


7-" 


— > T h.Q 


4.347e-3 


4.348e-3 


4.347e-3 


4.346e-3 


r+ 




w~ ^+ 


2.038e-3 


2.037e-3 


2.037e-3 


2.037e-3 


r+ 


T~ 


hr 


2.037e-3 


2.039e-3 


2.037e-3 


2.037e-3 


T+ 


T~ 


z z 


1.782e-l 


1.784e-l 


1.783e-l 


1.784e-l 


r+ 


T~ 




3.017C+0 


3.018e+0 


3.018e+0 


3.018e+0 


r+ 


7- 


— > 11^ n 


1.756C-1 


1.753e-l 


1.755e-l 


1.755e-l 


r+ 


7- 




1.453e-8 


1.453e-8 


1.453e-8 


1.452e-8 


r+ 


-7-" 


— > It 


1.452e-8 


1.453C--8 


1.453e-8 


1.452e-8 


r+ 


7- 


— ^ T 


7.423e+2 


7.421e+2 


7.421e+2 


7.422e+2 


r+ 


7" 


— ^ ct 


1.189e-7 


1.189e-7 


1.189e-7 


1.189e-7 


r+ 


r~ 


tc 


1.189e-7 


1.189e-7 


1.189e-7 


1.189e-7 


r+ 




tt 


2.690e-l 


2.686e-l 


2.687e-l 


2.688e-l 


r+ 


T~ 


ss 


1.431e-l 


1.432e-l 


1.431e-l 


1.433e-l 


r+ 


^_ 


sh 


5.252e+0 


5.252e-9 


5.249e-9 


5.242e-9 


r+ 


r~ 


hs 


5.248e+0 


5.249e-9 


5.249e-9 


5.242e-9 


r+ 




hi 


1.431e-l 


1.431e-l 


1.431e-l 


1.431e-l 



Table 11: Cross sections for a selection of r+r" initiated processes in the 2HDM. The built-in 
2HDM implementation in MadGraph is denoted MG-ST, while the FEYNRuLES-generated ones 
are MG-FR, CH-FR and SH-FR. The center-of-mass energy is fixed to 800 GeV, and a pt cut of 
20 GeV is applied to each final state particle. 
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Neutrino processes in the Two-Higgs-Doublet Model 



Process 


MG-FR 


MG-ST 


CH-FR 


SH-FR 


r+ 


u 






1.440e-2 


1.439e-2 


1.440e-2 


1.441e-2 


r+ 






/i+ ho 


4.277e-3 


4.277e-3 


4.278e-3 


4.277e-3 


r+ 


u 




h+ /i3 


4.063e-3 


4.066e-3 


4.064e-3 


4.065e-3 


r+ 






h+ Z 


1.731e-3 


1.731e-3 


1.731e-3 


1.731e-3 


r+ 








1.271e-3 


1.271e-3 


1.271e-3 


1.271C-3 


r+ 






hi 


9.814e-4 


9.803e-4 


9.808e-4 


9.807e-4 


r+ 








1.802e-2 


1.801e-2 


1.802e-2 


1.802e-2 


r+ 








1.076e-2 


1.076e-2 


1.077e-2 


1.077e-2 


r+ 






Z 


2.250e+0 


2.251e+0 


2.251e+0 


2.250e+0 


r+ 






vr 1 


1.439e+0 


1.439e+0 


1.439e+0 


1.440e+0 


r+ 






p'^ V^ 


1.816e-l 


1.816e-l 


1.816e-l 


1.816e-l 


r+ 








1.816e-l 


1.816e-l 


1.816e-l 


1.816e-l 


r+ 


Vr 




^+ Vr 


1.002e-8 


1.002e-8 


1.002e-8 


1.002e-8 


r+ 


Vr 




Vr 


9.179e+0 


9.180e+0 


9.180e+0 


9.180C+0 


r+ 


Ur 




u d 


5.174e-l 


5.175e-l 


5.174e-l 


5.175C-1 


r+ 


Vr 




u s 


2.733e-2 


2.733e-2 


2.733e-2 


2.734e-2 


r+ 


Vr 




cd 


2.733e-2 


2.734e-2 


2.733e-2 


2.734e-2 


r+ 


Vr 




cs 


5.175e-l 


5.174e-l 


5.174e-l 


5.184e-l 


r+ 


Vr 




cb 


1.202e-7 


1.202e-7 


1.202e-7 


1.202e-7 


r+ 


Vr 




ts 


4.363e-9 


4.364e-9 


4.364e-9 


4.364e-9 


r+ 


Vr 




tb 


5.097e-l 


5.097e-l 


5.098e-l 


5.099e-l 



Table 12: Cross sections for a selection of T~^Vr initiated processes in the 2HDM. The built-in 
2HDM implementation in MadGraph is denoted MG-ST, while the FEYNRuLES-generated ones 
are MG-FR, CH-FR and SH-FR. The center-of-mass energy is fixed to 800 GeV, and a pr cut of 
20 GeV is applied to each final state particle. 



-76- 



Weak boson processes in the Two-Higgs-Doublet Model 



Procoss 


MG-FE 


MG-ST 


CII-FR 


SII-FR 


W 




1.347e+3 


1.346e+3 


1.347e+3 


1.347e+3 


w 


~ ^ z z 


2.787e+2 


2.774e+2 


2.782e+2 


2.782e+2 


w 


' ^ z^ 


1.510e+2 


1.513e+2 


1.511e+2 


1.509e+2 


z z ^ 


z z 


1.616e+l 


1.615e+l 


1.616e+l 


1.616e+l 


w 


- ^ 1^+ /j- 


3.589e-l 


3.586e-l 


3.588e-l 


3.587e-l 


W 


" W~ h+ 


3.586e-l 


3.589e-l 


3.588e-l 


3.587C-1 


w 


- ^ Zhi 


1.788C-1 


1.787e-l 


1.787e-l 


1.788C-1 


w 




3.680C+1 


3.685C+1 


3.685C+1 


3.684C+1 


W+ W 

vv vv 




2.130e+l 


2.131e+l 


2.130e+l 


2.129e+l 


z z ^ 


Zhi 


5.934e-l 


5.943e-l 


5.940e-l 


5.939e-l 


w 




6.617e-l 


6.620e-l 


6.618e-l 


6.617e-l 


W+ W 

vv vv 


~" ^ hrt fin 


4.455e-l 


4.449e-l 


4.453e-l 


4.452e-l 


W+ W 

vv vv 




1.077e+0 


1.076e+0 


1.076e+0 


1.076e+0 


w+ w 

vv vv 


~ ^ h-^ h~ 

~ 1 V It 


2.446e+0 


2.444e+0 


2.443e+0 


2.443e+0 


Z - 


/i+ hi 


2.888e-l 


2.888e-l 


2.887e-l 


2.887e-l 


z - 


/;+ ho 

^ 1 V 11/2 


1.485e-2 


1.486e-2 


1.485e-2 


1.485e-2 


Z - 


/?+ h-i 


3.014e-2 


3.013C-2 


3.012C-2 


3.011e-2 


M/+ 'V - 

vV f 


/?+ /)! 

' 1 i / fc 


1.970e-2 


1.969e-2 


1.969e-2 


1.969e-2 


7 - 


/l+ /l9 

' 1 fj 1 ij £ 


7.925e-3 


7.927e-3 


7.926C-3 


7.925e-3 


7 - 




1.692e-2 


1.693e-2 


1.693e-2 


1.693e-2 


z z ~^ 


hi hi 


2.333e+0 


2.333e+0 


2.335e+0 


2.334e+0 


z z ^ 


ho ho 


6.636e-l 


6.627e-l 


6.630e-l 


6.629e-l 


zz ^ 


h2 /l3 


1.065e+0 


1.066e+0 


1.066e+0 


1.065e+0 


zz ^ 


h+ h- 


1.356e+0 


1.356e+0 


1.356e+0 


1.356e+0 


Zhi 


h2 hz 


2.929e+0 


2.928e+0 


2.930e+0 


2.929e+0 


W+ hi 


^h+ h2 


4.732e+0 


4.727e+0 


4.731e+0 


4.729e+0 


W+ h2 


^h+ h2 


2.998e+0 


2.997e+0 


2.999e+0 


2.997e+0 


W+ h- 


h2 hz 


2.294e+0 


2.292e+0 


2.293e+0 


2.292e+0 



Table 13: Cross sections for a selection of processes in the 2HDM with two weak bosons in the 

initial state. The built-in 2HDM implementation in MadGraph is denoted MG-ST, while the 
FEYNRuLES-generated ones are MG-FR, CH-FR and SH-FR. The center-of-mass energy is fixed to 
800 GeV, and a pt cut of 20 GeV is applied to each final state particle. 
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B.3 The most general Minimal Supersymmetric Standard Model 

In order to fully determine the MSSM Lagrangian at low energy scale, it is sufficient to 
fix the SM sector and the super symmetry-breaking scenario. We have chosen the typi- 
cal minimal supergravity point SPS la |182] which is completely defined once we fix the 
values of four parameters at the gauge coupling unification scale and one sign. The com- 
plete set of input parameters is shown in Table 14. We then use the numerical program 
SOFTSUSY |18J] to solve the renormalization group equations linking this restricted set 
of supersymmetry-breaking parameters at high-energy scale to the complete set of masses, 
mixing matrices and parameters appearing in £mssm at the weak scale. The output is 
encoded in a file following the SLHA conventions, readable by FeynRules after the use 
of an additional translation interface taking into account the small differences between 
the SLHA2 format and the one of our implementation. This interface is available on the 
FeynRules website, as well as an the corresponding SLHA2 output file. 

In Tables O, 16, 17 and 18, we give some examples of the numerical checks which we 



have performed in order to validate the implementation of the MSSM in FeynRules. We 
recah that the built-in MSSM implementation in MadGraph and CalcHep are denoted 
MG-ST and CH-ST, respectively, while the FEYNRuLES-generated ones MG-FR and CH- 
FR. A pt cut of 20 GeV is applied to all final-state charged leptons, photons, and jets 
(including 6-jets). The full list of results can be found on the FeynRules webpage. 



Parameter 


Symbol 


Value 


Inverse of the electromagnetic coupling 


aEwiMz) 


127.934 


Strong coupling 


as{Mz) 


0.1172 


Fermi constant 


Gp 


1.16637e-5 GeV'^ 


Z mass 


Mz 


91.1876 GeV 


b quark mass 


mb{mb) 


4.25 GeV 


t quark mass 


nit 


175 GeV 


r lepton mass 


m-r 


1.777 GeV 


Universal scalar mass 


mo 


100 GeV 


Universal gaugino mass 




250 GeV 


Universal trilinear coupling 


Ao 


-100 GeV 


Ratio of the two vevs 


tan/3 


10 


Off-diagonal Higgs mixing parameter 


sign(^) 


+ 



Table 14: Input parameters for the SPS la benchmark point for the MSSM. 
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Decay widths in the MSSM 



Process 


MG-FR MG-ST 


Process 


MG-FR MG-ST 


t ^ W+b 

W+ cs 
bb 


1.5608 1.5561 
6.6807e-l 6.6830e-l 
4.5479e-l 4.5479e-l 


Z ^ dd 
bb 

H+ tb 


3.6625e-l 3.6712e-l 
1.6268e-3 1.6268e-3 
4.3833e-l 4.3833e-l 


ui xt^ 
di Xid 


1.3661 1.3661 
2.8196e-l 2.8196e-l 
1.4545C-1 1.4545e-l 


U2 ^ XlU 

d5 ^ XiU 
h Xi^" 


1.1373 1.1373 
3.2220 3.2220 
2.1612e-l 2.1612e-l 


~g ^ ell b 

xl ^ xt w- 

Xt ^t^T 


5.5408e-l 5.5370e-l 
5.6624e-l 5.6460e-l 
1.5768e-2 1.5748e-2 


xl - li r+ 

xl ^ xt w- 
xt ^ xlw+ 


9.1581e-3 9.1507e-3 
6.4519e-l 6.4576e-l 
7.2945e-l 7.2998e-l 



Table 15: Widths (in GeV) of some of the allowed decay channels in the SPS la scenario. MG-FR 
amd MG-ST denote the FEYNRuLES-generated and built-in MadGraph implementations. 



- 79 - 



Higgs production in the MSSM 



Process 


MG-FR MG-ST 


CH-FR CH-ST 


e+ e- ^ ZhP 
e+ e- H+ H- 

T+ T- ^hPH"^ 

r+ r- ^ ^0 

T- Ur H- hP 

T~ W- 

T~ Vt ^ Z H~ 


8.787e-3 8.788e-3 
8.121e-3 8.121e-3 
1.610e-5 1.610e-5 
1.741e-5 1.741e-5 
6.245e-6 6.245e-6 
1.810e-2 1.811e-2 
3.125e-2 3.123e-2 


8.787e-3 8.787e-3 
8.119e-3 8.119e-3 
1.610e-5 1.610e-5 
1.741e-5 1.741e-5 
6.244e-6 6.243e-6 
1.810e-2 1.810e-2 
3.124e-2 3.124e-2 


uu ^ Z hP 
dd ZH^ 
bb^A^ A^ 
bb H+ H- 
bt^H- H° 
bt ZH- 


3.331e-3 3.331e-3 
4.213e-7 4.215C-7 
7.221e-5 7.218C-5 
9.240e-4 9.237e-4 
2.070e-3 2.070e-3 
2.590e-2 2.587e-2 


3.331e-3 3.331C-3 
4.214C-7 4.214C-7 
7.220C-5 7.214C-5 
9.237e-4 9.237e-4 
2.069e-3 2.069e-3 
2.592e-2 2.592e-2 


W+ W- hP 

w+ w- zhP 

W+ W- H+ H- 

ZZ^h"" 
Z-i H+ H- 
W Z ^W- A° 
W-j^W- hP 


1.109e-3 l.llOe-3 
8.217e+l 8.216e+l 
3.689e-2 3.686e-2 
7.829e+0 7.827e+0 

1.126e-2 1.124e-2 
4.375e-4 4.387e-4 
1.588e+l 1.589e+l 


l.llOe-3 l.llOe-3 
8.213e+l 8.213e+l 
3.685e-2 3.685e-2 
7.827e+0 7.827e+0 

1.124C-2 1.124e-2 
4.378e-4 4.378e-4 
1.589e+l 1.589e+l 



Table 16: Cross sections for a selection of Higgs production processes in the MSSM scenario SPS 
la. The built-in MSSM implementation in MadGraph and CalcHep are denoted MG-ST and 
CH-ST, respectively, while the FEYNRuLES-generated ones are MG-FR and CH-FR. The center- 
of-mass energy is fixed to 1200 GeV. 



-80- 



Supersymmetric particle production from fermions 



Process 


MG-FR 


MG-ST 


CH-FR 


CH-ST 


6+ 6~ 




1.944e- 


-1 


1.944e- 


-1 


l,944e- 


-1 


1.944c- 


-1 


6+ &~ 




4.862e- 


-1 


4.863e- 


-1 


4.862e- 


-1 


4.863c- 


1 


6+ e~ 


— > Ua U*4 


1.664e- 


-3 


1.662e- 


-3 


1.663e- 


-3 


1.663e- 


-3 


6+ 6~ 


— > dl 


2.858e- 


-4 


2.857e- 


-4 


2.858e- 


-4 


2.858e- 


-4 


T+ T~ 


^ ^1 ^6 


4.332e- 


-2 


4.326e- 


-2 


4.329e- 


-2 


4.382e- 


-2 


T~ Vj- 




1.206e- 


-1 


1.206e- 


-1 


1.206e- 


-1 


1.206e- 


-1 


6+ 6 


^ yS yS 


4.334e- 


-2 


4.334e- 


-2 


4.330e- 


-2 


4.330e- 


-2 


6+ 6 


— V vZ^ V-. 
' Al Al 


1.023e- 


-1 


1.023e- 


-1 


1.023e- 


-1 


1.023e- 


-1 




^ yS y2 

A3 A3 


9 280e- 


-5 


9.267e- 


-5 


9.272e- 


-5 


9.271e- 


-5 




— ^ Y^ yI7 

Al A2 


1.169e- 


-2 


1.170e- 


-2 


1.170e- 


-2 


1.170e- 


-2 




— > Yi Y? 
Al Al 


5.445e- 


-2 


5.444e- 


-2 


5.446e- 


-2 


5.446e- 


-2 




— > Yo Y? 
A2 A4 


5.306e- 


-2 


5.313e- 


-2 


5.306e- 


-2 


5.306e- 


-2 


uu — 


'3 '3 


2.123e- 


-3 


2.123e- 


-3 


2.123e- 


-3 


2.123e- 


-3 


uu — 


* U,r: ilX 


6.141e- 


-1 


6.142e- 


-1 


6.141e- 


-1 


6.141e- 


-1 


dd—> 




3.607e- 


-3 


3 606e- 


-3 


3.607e- 


-3 


3.Q07e- 


-3 


dd 




1.306e- 


-1 


1.307e- 


-1 


1.307e- 


-1 


1.307e- 


-1 


bb^ 


dxd\ 


3.403e- 


-1 


3.401e- 


-1 


3.401e- 


-1 


3.401e- 


-1 


bb^ 


u% u\ 


5.913e- 


-3 


5.915e- 


-3 


5.912e- 


-3 


5.270e- 


-3 


bt^ 


Iq Dl 


1.114e- 


-2 


1.114e- 


-2 


1.115e- 


-2 


l.llSe- 


-2 


bt^ 


d-i ul 


4.417e- 


-1 


4.419e- 


-1 


4.420e- 


-1 


4.420e- 


-1 


uu 




1.272c- 


-4 


1.272e- 


-4 


1.272c- 


-4 


1.272c- 


-4 


uu — 


' Xt X2 


1.321e- 


-2 


1.317e- 


-2 


1.319e- 


-2 


1.320c- 


-2 


bb^ 


X3 A4 


1.087e- 


-2 


1.087e- 


-2 


1.085e- 


-2 


1.085e- 


-2 


bb^ 


X2 


9.556e- 


-2 


9.545e- 


-2 


9.556c- 


-2 


9.556c- 


-2 


bt^ 


xr X2 


1.556e- 


-2 


1.557e- 


-2 


1.557c- 


-2 


1.557c- 


-2 


bi^ 


X2 X3 


3.981e- 


-2 


3.971e- 


-2 


3.977e- 


-2 


3.977e- 


-2 



Table 17: Cross sections for a selection of supersymmetric particle pair production processes in 
the MSSM scenario SPS la. The built-in MSSM implementation in MadGraph and CalcHep 
are denoted MG-ST and CH-ST, respectively, while the FEYNRULES-generated ones are MG-FR 
and CH-FR. The center-of-mass energy is fixed to 1200 GeV. 
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Supersymmetric particle production from gauge bosons 



Process 


MG-FR 


MG-ST 


CH-FR 


CH-ST 


W+ W- [j" 














1 fi«9p_ 


.9 


W+ W~ i>3 Z>3 




-2 




-2 




-2 




-2 


W+ W- ui ul 




-2 


5 34Qp- 


-2 




-2 


q 1 83p- 


-1 


Z Z ^viul 


7.678e- 


-2 


7.686e- 


-2 


7.686e- 


-2 


7.686e- 


-2 


Z Z ^ d^d^ 


4.693e- 


-2 


4.695e- 


-2 


4.693e- 


-2 


4.693e- 


-2 


Z ^ —>■ u \ ui 


3.283e- 


-2 


3.285e- 


-2 


3.286e- 


-2 


3.286e- 


-2 


Z 'y ^ lo it 


1.712e- 


-2 


1.711e- 


-2 


1.712c- 


-2 


1.712c- 


-2 


W- Z ^ 17 z>? 

O 


3.952c- 


-2 


3.941c- 


-2 


3.950c- 


-2 


3.950c- 


2 


W- Z ^(h ul 


2.690e- 


-2 


2.689e- 


-2 


2.690c- 


-2 


2.690c- 


-2 


J d2ul 


3.618e- 


-4 


3.618e- 


-4 


3.618c- 


-4 


3.618c- 


-4 




1.129e- 


-1 


1.129e- 


-1 


1.129e- 


-1 


1.129e- 


■1 


g Z ^ d4:d\ 


4.637e- 


-3 


4.633e- 


-3 


4.634e- 


-3 


4.634e- 


-3 


gW- ^dr, ul 


2.569e- 


-1 


2.566e- 


-1 


2.566e- 


-1 


2.566e- 


-1 


gW^^dluQ 


2.208e- 


-2 


2.206e- 


-2 


2.206e- 


-2 


2.206e- 


-2 


99 — ll-A "3 


i.SGoc- 


-1 


1.8600- 


-1 


1.866('- 


-1 


I.8G60- 


1 


w+ w-^xUs 


6.514e- 


-1 


6.515e- 


-1 


6.515e- 


-1 


6.515e- 


-1 


w+w-^xtxT 


1.836e+0 


1.837e+0 


1.835e+0 


1.835e+0 


11 ^xt Xi 


6.250e- 


-1 


6.263e- 


-1 


6.257e- 


-1 


6.257e- 


■1 


w-z^xlxi 


4.738e- 


-1 


4.751e- 


-1 


4.746c- 


-1 


4.746e- 


-1 


W-^^X2 X2 


5.235e- 


-2 


5.235e- 


-2 


5.236e- 


-2 


5.236e- 


■2 



Table 18: Cross sections for a selection of supersymmetric particle pair production processes in 
the MSSM scenario SPS la. The built-in MSSM implementation in MadGraph and CalcHep 
are denoted MG-ST and CH-ST, respectively, while the FEYNRuLES-generated ones are MG-FR 
and CH-FR. The center-of-mass energy is fixed to 1200 GeV. 
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Parameter 


Symbol 


Value 


Inverse of the electromagnetic coupling 




127.9 


Strong coupling 




0.1172 


Fermi constant 


Gp 


1.16637e-5 GeV-^ 


Z mass 


Mz 


91.1876 GeV 


W mass 


Mw 


80.398 GeV 


W' mass 


Mw' 


500 GeV 


Heavy fermion mass 


Mp 


4 TeV 



Table 19: Input parameters for the MHM. 



B.4 The Minimal Higgsless Model 

The external parameters used for the validation of the Minimal Higgsless Model are shown 
in Table 19, and all the widths are set to zero. We employ a center-of-mass energy of 600 
GeV and a pp cut of 20 GeV if only SM particles are present, a center-of-mass energy of 
1200 GeV and a px cut of 200 GeV if heavy vector bosons are present but heavy fermions 
not, and a center-of-mass energy of 10000 GeV and a pp cut of 2000 GeV if heavy fermions 
are present. 
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Strong Processes in the Minimal Higgsless Model 



Process 


CH-LH-F 


CH-FR-F 


CH-FR-U 


SH-FR-U 


MG-FR-U 


G,G^ 


G,G 


1.143e+5 


1.143e+5 


1.143e+5 


1.143e+5 


1.143e+5 


u,u ^ 


G,G 


1.705e+2 


1.705e+2 


1.705e+2 


1.705e+2 


1.706e-h2 


U,u 


G,G 


O.OOOe-1 


O.OOOe-1 


O.OOOe-1 


O.OOOe-1 


O.OOOe-1 


U,Lf ^ 


G,G 


8.696e-2 


8.696e-2 


8.696e-2 


8.696C-2 


8.690e-2 


u,u — 


> t,t 


8.116C+0 


8.116C+0 


8.116e+0 


8.117C+0 


8.100C+0 


U,u- 


> t,t 


6.464C-2 


6.464e-2 


6.464C-2 


6.467e-2 


6.466e-2 


U,u- 


T,t 


6.592e-l 


6.592e-l 


6.592e-l 


6.594e-l 


6.597e-l 


U,u^ 


t,f 


6.592e-l 


6.592e-l 


6.592e-l 


6.594e-l 


6.596e-l 


U,U - 


^t,t 


5.033e-2 


5.033e-2 


5.033e-2 


5.035e-2 


5.044e-2 


U,U^T,t 


1.366e-l 


1.366e-l 


1.366e-l 


1.367e-l 


1.365e-l 


u, u - 


T,T 


1.589('-1 


1.589e-i 


i.589c'-l 


1.589e-l 


1.589C-1 



Table 20: Cross sections for a selection of strong processes in the Minimal Higgsless Model. The 
LANHEP-generated MHM implementation in CalcHep is denoted CH-LH, while the FeynRules- 
generated ones are MG-FR, CH-FR and SH-FR. F means the calculation was done in Feynman 
gauge while U means it was done in unitary gauge The center-of-mass energy is fixed to 600 GeV 
and a pT cut of 20 GcV is applied to each final state particle, if only SM particles are present. A 
center-of-mass energy of 1200 GeV and a Pt cut of 200 GeV is used if heavy vector bosons are 
present but heavy fermions not, and a center-of-mass energy of 10000 GeV and a pt cut of 2000 
GeV if heavy fermions are present. 
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Charged Electroweak Processes in the Minimal Higgsless Model 



Process 


/^TT T TT TTl 


CH-i' K-r 




SH-l^K-U 


IV TTi"r> T T 

MCj-i' K- U 


z,w+ w+,z 


3.121e+2 


3.121e+2 


3.121e+2 


3.119e+2 


3.118e+2 


w'+,z ^w+,z 


4.084e+0 


4.084e+0 


4.084e+0 


4.077e+0 


4.080e+0 


z',w+ w+,z 


4.214e+0 


4.214e+0 


4.214e+0 


4.210e+0 


4.214e+0 


z',w'+ w+,z 


3.072e+l 


3.072e+l 


3.072e+l 


3.071e+l 


3.074e+l 


w'+,z ^ w'+,z 


2.289e+l 


2.289C+1 


2.289C+1 


2.289C+1 


2.288C+1 


w'+,z ^ w+,z' 


1.303e+2 


1.303e+2 


1.303e+2 


1.302e+2 


1.304e+2 


z ,w ^ ^ ,z 


8.098e+0 


8.098e+0 


8.098e+0 


8.095e+0 


8.111e+0 


z\w'^ w+,z' 


1.914e+l 


1.914e+l 


1.914e+l 


1.913e+l 


1.915e+l 


z ,w + + ,z 


6.967e+2 


6.967e+2 


6.967e+2 


6.963e+2 


6.965e+2 


u,d^ Z, W+ 


1.112e+0 


1.112e+0 


1.112e+0 


l.llOe+0 


1.112e+0 


u,d^ Z, W'+ 


8.639e-2 


8.639e-2 


8.639e-2 


8.638e-2 


8.647e-2 


u,d^ Z',W+ 


8.463e-2 


8.463e-2 


8.463e-2 


8.460e-2 


8.456e-2 


U,d^ Z,W+ 


7.987e-2 


7.987e-2 


7.987e-2 


7.986e-2 


7.996e-2 


U,D ^ Z, W+ 


1.794e+0 


1.794e+0 


1.794e+0 


1.793e+0 


1.793e+0 


U,d^ Z',W+ 


1.464e-l 


1.464e-l 


1.464e-l 


1.464e-l 


1.463e-l 


U,d^ Z,W'+ 


1.528e-l 


1.528e-l 


1.528e-l 


1.529e-l 


1.529e-l 


U,D ^ Z',W+ 


2.499e+0 


2.499C+0 


2.499e+0 


2.499e+0 


2.499e+0 


U,D ^ Z, W'+ 


2.564e+0 


2.564e+0 


2.564e+0 


2.563e+0 


2.562e+0 


U,d^ Z',W'+ 


2.512C-1 


2.512e-l 


2.512C-1 


2.512e-l 


2.513e-l 


U,D ^ Z',W'+ 


2.287C+0 


2.287e+0 


2.287C+0 


2.286e+0 


2.286e+0 


e, v'l b,t 


9.514C-1 


9.514e-l 


9.514C-1 


9.515e-l 


9.514C-1 


e, ui ^ B,t 


5.346e-5 


5.346e-5 


5.346e-5 


5.347e-5 


5.339e-5 


E-,ui ^ B,i 


8.287e+0 


8.287e+0 


8.287e+0 


8.289e+0 


8.287e+0 


e,Pi ^ B,T 


1.106e-3 


1.106e-3 


1.106e-3 


1.106e-3 


1.107e-3 


E-,Ni ^b,t 


8.071e-2 


8.071e-2 


8.071e-2 


8.075e-2 


8.054e-2 


E-,vi B,f 


5.904e-l 


5.904e-l 


5.904e-l 


5.906e-l 


5.905e-l 


E-,Ni ^ B,t 


1.717e+0 


1.717e+0 


1.717e+0 


1.717e+0 


1.715e+0 


E-,Ni ^ B,f 


1.732e+0 


1.732e+0 


1.732e+0 


1.733e+0 


1.733e+0 



Table 21: Cross sections for a selection of charged electroweak processes in the Minimal Higgsless 
Model. The LANHEP-generated MHM implementation in CalcHep is denoted CH-LH, while the 
FEYNRuLES-generated ones are MG-FR, CH-FR and SH-FR. F means the calculation was done in 
Feynman gauge while U means it was done in unitary gauge The centcr-of-mass energy is fixed to 
600 GeV and a px cut of 20 GeV is applied to each final state particle, if only SM particles are 
present. A center-of-mass energy of 1200 GeV and a pr cut of 200 GeV is used if heavy vector 
bosons are present but heavy fermions not, and a center-of-mass energy of 10000 GeV and a pt cut 
of 2000 GeV if heavy fermions are present. 
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Neutral Electroweak Processes in the Minimal Higgsless Model 







T T 17 O T7 

V iJL-Jb JK-r 


C I L-r K-L 


bil-t Jrv-L 


A ir^ iro TT 

AilT-_b JK- U 


w+,w- - 


► w+,w- 


1.407e+3 


1.407e+3 


1.407e+3 


1.407e+3 


1.407e-F3 


w+,w- 


w+,w'- 


2.917e+0 


2.917e+0 


2.917e+0 


2.917e+0 


2.918e-F0 


w+,w'- - 


> w+,w- 


5.313e+0 


5.313e+0 


5.313e+0 


5.313e+0 


5.316e-F0 


w+,w- 


w'+,w'~ 


7.078C+0 


7.078C+0 


7.078C+0 


7.074C+0 


7.078C+0 


w+,w'- - 


> w+,w'- 


1.406C+2 


1.406C+2 


1.406C+2 


1.405e+2 


1.404c-h2 


w'+,w'- - 


w+,w- 


3.123C+1 


3.123C+1 


3.123C+1 


3.123C+1 


3.126C+1 


w'+,w'- - 


> w+,w'- 


1.948C+1 


1.948C+1 


1.948C+1 


1.947e+l 


1.944C+1 


w+,w'- 


w'+,w'- 


6.991e+0 


6.991C+0 


6.991C+0 


6.995C+0 


6.974C-F0 


w'+,w'- - 


y W'+,W'~ 


7.245e+2 


7.245e+2 


7.245e+2 


7.241e-i-2 


7.256e-F2 


u,u ^ 


z,z 


1.957e-l 


1.957e-l 


1.957e-l 


1.958e-l 


1.955e-l 


u,u ^ 


z,z' 


9.367e-5 


9.367e-5 


9.367e-5 


9.362e-5 


9.368e-5 


U,u^ 


z,z 


4.497e+0 


4.497e+0 


4.497e+0 


4.497e+0 


4.494e-F0 


U,U - 


y z,z 


6.302e-l 


6.302e-l 


6.302e-l 


6.302e-l 


6.304e-l 


U,u^ 


z, z' 


9.733e+0 


9.733e+0 


9.733e+0 


9.732e+0 


9.744e-F0 


u,u ^ 


z\z' 


1.196e-5 


1.196e-5 


1.196e-5 


1.196e-5 


1.196e-5 


U,U 


z,z' 


1.371e+0 


1.371e+0 


1.371e+0 


1.371e+0 


1.371e-F0 




z\z' 


5.268e+0 


5.268C+0 


5.268e+0 


5.268C+0 


5.266c-h0 


u,u ^ 


z\z' 


8.832e-l 


8.832C-1 


8.832e-l 


8. 832c- 1 


8.828e-l 


e,e — 


> t,t 


4.480e-l 


4.480c- 1 


4.480e-l 


4.480c- 1 


4.478e-l 


e,e — 


T,i 


1.388e-5 


1.388e-5 


1.388e-5 


1.388C-5 


1.389e-5 


E-,e- 




1.978e+0 


1.978e+0 


1.978e+0 


1.978e-F0 


1.978e-F0 


e, e — > 


T,T 


4.093e-4 


4.093e-4 


4.093e-4 


4.093e-4 


4.093e-4 


E-,E 


^t,t 


4.155e-2 


4.155e-2 


4.155e-2 


4.157e-2 


4.153e-2 


E-,e^T,T 


1.359e-l 


1.359e-l 


1.359e-l 


1.359e-l 


1.359e-l 


E-,E- 




4.097e-l 


4.097e-l 


4.097e-l 


4.098e-l 


4.095e-l 


E-,E^ T,T 


4.207e-l 


4.207e-l 


4.207e-l 


4.207e-l 


4.202e-l 



Table 22: Cross sections for a selection of neutral electroweak processes in the Minimal Higgsless 
Model. The LANHEP-generated MHM implementation in CalcHep is denoted CH-LH, while the 
FEYNRuLES-generated ones are MG-FR, CH-FR and SH-FR. F means the calculation was done in 
Feynman gauge while U means it was done in unitary gauge The center-of-mass energy is fixed to 
600 GeV and a pr cut of 20 GeV is applied to each final state particle, if only SM particles are 
present. A center-of-mass energy of 1200 GeV and a px cut of 200 GeV is used if heavy vector 
bosons are present but heavy fermions not, and a center-of-mass energy of 10000 GeV and a pt cut 
of 2000 GeV if heavy fermions are present. 
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Parameter 


Symbol 


Value 


Inverse of the electromagnetic coupling 




128 


Strong coupling 




0.1172 


Fermi constant 


Gp 


1.16639e-5 GeV-^ 


Z pole mass 


Mz 


91.1876 GeV 


c quark mass 


rric 


1.42 GeV 


h quark mass 


TJlb 


4.2 GeV 


t quark mass 




175 GeV 


r lepton mass 


rUr 


1.777 GeV 


Higgs mass 


niH 


120 GeV 


Sine of the electroweak mixing angle 


Sw 


0.48076 


Parameter of the CKM matrix 


Sl2 


0.221 


Parameter of the CKM matrix 


S23 


0.041 


Parameter of the CKM matrix 


Sl3 


0.0035 


Extra dimensional radius 


R 


0.002 GeV"^ 


Cutoff scale 


A 


10 TeV 



Table 23: Input parameters for our MUED benchmark scenario. 



B.5 Extra Dimensional models 

The benchmark point which we have used for our numerical comparison of the MUED 
implementation in FeynRules to the existing one in CalcHep is defined through the 
various external parameters shown in Table 23. The masses of the first Kaluza-Klein 
excitations are computed via one-loop calculations |16£ , 189 1. We obtain for the excitations 
of the quarks, 



mi 



m„i 

*I3 



rrifi 



560.4622 GeV and m.i 



^4 



"^4 



m.i 



586.2638 GeV and m^i 



573.3793 GeV, 
558.9203 GeV, 
562.0523 GeV, 
560.2356 GeV, 
560.2514 GeV, 



for those of the leptons 



(B.l) 



m 1- 



m 1- 



mi- 



m 1- 



514.9604 GeV, 
505.4502 GeV, 



and for those of the gauge bosons 
m^i = 603.3141 GeV, 

m^i = 535.4923 GeV, m^^ii = 500.8931 GeV and m^i = 500.8931 GeV. 



(B.2) 



(B.3) 
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MUED processes with gauge boson excitations 



Process 


MG-FR 


CH-FR 


CH-ST 


SH-FR 


^1 ^1 ^ w- W+ 


2.856e+l 


2.854e+l 


2.855e+l 


2.855e+l 


^i+^i- ^ zz 


8.400e+0 


8.408e+0 


8.408e+0 


8.408e+0 




5.077e+0 


5.074e+0 


5.074e+0 


5.074e+0 




8.707e+0 


8.714e+0 


8.714e+0 


8.713e+0 




7.648e-l 


7.656e-l 


7.656e-l 


7.657e-l 


7I7I ^ it 


8.988e-2 


8.985e-2 


8.985e-2 


8.985e-2 


Z^ Z^ dd 


3.553e-l 


3.556e-l 


3.556e-l 


3.557e-l 


W^- e- e+ 


1.560e-l 


1.557e-l 


1.557e-l 


1.557e-l 


GG 


7.890e+l 


- 


7.858e+l 


7.854e+l 


G^ Z^ -^cc 


6.890e-l 


6.881e-l 


6.882e-l 


6.882e-l 




1.096e-l 


1.098e-l 


1.098e-l 


1.098e-l 


G^ G^ G^ G^ 


1.459e+5 


1.462e+5 


1.462e+5 


1.462e+5 


G^G^G^G 


9.534e+4 


9.542e+4 


9.539e+4 


9.535e+4 


Z7 ^ w^- 


3.185e+0 


3.185e+0 


3.185e+0 


3.185e+0 


W+W^~ ^ ZI7 


1.279e+2 


1.278e+2 


1.278e+2 


1.280e+2 


zz^ ^ w'^w^- 


4.944e+2 


4.944e+2 


4.944e+2 


4.952e+2 


w'^w^- z z^ 


4.930e+2 


4.926e+2 


4.926e+2 


4.932e+2 




2.199e+3 


2.202e+3 


2.202e+3 


2.203e+3 


2-1 Z^ W^+W^- 


1.905e+3 


1.905e+3 


1.905e+3 


1.904e+3 



Table 24: Cross sections for a selection of processes in MUED with two gauge bosons in the initial 
state. The existing MUED implementation in CalcHep is denoted CH-ST, while the FeynRules- 
generated ones in MadGraph CalcHep and Sherpa are denoted MG-FR, CH-FR and SH-FR. 
The center-of-mass energy is fixed to 1400 GeV, and a pt cut of 20 GeV is applied to each final 
state particle. 



Some examples of the obtained results for the calculation of cross sections of some 2 — > 2 
processes relative to the production of two Standard Model particles or two Kaluza-Klein 
excitations are shown in Tables 24 and 25. We set the center-of-mass energy to 1400 GeV 
and we applied a pr of 20 GeV on each final state particle. 
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MUED processes with fermion excitations 



Process 


MG-FR 


CH-FR 


CH-ST 


SH-FR 




1.109C-1 


1.109e-l 


1.109e-l 


l.llOe-1 


^ s s 


1.071e+0 


1.071e+0 


1.071e+0 


1.071e+0 


elr u]~^ — > e~ 
s ^ 


4.766e-l 


4.768e-l 


4.768e-l 


4.768e-l 


elr e]^ — > ^ ^ 
s S 1 ' 


2.078e-l 


2.079e-l 


2.079e-l 


2.079e-l 


Z^g Z^P — > u u 


1.635e-l 


1.635e-l 


1.635e-l 


1.636e-l 


u} u} W+ W~ 

e e 


5.905e-l 


5.901e-l 


5.901e-l 


5.897e-l 


uu 


2.298e-l 


2.298e-l 


2.298e-l 


2.299e-l 


el- el+ ^ e- e+ 

D D 


2.496e-l 


2.498e-l 


2.498e-l 


2.498e-l 


VgX — ^ du 


6.399e-l 


6.400e-l 


6.400e-l 


6.403e-l 


elr — > e~ Vp 

Ij *^ e ^ c 


1.052e+0 


1.052e+0 


1.052e+0 


1.052e+0 




l.llOe-1 


1.109e-l 


1.109e-l 


l.llOe-1 


rl- rl^ ^ r- r+ 

s s 


2.553e-l 


2.554e-l 


2.554e-l 


2.553e-l 




6.585e-l 


6.582e-l 


6.582e-l 


6.582e-l 


el- ^ e- r+ 


4.765e-l 


4.768c- 1 


4.768e-l 


4.768e-l 


1 1 ^ 1 1 

T T 


1.502e-l 


1.504e-l 


1.504e-l 


1.504C-1 


u}, z/,1 ^ c c 


1.634e-l 


1.635e-l 


1.635e-l 


1.636e-l 


ulul Z Z 

T T 


4.141e-l 


4.135e-l 


4.135e-l 


4.133e-l 


t}^t]+ -^bb 

D JJ 


1.426e-l 


1.427e-l 


1.427e-l 


1.428e-l 


tIt vl bt 

U T 


6.557e-l 


6.560e-l 


6.560e-l 


6.563e-l 


JJ T ^ 


1.053e+0 


1.052e+0 


1.052e+0 


1.052e+0 


7I r^.- - T- 


1.G380-1 


l.G39('-l 


i.G39('-l 


1.G390-1 


Zi ^ Gn 


2.065e+0 


2.068e+0 


2.068e+0 


2.072e+0 


d}^ d\) — > dd 


9.147e+0 


9.136e+0 


9.139e+0 


9.143e+0 


D D 


7.976e+0 


7.986e+0 


7.989e+0 


7.992C+0 


u\ u\ — > uu 


7.149e+0 


7.147e+0 


7.150e+0 


7.152e+0 


d\dh: ^ dd 




9.098e+0 


9.100e+0 


9.103e+0 


9.107e+0 


fl fl -^if 


8.471e+0 


8.474C+0 


8.477C+0 


8.474e+0 


An -^ff 


7.502e+0 


7.492C+0 


7.495C+0 


7.492C+0 


^ s ^ s ^ ^ ^ 


7.672e+0 


7.674e+0 


7.677e+0 


7.674e+0 


bl bl bb 


5.853e+0 


5.858e+0 


5.860e+0 


5.862e+0 


^ s ^ s ^ ^ ^ 


9.208e+0 


9.211e+0 


9.215e+0 


9.211e+0 


6^6^ ^ bb 


9.092e+0 


9.100e+0 


9.104e+0 


9.107e+0 



Table 25: Cross sections for a selection of fcrmionic processes in MUED. The existing MUED im- 
plementation in CalcHep is denoted CH-ST, while the FEYNRuLES-generated ones in MadGraph 
CalcHep and Sherpa are denoted MG-FR, CH-FR and SH-FR. The center-of-mass energy is fixed 
to 1400 GeV, and a px cut of 20 GeV is applied to each final state particle. 
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